它可以运行的. 

如果要算个偏导数: d(Prob)/d(spot) 如何弄,好像很麻烦啊,除非我把表达式直接都写出来,然后用符号运算来求偏导的表达式?  另外那里可以下载到一个小一点的matlab软件? 700M的没法下载啊

    T=1/4;
    t=0;
    miu=0.004505553;
    sigma=0.21;
    spot=1.3657;
    l=log((spot+0.1)/(spot-0.1));
    x=log(spot/(spot-0.1));
    k=1:100;
    k=k';
    lambda=(miu^2/(sigma^2)+(pi^2)*(sigma^2)*(k.^2)/(l^2))/2;
    s=pi*sum((((exp(-lambda)*(T-t))./lambda).*k).*sin(k*pi*(l-x)/l));
    s0=(i*sin(i*miu*x/sigma^2))/(i*sin(i*miu*l/sigma^2));
    prob=exp(miu*(l-x)/sigma^2)*(s0-sigma^2*s/l^2);

[此贴子已经被作者于2009-1-9 21:01:58编辑过]