wanghao023
2018-12-01, 21:11
原来的代码
syms d c1 c2 v1 v2 p;
f=d/4*(1+3*(1-(p/2)^(1/3)))-v2;
g=c1*d*v2-p;
h=1.483+25.306*log10(1+v2/5.19)-d;
m=417*(c2^6.29-c1^6.29)-p+2;
l=0.25*((p-2)/(1/c1-1/c2))-d^2;
o=417*10^6*(c2^6.29-c1^6.29)/c1*(1-c1^2/c2^2)-v2;
t=299*10^6*((c2^2/c1)^7.15-1)-p;
[d,c1,c2,v1,v2,p]=solve(g,h,m,l,o,t)
请问各位是我哪里写错了，电脑计算了十来分钟也没个结果。
感谢各位了。
syms d c1 c2 v1 v2 p;
f=d/4*(1+3*(1-(p/2)^(1/3)))-v2;
g=c1*d*v2-p;
h=1.483+25.306*log10(1+v2/5.19)-d;
m=417*(c2^6.29-c1^6.29)-p+2;
l=0.25*((p-2)/(1/c1-1/c2))-d^2;
o=417*10^6*(c2^6.29-c1^6.29)/c1*(1-c1^2/c2^2)-v2;
t=299*10^6*((c2^2/c1)^7.15-1)-p;
[d,c1,c2,v1,v2,p]=solve(g,h,m,l,o,t)
请问各位是我哪里写错了，电脑计算了十来分钟也没个结果。
感谢各位了。