Math-Lsoda
view release on metacpan or search on metacpan
int liw = 20 + dim;
int itol = 2;
int itask = 1;
int istate = 1;
int iopt = 0;
int jt = 2;
double a_x[dim], a_rtol[dim], a_atol[dim];
double rwork[lrw];
int iwork[liw];
FILE *fp = PerlIO_findFILE(stream);
int maxStep = required_size(t, tout, dt);
double t1 = t + dt;
for(i=0;i<dim;i++){
SV** pv = av_fetch(x,i,0);
a_x[i] = SvNV(*pv);
pv = av_fetch(rtol,i,0);
a_rtol[i] = SvNV(*pv);
pv= av_fetch(atol,i,0);
a_atol[i] = SvNV(*pv);
}
( run in 0.343 second using v1.01-cache-2.11-cpan-0d8aa00de5b )