Hello, I’m new to the BoundaryValueDiffEq package and I’ve seen some examples of its use for solving second order differential equations. However, I was wondering whether it’s posible to solve first order BVPs too.
I tried to initialize the following system but the solver throws retcode: failure
function ode!(du,u,p,t) du = sqrt((1 + p - cos(u)) * (1 + cos(u / p))) end function bc2!(residual, u, p, t) residual = u +p*pi residual = u[end]-p*pi end tspan=(-50,50) bvp2 = TwoPointBVProblem(ode!, bc2!,[-4*pi,4*pi] ,tspan,[0.1,4]) sol2 = solve(bvp2, MIRK4(), dt=0.5)