% AsianMC.m
function [P,CI] = AsianMC(S0,X,r,T,sigma,NSamples,NRepl)
Payoff = zeros(NRepl,1);
for i=1:NRepl
   Path=AssetPaths1(S0,r,sigma,T,NSamples,1);
   Payoff(i) = max(0, mean(Path(2:(NSamples+1))) - X);
end
[P,aux,CI] = normfit( exp(-r*T) * Payoff);