julia
Continuation of 10.1.5. Now compute \(u'(0) - u'(1)\) for the three values of \(\epsilon\).
for epsilon in [0.2, 0.02, 0.002] x, u, du_dx = FNC.shoot(phi, (a,b), g1, g2, [-1.0, 0.0]); @show du_dx[1] - du_dx[end] end