function f=maxwell(a,x) f=a(1)*(x.^.5).*exp(-x*a(2)); %f=a(1)*(x.^2).*exp(-x.^2.*a(2));