Here is if I increase the number:
Code:
GRma:2 GR[]:4308.194000 6.090400
54GR:4308.194000 6.090400
61GR:4308.194000 6.090400
54GR:4308.194000 6.090400
61GR:0.000000 6.090400
Here is if I decrease the value:
Code:
GRma:2 GR[]:4.308194 6.090400
54GR:4.308194 6.090400
61GR:4.308194 6.090400
54GR:4.308194 6.090400
61GR:0.000000 6.090400
Here is the pertinent declaration in my main() code:
Code:
int main()
{
* * *
double GRbeta[GRma][0], GRalpha[GRma][GRma];
double **beta, **alpha, *a;
double GRX, GRN;
// GRX = 430.8194;
GRX = 4.308194;
GRN = 6.0904;
GR[0] = GRX;
GR[1] = GRN;
alpha = dmatrix(1,GRma,1,GRma);
beta = dmatrix(1,GRma,0,0);
a = dvector(1,GRma);
for(i=1; i<=GRma; i++)
{
alpha[i] = GRalpha[i];
beta[i] = GRbeta[i];
}
a = GR;
oFile << "GR:" << a[0] << " " << a[1] << endl;
fGR(x, y, sig, alpha, beta, GRma, GRchi, GRnu, GRQ, a);
* * *
}
Here is the function as of now:
Code:
void fGR(double x[NSTEPS], double y[NSTEPS], double sig[NSTEPS],
double **GRalpha, double **GRbeta, int ma, double chi,
int nu, double Q, double *GRa)
{}
That is everything relevant to this problem. If there's anything I've missed, I'm basing this on the Numerical Recipes' gaussj and mrqmin functions.