for (double r1 = 0.0; r1 < 3000.0; r1 += pdf_step_size) {
double z = (r1 - mu1) / sigma1;
double gaussian = exp(-(z*z/2.0));
- curve.push_back(make_pair(r1, gaussian * simpson_integrate(OuterIntegralEvaluator(r1,mu2,sigma2,mu3,sigma3,mu4,sigma4,score1,score2,-0.5), 0.0, 3000.0, int_step_size)));
+ curve.push_back(make_pair(r1, gaussian * simpson_integrate(OuterIntegralEvaluator(r1,mu2,sigma2,mu3,sigma3,mu4,sigma4,score1,score2,-1.0), 0.0, 3000.0, int_step_size)));
}
} else {
for (double r1 = 0.0; r1 < 3000.0; r1 += pdf_step_size) {
double z = (r1 - mu1) / sigma1;
double gaussian = exp(-(z*z/2.0));
- curve.push_back(make_pair(r1, gaussian * simpson_integrate(OuterIntegralEvaluator(r1,mu2,sigma2,mu3,sigma3,mu4,sigma4,score2,score1,0.5), 0.0, 3000.0, int_step_size)));
+ curve.push_back(make_pair(r1, gaussian * simpson_integrate(OuterIntegralEvaluator(r1,mu2,sigma2,mu3,sigma3,mu4,sigma4,score2,score1,1.0), 0.0, 3000.0, int_step_size)));
}
}