+ double mu, sigma, probability;
+ compute_new_double_rating(mu1, sigma1, mu2, sigma2, mu3, sigma3, mu4, sigma4, score1, score2, &mu, &sigma, &probability);
+ printf("%f %f %f\n", mu, sigma, probability);
+ } else if (argc > 8) {
+ double mu3 = atof(argv[5]);
+ double sigma3 = atof(argv[6]);
+ double mu4 = atof(argv[7]);
+ double sigma4 = atof(argv[8]);
+ int k = atoi(argv[9]);
+
+ // assess all possible scores
+ for (int i = 0; i < k; ++i) {
+ double newmu1_1, newmu1_2, newmu2_1, newmu2_2;
+ double newsigma1_1, newsigma1_2, newsigma2_1, newsigma2_2;
+ double probability;
+ compute_new_double_rating(mu1, sigma1, mu2, sigma2, mu3, sigma3, mu4, sigma4, k, i, &newmu1_1, &newsigma1_1, &probability);
+ compute_new_double_rating(mu2, sigma2, mu1, sigma1, mu3, sigma3, mu4, sigma4, k, i, &newmu1_2, &newsigma1_2, &probability);
+ compute_new_double_rating(mu3, sigma3, mu4, sigma4, mu1, sigma1, mu2, sigma2, i, k, &newmu2_1, &newsigma2_1, &probability);
+ compute_new_double_rating(mu4, sigma4, mu3, sigma3, mu1, sigma1, mu2, sigma2, i, k, &newmu2_2, &newsigma2_2, &probability);
+ printf("%u-%u,%f,%+f,%+f,%+f,%+f\n",
+ k, i, probability, newmu1_1-mu1, newmu1_2-mu2,
+ newmu2_1-mu3, newmu2_2-mu4);
+ }
+ for (int i = k; i --> 0; ) {
+ double newmu1_1, newmu1_2, newmu2_1, newmu2_2;
+ double newsigma1_1, newsigma1_2, newsigma2_1, newsigma2_2;
+ double probability;
+ compute_new_double_rating(mu1, sigma1, mu2, sigma2, mu3, sigma3, mu4, sigma4, i, k, &newmu1_1, &newsigma1_1, &probability);
+ compute_new_double_rating(mu2, sigma2, mu1, sigma1, mu3, sigma3, mu4, sigma4, i, k, &newmu1_2, &newsigma1_2, &probability);
+ compute_new_double_rating(mu3, sigma3, mu4, sigma4, mu1, sigma1, mu2, sigma2, k, i, &newmu2_1, &newsigma2_1, &probability);
+ compute_new_double_rating(mu4, sigma4, mu3, sigma3, mu1, sigma1, mu2, sigma2, k, i, &newmu2_2, &newsigma2_2, &probability);
+ printf("%u-%u,%f,%+f,%+f,%+f,%+f\n",
+ i, k, probability, newmu1_1-mu1, newmu1_2-mu2,
+ newmu2_1-mu3, newmu2_2-mu4);
+ }