function ggprime=bvp1(t,gg,opt,a)

ggprime=[gg(2)
    gg(3)
    -a*gg(2)-sin(gg(1))+cos(t)*sin(t)];