Hello Sheehan,
I get a failure with the following example, do you have an idea about the
why ?:
*# solves u" = phi²*sinh(u)-2u'/(x+gamma) , u'(a)=-xi, u'(R)=0*
*a= 3.514457e-07*
*R= 7.60773e-07*
*x=Fun(identity, Interval(a,R))*
*d=x.domain*
*B=neumann(d)*
*D=diff(d)*
*# Solves Lu + g(u) == 0*
*phi=1.341211*
*gamma=0.8585931*
*L = D^2 + 2/(x.+gamma)*D*
*g = u -> -phi^2*(exp(u)-exp(-u))/2; gp = u -> -phi^2*(exp(u)+exp(-u))/2*
*u=0.x #initial guess *
*xi=9.403218*
*for k=1:5*
* u=u-[B,L+gp(u)]\[diff(u)[a]+xi,diff(u)[R],L*u+g(u)];*
*end*
*julia> u*
*IFun{Float64,Interval{Float64}}([NaN,NaN,NaN,NaN,NaN,NaN,NaN,NaN,NaN,NaN
…
NaN,NaN,NaN,NaN,NaN,NaN,NaN,NaN,NaN,NaN],Interval{Float64}(3.514457e-7,7.60773e-7))*
