function y=Corr_CV_MC_call(S0,K,T,r,sigma,M)
    X=rand(1,M,'normal');
    S=S0*exp(sigma*X*sqrt(T)+(r-1/2*sigma^2)*T);
    C=exp(-r*T)*max(S-K,0);
    m=1/M*sum(S);
    P=1/M*sum(C);
    VarEst1=1/(M-1)*sum((S-m)^2);
    VarEst2=1/(M-1)*sum((C-P)^2);
    CoVarEst=1/(M-1)*sum((S-m).*(C-P));
    CorrEst=CoVarEst/sqrt(VarEst1*VarEst2);
    y=CorrEst;Co
endfunction