BVP，NDSolve：与方程相比，因变量更少，因此系统超定

Question

我是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}]