function [err,stat] = WeakEulerBS(S0, sigma, K, T, N, M)
    truePrice = BS_call(S0,K,T,0,sigma);
    function y = Vdr(x)
        y = zeros(x);
    endfunction
    function y = Vdi(x)
        y = sigma*x;
    endfunction
    
    mu = 0;
    mu2 = 0;
    
    for i=1:M
        [Xbar, dB] = EulerTrajectory(S0, Vdr, Vdi, 1, N, T);
        mu = mu + max(Xbar-K,0);
        mu2 = mu2 + max(Xbar-K,0)^2;
    end
    err = abs(mu/M - truePrice);
    sd = sqrt(mu2/M-(mu/M)^2);
    stat = 1.65*sd/sqrt(M);
endfunction