function [J] = cassini1(t,problem)
[J] = mga(t,problem);
