新装的gauss8,运行别人书中给出的程序,但是不知道为什么不显示出结果,下面附上结果,求能人解答或者告诉我怎么运行能出结果,谢谢大家
proc steady(p);
local alp, bet, del, rh, sige, ph, ps, _1_alp, ar, th, et, ybar, cbar, ibar, nbar1, nbar,
fvc, jc, tcode, lbar, kbar, abar;
alp = p[1];
bet = p[2];
del = p[3];
rh = p[4];
sige= p[5];
ph = p[6];
ps = p[7];
_1_alp = 1-alp;
ar = 1/bet - 1;
th = (alp/(ar+del))^(1/_1_alp);
et = th^alp;
nbar1 = 1 + (1/_1_alp)*((1-ps)/ps)*(1-del*(th^_1_alp));
nbar = 1/nbar1;
lbar = 1-nbar;
cbar = (et - del*th)*nbar;
kbar = th*nbar;
ybar = et*nbar;
ibar = ybar - cbar;
abar = 1;
retp(ybar|cbar|ibar|nbar|lbar|kbar|abar);
endp;


雷达卡




京公网安备 11010802022788号







