我是mathematica的新手。请帮我解决这个问题。
形式
NDSolve[ode==0,bc1==0,bc2==0,f,{R,1,2}]
被修改
错误：Power :: infy：“无限表达式1/0。^（11/6）遇到” NDSolve :: ndnum：在R == 1
处遇到导数的非数值
BC1 = -1 + (0.8*e5*(-1 - f[R]^2/R^2))/((f[R]^2*Derivative[1][f][R]^2)/R^2)^(5/6)
+ 2*e5*(-1 + Sqrt[(f[R]^2*Derivative[1][f][R]^2)/R^2]) /.R -> 2
BC2 = 10 - (0.8*e5*(R^2 + f[R]^2))/(R^2*((f[R]^2*Derivative[1][f][R]^2)/R^2)^(5/6))
+ 2*e5*(-1 + Sqrt[(f[R]^2*Derivative[1][f][R]^2)/R^2]) /. R -> 1
NDSolve[{
-((f[R]*((0.8*e5*(-1 - f[R]^2/R^2))/((f[R]^2*Derivative[1][f][R]^2)/R^2)^(5/6) + 2*e5*(-1 + Sqrt[(f[R]^2*Derivative[1][f][R]^2)/R^2])))/R^2) +
(Derivative[1][f][R]*((0.8*e5*(-1 - f[R]^2/R^2))/((f[R]^2*Derivative[1][f][R]^2)/R^2)^(5/6) + 2*e5*(-1 + Sqrt[(f[R]^2*Derivative[1][f][R]^2)/R^2])))/R +
((f[R]*((0.8*e5*(-1 - f[R]^2/R^2))/((f[R]^2*Derivative[1][f][R]^2)/R^2)^(5/6) + 2*e5*(-1 + Sqrt[(f[R]^2*Derivative[1][f][R]^2)/R^2])))/R -
Derivative[1][f][R]*((0.8*e5*(-1 - Derivative[1][f][R]^2))/((f[R]^2*Derivative[1][f][R]^2)/R^2)^(5/6) + 2*e5*(-1 + Sqrt[(f[R]^2*Derivative[1][f][R]^2)/R^2])))/R +
(1/R)*(f[R]*((0.8*e5*((2*f[R]^2)/R^3 - (2*f[R]*Derivative[1][f][R])/R^2))/((f[R]^2*Derivative[1][f][R]^2)/R^2)^(5/6) -
(0.6666666666666667*e5*(-1 - f[R]^2/R^2)*(-((2*f[R]^2*Derivative[1][f][R]^2)/R^3) + (2*f[R]*Derivative[1][f][R]^3)/R^2 + (2*f[R]^2*Derivative[1][f][R]*Derivative[2][f][R])/R^2))/
((f[R]^2*Derivative[1][f][R]^2)/R^2)^(11/6) + (e5*(-((2*f[R]^2*Derivative[1][f][R]^2)/R^3) + (2*f[R]*Derivative[1][f][R]^3)/R^2 + (2*f[R]^2*Derivative[1][f][R]*Derivative[2][f][R])/R^2))/
Sqrt[(f[R]^2*Derivative[1][f][R]^2)/R^2])) == 0,
BC1 == 0,BC2 ==0
},
f, {R, 1, 2}]