function y=rexp(b,x) a=b(1); s=b(2); c1=b(3);