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

lambeef

具体代码如下：
; X. x# D" B* y+ E: s2 o! y1 ml := 1000*10^(-9);
$ e! i. L; c1 J+ s; b' q' iw := 50*10^(-9);
; U1 v0 l( t2 s& [  Wk := (1/12)*w^3;
% P- S2 Z: L* Z" x. GE := 150*10^9;
d := 20*10^(-9);
$ H' e" W; }0 f4 C: J v := 2.3;
2 m4 r1 d2 V4 C. B3 m( y) e* M% ielectrostatic := 8.854*10^(-12)*v^2*(1/2);
+ E9 `5 w7 g6 q1 @ode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
3 x9 {- Z* t4 F9 H# ^1 i/ @, q( n问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！

eigen

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