matlab解方程组[c1c2r]=solve('13.622=c1^0.5+0.95*c2^0.5','c2=105-(1+r)*c1','r+1=(1/0.95^2)*(13.622*c1^(-0.5)-1)')为什么解出来明显不对.跪谢

Question

matlab解方程组
[c1 c2 r]=solve('13.622=c1^0.5+0.95*c2^0.5','c2=105-(1+r)*c1','r+1=(1/0.95^2)*(13.622*c1^(-0.5)-1)')为什么解出来明显不对.跪谢

Answer 1

s=solve('13.622=c1^0.5+0.95*c2^0.5','c2=105-(1+r)*c1','r+1=(1/0.95^2)*(13.622*c1^(-0.5)-1)')
c1=s.c1
c2=s.c2
r=s.r
用maple作符号运算的版本能解.
如用5.3版、7.1版均能解出如下结果：
c1 = 44.427855836185434268940742282114
c2 = 53.622129458377212486360933276572
r = .15643371877012194670660012187269