-> alpha=0.001184048, beta=0.014709728
 alpha  =

   0.001184

 beta  =

   0.0147097


--> n=1259; m=5;

--> t=linspace(0,1,n+1)';

--> h=diff(t(1:2));

--> s=ones(n+1, m);

--> for j=1:n
  > dw=rand(1,m)*sqrt(h);
  > s(j+1,:)=s(j,:)+alpha*s(j,:)*h+beta*s(j,:).*dw;
  > end