function rhs = bvp_shooting_rhs(t,y) % y'' = -sin(y) % % Change variables: % y1 = y % y2 = y' % % Equation becomes 1st-order system: % y1' = y2 % y2' = -sin(y1) rhs = [y(2) -sin(y(1))];