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

lambeef

具体代码如下：
: ~; |6 D- [/ r! J( N+ l. Ll := 1000*10^(-9);
5 T. z7 G2 K+ k$ U2 m1 N6 M. dw := 50*10^(-9);
: H' [* G5 V1 I1 j" \7 Rk := (1/12)*w^3;
E := 150*10^9;
2 B4 R9 ?" U+ B$ Z) j: Jd := 20*10^(-9);
, G1 y* ]4 _# h  [2 Y/ T% N- s v := 2.3;
6 ^& R6 d: z/ v- Gelectrostatic := 8.854*10^(-12)*v^2*(1/2);
+ ~6 L1 V( ]' F+ x5 A+ y' Sode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
* Q+ J0 u) W2 z0 @bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
6 s4 Q6 u( E9 U+ L1 j% g$ C! cansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
6 @% c2 ^8 Z* ?5 a' Q- V9 v2 ^4 C4 n1 f问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！

eigen

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