M = 7;
T = 2;
N = 500; 
h = T/N;
t = (0:N)*h;
Z = randn(M,N); 

W = zeros(M,N+1);
for i=1:N
    W(:,i+1)=W(:,i)+sqrt(h)*Z(:,i);
end

% monochrome plot
% hold on
% for i=1:M
%    plot(t,W(i,:))
% end

% color plot - for up to 7 colors only! 
hold on
plot(t,W(1,:),'b')
plot(t,W(2,:),'r')
plot(t,W(3,:),'g')
plot(t,W(4,:),'y')
plot(t,W(5,:),'k')
plot(t,W(6,:),'m')
plot(t,W(7,:),'c')