求助：关于‘newton iteration not converging’问题

lambeef

具体代码如下：
l := 1000*10^(-9);
  v. Z5 n6 Q3 {w := 50*10^(-9);
3 d4 A' t3 y8 c/ t6 F& W9 P; H: D! Pk := (1/12)*w^3;
2 b+ \7 A5 L4 [. }  ^; m% xE := 150*10^9;
1 L9 q3 }0 ~  |1 S+ t1 }+ k3 @9 Md := 20*10^(-9);
& I6 A* V( M: \& M6 N* s: G v := 2.3;
electrostatic := 8.854*10^(-12)*v^2*(1/2);
* O2 N( l) u0 vode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
- o& ^7 L; ^" V# \  Ibcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
! f5 ~& \) N: e- D9 v) T- R6 @! D$ ~ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！

eigen

应该是迭代的初值没有选取好，或者这个问题本身不稳定。建议用Matlab软件再算一算。