function W = BM(T, nSteps, nTrajectories)

    dt = T/nSteps;
    dW = sqrt(dt)*randn(nTrajectories, nSteps);
    W = [zeros(nTrajectories, 1), cumsum(dW, 2)];

end
