edongliang
2010-06-01, 21:50
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的解下面就不动了，像死机了一样的！！
这是什么原因？
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的解下面就不动了，像死机了一样的！！
这是什么原因？