[u,l,xx]=gausstype2(x,P);

%% gausstype2.m

%% Computes upper and lower membership functions for a Gaussian
%% primary membership functions with uncertain means.

function [u,l,xx]=gausstype2(x,P);

sigma=P(:,1);
m1=min([P(:,2),P(:,3)]);
m2=max([P(:,3),P(:,2)]);

L=length(x);

u=[];
l=[];

for i=1:L
if x(i)>=m1 & x(i)<=(m1+m2)/2
mu1=1; 
mu2=gaussmf(x(i),[sigma,m2]);

elseif  x(i)>(m1+m2)/2 & x(i)<=m2
mu1=1;
mu2=gaussmf(x(i),[sigma,m1]);

else
mu1=gaussmf(x(i),[sigma,m1]);
mu2=gaussmf(x(i),[sigma,m2]);
end

u=[u,max([mu1,mu2])];
l=[l,min([mu1,mu2])];
end

xx=x;