function [z,s]=simbm(dim,t)
% not completed
%simulates a standard brownian motion on the interval t(1) t(2) of dimension dim
%[z,s]=simbm(2,[0 1]);
z=[];
delta=(t(2)-t(1))/10000;
s=t(1):delta:t(2);
for i=1:dim
z=[z;[0 cumsum(normrnd(0,sqrt(delta),1,10000))]];
end