/* simeq_newton3.h  solve up to third order system of nonlinear equations */
void simeq_newton3(int n,        /* number of unknowns, linear variables */
                   int nlin,     /* number of nonlinear terms */
		   double A[n][nlin], /* [n][n+nlin] coefficients */
		   double Y[],  /* [n] right hand side */
		   int var1[],  /* [n+nlin] variables first, then nonlinear */
		   int var2[],  /* [n+nlin] -1 for variable, else variable  */
		   int var3[],  /* [n+nlin] -1 for variable, else variable  */
		   double X[],  /* [n+nlin] initial guess and returned solution */
		   int maxiter, /* maximum number of iterations */
                   double eps,  /* convergence limit */
		   int monitor); /* zero none, 1 just residuals, 3 Jacobian */