#include  "trs5-inc.mm"

Matrix numerical_jacobi();

Func Matrix Dftrs5(t, x, y)
     // Dftrs5 - 微分方程式の微分値と，変分系の微分値
     Real t;
     Matrix x, y;
{

//       Real epsi;
//       Integer dim, k, l;
//       Matrix epsm, jacobian;

//       epsi= 0.0001;
//       dim= length(x);
//       epsm = epsi*I(dim);

//       jacobian = Z(dim);
//       for (k=1; k<=dim; k++){
//  	  for (l=1; l<=dim; l++){
//  	       jacobian(k,l)=
//  		    [ftrs5(t,x+epsm(:,l))-ftrs5(t,x)](k) /epsi;
//  	  }
//       }
//       return jacobian*y;

     return numerical_jacobi(ftrs5, t, x, xbase*1.0e-4)*y;
}