diff -Naur a/minpack.h b/minpack.h
--- a/minpack.h	1969-12-31 21:00:00.000000000 -0300
+++ b/minpack.h	2024-04-16 10:24:21.429666175 -0300
@@ -0,0 +1,113 @@
+/* Declarations for minpack */
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+/* find a zero of a system of N nonlinear functions in N variables by
+   a modification of the Powell hybrid method (Jacobian calculated by
+   a forward-difference approximation) */
+void hybrd1_ ( void (*fcn)(int *n, double *x, double *fvec, int *iflag ), 
+	       int *n, double *x, double *fvec, double *tol, int *info,
+	       double *wa, int *lwa );
+
+/* find a zero of a system of N nonlinear functions in N variables by
+   a modification of the Powell hybrid method (Jacobian calculated by
+   a forward-difference approximation, more general). */
+void hybrd_ ( void (*fcn)(int *n, double *x, double *fvec, int *iflag ), 
+	      int *n, double *x, double *fvec, double *xtol, int *maxfev,
+	      int *ml, int *mu, double *epsfcn, double *diag, int *mode,
+	      double *factor, int *nprint, int *info, int *nfev,
+	      double *fjac, int *ldfjac, double *r, int *lr, double *qtf,
+	      double *wa1, double *wa2, double *wa3, double *wa4);
+  
+/* find a zero of a system of N nonlinear functions in N variables by
+   a modification of the Powell hybrid method (user-supplied Jacobian) */
+void hybrj1_ ( void (*fcn)(int *n, double *x, double *fvec, double *fjec,
+			   int *ldfjac, int *iflag ), int *n, double *x,
+	       double *fvec, double *fjec, int *ldfjac, double *tol,
+	       int *info, double *wa, int *lwa );
+          
+/* find a zero of a system of N nonlinear functions in N variables by
+   a modification of the Powell hybrid method (user-supplied Jacobian,
+   more general) */
+void hybrj_ ( void (*fcn)(int *n, double *x, double *fvec, double *fjec,
+			  int *ldfjac, int *iflag ), int *n, double *x,
+	      double *fvec, double *fjec, int *ldfjac, double *xtol,
+	      int *maxfev, double *diag, int *mode, double *factor,
+	      int *nprint, int *info, int *nfev, int *njev, double *r,
+	      int *lr, double *qtf, double *wa1, double *wa2,
+	      double *wa3, double *wa4 );
+
+/* minimize the sum of the squares of nonlinear functions in N
+   variables by a modification of the Levenberg-Marquardt algorithm
+   (Jacobian calculated by a forward-difference approximation) */
+void lmdif1_ ( void (*fcn)(int *m, int *n, double *x, double *fvec,
+			   int *iflag ),
+	       int *m, int *n, double *x, double *fvec, double *tol,
+	       int *info, int *iwa, double *wa, int *lwa );
+
+/* minimize the sum of the squares of nonlinear functions in N
+   variables by a modification of the Levenberg-Marquardt algorithm
+   (Jacobian calculated by a forward-difference approximation, more
+   general) */
+void lmdif_ ( void (*fcn)(int *m, int *n, double *x, double *fvec,
+			  int *iflag ),
+	      int *m, int *n, double *x, double *fvec, double *ftol,
+	      double *xtol, double *gtol, int *maxfev, double *epsfcn,
+	      double *diag, int *mode, double *factor, int *nprint,
+	      int *info, int *nfev, double *fjac, int *ldfjac, int *ipvt,
+	      double *qtf, double *wa1, double *wa2, double *wa3,
+	      double *wa4 );
+
+/* minimize the sum of the squares of nonlinear functions in N
+   variables by a modification of the Levenberg-Marquardt algorithm
+   (user-supplied Jacobian) */
+void lmder1_ ( void (*fcn)(int *m, int *n, double *x, double *fvec,
+			   double *fjec, int *ldfjac, int *iflag ),
+	       int *m, int *n, double *x, double *fvec, double *fjec,
+	       int *ldfjac, double *tol, int *info, int *ipvt,
+	       double *wa, int *lwa );
+
+/* minimize the sum of the squares of nonlinear functions in N
+   variables by a modification of the Levenberg-Marquardt algorithm
+   (user-supplied Jacobian, more general) */
+void lmder_ ( void (*fcn)(int *m, int *n, double *x, double *fvec,
+			  double *fjec, int *ldfjac, int *iflag ),
+	      int *m, int *n, double *x, double *fvec, double *fjec,
+	      int *ldfjac, double *ftol, double *xtol, double *gtol,
+	      int *maxfev, double *diag, int *mode, double *factor,
+	      int *nprint, int *info, int *nfev, int *njev, int *ipvt,
+	      double *qtf, double *wa1, double *wa2, double *wa3,
+	      double *wa4 );
+
+/* minimize the sum of the squares of nonlinear functions in N
+   variables by a modification of the Levenberg-Marquardt algorithm
+   (user-supplied Jacobian, minimal storage) */
+void lmstr1_ ( void (*fcn)(int *m, int *n, double *x, double *fvec,
+			   double *fjrow, int *iflag ), int *m, int *n,
+	       double *x, double *fvec, double *fjac, int *ldfjac,
+	       double *tol, int *info, int *ipvt, double *wa, int *lwa );
+
+/* minimize the sum of the squares of nonlinear functions in N
+   variables by a modification of the Levenberg-Marquardt algorithm
+   (user-supplied Jacobian, minimal storage, more general) */
+void lmstr_ ( void (*fcn)(int *m, int *n, double *x, double *fvec,
+			  double *fjrow, int *iflag ), int *m,
+	      int *n, double *x, double *fvec, double *fjac,
+	      int *ldfjac, double *ftol, double *xtol, double *gtol,
+	      int *maxfev, double *diag, int *mode, double *factor,
+	      int *nprint, int *info, int *nfev, int *njev, int *ipvt,
+	      double *qtf, double *wa1, double *wa2, double *wa3,
+	      double *wa4 );
+ 
+void chkder_ ( int *m, int *n, double *x, double *fvec, double *fjec,
+	       int *ldfjac, double *xp, double *fvecp, int *mode,
+	       double *err  );
+
+double dpmpar_ ( int *i );
+
+double enorm_ ( int *n, double *x );
+
+#ifdef __cplusplus
+}
+#endif