%
% M = number of sample paths
%
function [S] = SimPath(M, N, T, spot, mu, sigma)


    % model parameters
    
    dt  = T/N;
    
    z = randn(M,N);
    S = spot * ones(M,N+1);
    
    for i = 2 : N +1
        S(:,i) = S(:,i-1) .* exp( (mu-0.5*sigma^2)*dt + sqrt(dt)*sigma*z(:,i-1) );
    end
    
end