#ifndef dfridr_h
#define dfridr_h

#include <default.h>
#include <matrix.h>

static double maxarg1,maxarg2;
#define DMAX(a,b) (maxarg1=(a),maxarg2=(b),(maxarg1) > (maxarg2) ? (maxarg1) : (maxarg2))

double dfridr(double (*func)(double), double x, double h, double *err);

#endif