#pragma warning(error : 15618)
#include <origin.h>

// Add your special include files here.
// For example, if you want to fit with functions from the NAG library, 
// add the header file for the NAG functions here.

// Add code here for other Origin C functions that you want to define in this file,
// and access in your fitting function.

// You can access C functions defined in other files, if those files are loaded and compiled 
// in your workspace, and the functions have been prototyped in a header file that you have
// included above. 

// You can access NLSF object methods and properties directly in your function code.

// You should follow C-language syntax in defining your function. 
// For instance, if your parameter name is P1, you cannot use p1 in your function code. 
// When using fractions, remember that integer division such as 1/2 is equal to 0, and not 0.5
// Use 0.5 or 1/2.0 to get the correct value.

// For more information and examples, please refer to the "User-Defined Fitting Function" 
// section of the Origin Help file.


//----------------------------------------------------------
// 
void _nlsfThomasHopfield(
// Fit Parameter(s):
double N, double W,
// Independent Variable(s):
double x, double r,
// Dependent Variable(s):
double& y)
{
	const double a=1;
	const double d=1;
	// Beginning of editable part
	#include <origin.h>
	#include <oc_nag8.h>
	
	struct user
	{
	 double Nkonz,Wmax,fitX;
	};
	
	static double NAG_CALL f_callback(double r, Nag_User *comm)
	{
	 struct user *sp = (struct user *) (comm->p);
	
	double aa,bb,fitX;
	aa=sp->Nkonz;
	bb=sp->Wmax;
	fitX=sp->fitX;
	
	return (exp(-Wmax*exp(-2fitX/a))-1)*fitX*fitX;
	}
	
	void nlsfnag_integration_fitting
	(
	double W, double N;
	double x,r;
	double& y
	)
	
	{
	 double epsabs=0.0, epsrel=0.0001;
	Integer max_num_subint=400;
	double result, abserr;
	
	NagQuadProgress qp;
	
	static NagError fail;
	
	Nag_User comm;
	struct user s;
	s.Nkonz = N;
	s.Wmax= W;
	s.fitX=r;
	comm.p = (;Pointer); &s;
	
	d01smc(f_callback, Nag_UpperSemiInfinite, r, epsabs, epsrel, max_num_subint, &result, &abserr, &qp, &comm, &fail);
	if (fail.code != NE_INT_ARG_LT && fail.code != NE_BAD_PARAM && fail.code != NE_ALLOC_FAIL)
	        {
	                NAG_FREE(qp.sub_int_beg_pts);
	                NAG_FREE(qp.sub_int_end_pts);
	                NAG_FREE(qp.sub_int_result);
	                NAG_FREE(qp.sub_int_error);
	        }
	y=W*exp(-2.0*d/a)*exp(-x*W*exp(-2.0*d/a))*exp(4.0*PI*result);
	}
	// End of editable part
}

quote:

---

Originally posted by larry_lan

Hi:

Please define the structures outside the fitting function. Besides, you haven't use any independent variables in your callback function. Are you sure?

I made the code compiled successful, but not sure if the function work. You'd better check the logic by yourself.

---