hold off;f=5.11;e=1;
 m=400.3/e;
 a=m*f/(f-1);b=m-a;
 sigma=34.22*(e^(-.5));
 lamda=m*f/((f-1)*sigma);
 d=lamda^2;
x=-200:.25:1000;k=max(size(x));n=1;
for i=1:k,u(i)=1+(x(i)-m)/a;
  if(u(i)>0)
    f1(i)=(2*pi/n)^(-.5)*sigma^(-1)*u(i)^(-1.5)*exp(-n*((x(i)-m)^2)/(2*sigma^2*u(i)));
  else
    f1(i)=0;
  end;
end;
t=14.68*n^(-.5);
g=(2*pi)^(-.5)*t^(-1)*exp(-x.^2./(2*t^2));
t1=(t^2+sigma^2)^.5
g1=(2*pi)^(-.5)*t1^(-1)*exp(-(x-m).^2./(2*t1^2));
semilogy(x,f1,'y'); hold on;
s=conv(f1,g);
%for j=2390:1:2410,v(j)=f1(j)-(.0625)*s(j+800)
%end
sum=0;k=1
while(k<=length(s))
      sum=sum+s(k);
     k=k+1;
end
sum1=0;
for j=1:2374;
sum1=sum1+s(j);
end
semilogy(x,g1,'g');

t2=0:length(s)-1;
p=.25
t3=(t2+2*x(1)/p)*p;
semilogy(t3,s/16,'b'),

axis([40,400,10^(-12.5),10^(0)]);


e=10;f=5.11;
 m=400.3/e;
 a=m*f/(f-1);b=m-a;
 sigma=34.22*(e^(-.5));
 lamda=m*f/((f-1)*sigma);
 d=lamda^2;
x=-200:.25:1000;k=max(size(x));n=1;
for i=1:k,u(i)=1+(x(i)-m)/a;
  if(u(i)>0)
    f10(i)=(2*pi/n)^(-.5)*sigma^(-1)*u(i)^(-1.5)*exp(-n*((x(i)-m)^2)/(2*sigma^2*u(i)));
  else
    f10(i)=0;
  end;
end;
t=14.68*n^(-.5);
g=(2*pi)^(-.5)*t^(-1)*exp(-x.^2./(2*t^2));
t0=(t^2+sigma^2)^.5
gi=(2*pi)^(-.5)*t0^(-1)*exp(-(x-m).^2./(2*t0^2));
semilogy(x,f10,'y'); hold on;
s0=conv(f10,g);
%for j=2390:1:2410,v(j)=f10(j)-(.0625)*s(j+800)
%end
sum=0;k=1
while(k<=length(s0))
      sum=sum+s0(k);
     k=k+1;
end
semilogy(x,gi,'g'),
sum0=0;
for i=2374:length(s0);
sum0=sum0+s0(i);
end
t2=0:length(s0)-1;
p=.25
t3=(t2+2*x(1)/p)*p;
semilogy(t3,s0/16,'r'),
axis([40,400,10^(-12.5),10^(0)]);
p=(sum0+sum1)/32
[xt,yt]=ginput;
a=xt
