希望xuehe老师帮忙，非常感谢！

xuehe似乎不在，如下这样可以么？

new;
proc (1)=fun(par);
local VA,sigmaA,d1,d2,f;
local X,VE,r,sigmaE,T;
X=10;VE=3;r=0.05;T=1;sigmaE=0.4;
VA=par[1];sigmaA=par[2];
d1=(ln(VA/X)+(r+sigmaA^2/2)*T)/(sigmaA*sqrt(T));
d2=d1-sigmaA*sqrt(T);
f=zeros(2,1);
f[1]=VA*pdfn(d1)+exp(-r*T)*X*pdfn(d2)-VE;
f[2]=VA/VE*pdfn(d1)*sigmaA-sigmaE;
retp(f);
endp;

par0={38 ,0.6};
eqsolve(&fun,par0);


||F(X)|| at final solution:                                          19.484819
--------------------------------------------------------------------------------
Termination Code = 1:
Norm of the scaled function value is less than __Tol;
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
VARIABLE              START                   ROOTS                 F(ROOTS)
--------------------------------------------------------------------------------
X1                  38.00000                  38.961426           1.463655e-006
X2                   0.60000                 0.80000001          1.9773507e-007
--------------------------------------------------------------------------------

       38.961426
      0.80000001
       1.0000000

iooo老师，就是这样，非常非常感谢！

iooo老师很不错，上次帮台湾人解决问题，这次有居功至伟。
Matlab有很多工具箱，这是它的长处，但因为挂的太多，启动慢。GAUSS轻便灵活，这是它的长处。一般来说，matlab写的，gauss搞熟了，都可以改写。我们这里是中文GAUSS很不错的平台。

非常同意xuehe老师的指导，非常非常感谢！