Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
初级会员
注册日期: 2010-06-01
年龄: 37
帖子: 6
声望力: 0 ![]() |
![]()
for ii=0:0.05:1
x=vpa(double(solve(9*u^2/((1/3*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)+1/3*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)*(1/3*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2)-1/3*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)))^(1/2)/(1/4*3^(1/2)+2*((1/2+1/6*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)+1/12*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(-1/2+1/6*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)+1/12*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(1/2-1/6*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)+1/12*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(1/2+1/6*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)-1/12*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2)))^(1/2)+1/2*((1/2+1/6*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(-1/2+1/6*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2)))^(1/2))-ii)),8) // x=vpa(x,8)表示对X小数点后保留8位 end x 上面的红色部分是matlab解出来的 只有x=0的解下面就不动了,像死机了一样的!! 这是什么原因? |
![]() |
![]() |