function [J] = sagas(t,problem)
[J] = mga_dsm(t,problem);