SharedModel;
Variable x,y(2);//y1,y2;
Parameter a,b;
Constant epsilon=0.3921,h0=6.3,R=125,Sf=0.03143,VR=1308.33,h1=3.82,Re=133,z0=134.91;
ConstStr z=z0-sqrt(Re^2-x^2),Vx=h1*VR*(1+Sf)/(2*z),delta_f=15+9.2*epsilon^0.11,p1=exp(2*a*x/z),
           p2=(Vx-VR*(1+Sf))*2*b*x/(z^2)-2*delta_f*sqre(Re^2-x^2)/z ,p3=-(Vx-VR*(1+Sf))*2*a*b*x^2/(z^3),
           p4=2*delta_f*a*x*sqrt(Re^2-x^2)/(z^2)-2*delta_f*a*Re^2*arcsin(x/Re)/(z^2),p5=(1+4*Re/h1)*delta_f;
Function y1=p1*(p2+p3+p4+p5);
             y2=2*a*y1+delta_f*sqrt(Re^2-x^2)/z+(Vx-VR*(1+Sf))*2*b/z;
Data;
//x,y1,y2
1.12891986062718e+000    6.40830449826990e+001   3.19723183391004e+001
2.25783972125436e+000    1.00069204152249e+002   3.47404844290657e+001
3.76306620209059e+000    8.56747404844291e+001   4.13840830449827e+001
5.01742160278746e+000    8.06920415224914e+001   2.58823529411765e+001
6.39721254355401e+000    7.62629757785467e+001   2.25605536332180e+001
7.77700348432056e+000    6.68512110726644e+001   1.75778546712803e+001
8.96864111498258e+000    5.96539792387543e+001   9.68858131487885e-001
1.02229965156794e+001    5.57785467128028e+001  -1.61937716262976e+001
1.16655052264808e+001    5.35640138408304e+001  -2.83737024221453e+001
1.28571428571429e+001    5.35640138408304e+001  -3.05882352941176e+001
1.42369337979094e+001    5.46712802768166e+001  -3.44636678200692e+001
1.54285714285714e+001    4.85813148788927e+001  -1.95155709342561e+001
1.67456445993031e+001    4.30449826989619e+001  -1.01038062283737e+001
1.75609756097561e+001    4.35986159169550e+001  -4.01384083044982e+000

運(yùn)行提示：在程序所有行的左邊沒有發(fā)現(xiàn)y(2)�。。。�

我上面function不是有y2嗎，為什么還提示這個(gè)？
y2=2*a*y1+delta_f*sqrt(Re^2-x^2)/z+(Vx-VR*(1+Sf))*2*b/z;
像這個(gè)y1作為y2輸入的情況能否這樣寫代碼？？？
跪求指導(dǎo)�。。�！ 返回小木蟲查看更多