#include #include "nr.h" #define FREERETURN {free_matrix(fjac);free_vector(fvec);\ free_vector(p);free_ivector(indx);return;} void mnewt(int ntrial, double x[], int n, double tolx, double tolf) { int k,i,*indx; double errx,errf,d,*fvec,**fjac,*p; indx=ivector(n); p=vector(n); fvec=vector(n); fjac=matrix(n,n); for (k=0;k