附件中的联立方程如何用gauss编程，求助，在附件中，非常感谢！

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

附件中的联立方程如何用gauss编程，求助，在附件中，非常感谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：GAUSS 联立方程 非常感谢 如何用 Aus 求助 方程 附件 感谢 GAUSS

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

zhangtao 发表于 2011-3-15 19:20
附件中的联立方程如何用gauss编程，求助，在附件中，非常感谢！

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老师的指导，非常非常感谢！