hold off;f=5.11;e=2;
 m=400.3/e;
 a=m*f/(f-1);b=m-a;
 sigma=102.66*(e^(-.5));
 lamda=m*f/((f-1)*sigma);
 d=lamda^2;
x=-200:.125:500;k=max(size(x));n=49000;
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=5*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=3190:3210,v(j)=f1(j)-(.0625/4)*s(j+1600)
end
sum=0;k=1
while(k<=length(s))
      sum=sum+s(k);
     k=k+1;
end
t2=0:length(s)-1;
p=.125
t3=(t2+2*x(1)/p)*p;
semilogy(t3,s/64,'r'),
%axis([130,707,10^(-295),10^(05)]);
semilogy(x,g,'g')