source: Sophya/trunk/SophyaExt/CodeMinuit/cfortran/Examples/minuit-main.c@ 3892

/*
 * This example demonstrate how to call MINUIT from C.
 *
 * For question/problems contact: Heplib.Support@cern.ch
 *
 */

#if defined(__hpux) || defined(_IBMR2) || defined(__alpha) || defined(__apollo) || defined(__ultrix) || defined(__sgi) || defined(__sun) 
#define extname
#endif 

#ifdef extname
#define FCN fcn_
#else
#define FCN fcn
#endif

#include "cfortran.h"
#include "minuit.h"
#define Ncont 20
int main()
{
        int error_flag=0;
        void *fcn();
        struct {
                double x[Ncont];
                double y[Ncont];
                int n;
                } pts;
        double f_null=0.;
        

        MNINIT(5,6,7);   /*  initialise  */

        MNSETI(" Minuit Example ");   /* set title */

        MNPARM(1,"X",0.,.1,f_null,f_null,error_flag);

        MNPARM(2,"-Y-",0.,.01,f_null,f_null,error_flag);

        MNEXCM(fcn,"MIGRAD",0,0,error_flag,0);

        MNEXCM(fcn,"MINOS",0,0,error_flag,0);

        MNCONT(fcn,1,2,Ncont,pts.x[0],pts.y[0],pts.n, 0);

}