void fill_dym_dxn(VectorType const & V, ScalarType * res) const
{
for(std::size_t m = 0 ; m < M_ ; ++m){
ScalarType t = (ScalarType)(m+1)/5;
double a = V[0] + t*V[1] - exp(t);
double b = V[2] + V[3]*sin(t) - cos(t);
get(res,m,0) = 2*a;
get(res,m,1) = 2*t*a;
get(res,m,2) = 2*b;
get(res,m,3) = 2*sin(t)*b;
}
}
