m=5;

dt= 0.0001;
n= 1/dt;

itoerr=0;

R= sqrt(dt)*randn(n,m);
E=zeros(1,m);
B=zeros(1,m);
for dB=R';
  E=E + B .* dB';
  B=B+dB';
end
Ito=0.5*B.^2-0.5*n*dt
Real=E



  