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