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));
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
t2=0:length(s)-1;
p=.25
t3=(t2+2*x(1)/p)*p;
semilogy(t3,s/16,'r'),
axis([40,400,10^(-12.5),10^(0)]);
semilogy(x,g,'g')