sols = NDSolve[{x''[t] + 3/t*x'[t] - t*x[t] + 10.1406 (x[t])^2 -
9.1406*(x[t])^3 == 0, x'[0] == 0, x[3] == 0}, x, t,
Method -> {"Shooting",
"StartingInitialConditions" -> {x[0] == 0.7, x'[0] == 0}}]
Plot[Evaluate[x[t] /. sols], {t, 0, 3}, PlotStyle -> { Black}]
各位大佬，怎么操作才可以不考虑t=0处的发散？同时不影响结果

……如果是可去间断点，那么把0点稍微错开那么一点是最简单的做法。这问题真正麻烦的是非线性边值问题的求解。稍微试了一下，似乎不好办。