[Fractint] EHNewt

An attempt at solving a second order non linear differential equation has produced this formula , using Euler Heun and Newton's method ; where : x method , with complex divisor %i*s , for a second order

non linear differential equation.

Select zmag for outer .

} EHNewt{; EH Newton , w1 = pixel, sn=-1, Pd = 1*(10^(-13)), a5= 42.469607686999 , b1= 23546.24999999999 , b2= 9.343749999999999*(10^(8)) , b3= 96.8 , s=0.01 p=0.39757,

h=0.71743,

x = 0.1, y = 0, tol = 1:

u = y w = ( sn*a5*y*y+sn*(b2*x + b3)*y +sn*( sn*b2*w1 + sn*b1)*x +sn*b1)/s w = sn*log(w)/a5

a = x + p*u b = y + p*w

c = x + h*(u + b) d = ( sn*a*( sn*b2*w1 +sn*b1 ) +sn*b*( b3 + a*b2 ) + sn*b1 + sn*a5*b*b)/s

d = y + h*(w + d)

d = sn*log(d)/a5

x = c y = d z = y

tol = |u-y| tol > Pd

}