#include <stdio.h>
#include <malloc.h>

#include "atom.h"
#include "nrutil.h" 
#include "integrate.h"

#define n 3
void jacobn(double z, double *y, double *dfdz, double **fjac, int nv)
{
 int choice = 1;
 double *dummy;
 dummy = (double *) malloc((n+1)*sizeof(double));
 evaluate(choice, z, y, dummy, dfdz, fjac, n);
 free(dummy);
}

void derivs(double z, double y[], double fvec[])
{
  int choice = 2;
  double *dfxx, **fxac;
  dfxx = (double *) malloc((n+1)*sizeof(double));
  fxac = dmatrix(1,n,1,n);
  evaluate(choice, z,  y, fvec, dfxx, fxac, n);
  free(dfxx);   free_dmatrix(fxac,1,n,1,n);

}