coefficients( num_nodes, 0.0);
//цикл
for (int i = 0; i < num_basis_functions ; ++i) {
for (int j = 0; j < num_basis_functions; ++j) {
double integral = 0.0;
for (int k = 0; k < num_nodes; ++k) {
double x_k = interval_start + k * (interval_end - interval_start) / (num_nodes - 1);
integral += (a * phi1(x_k) * phi1(x_k) + b * phi1(x_k) * phi2(x_k) + c * phi2(x_k) * phi2(x_k));
}
coefficients[i] += integral;
}
}
//негізгі функциялардың интегралдарын есептеу арқылы Галеркин әдісі үшін есептеулер жүргізеді
double solution = 0.0;
for (int i = 0; i < num_basis_functions; ++i) {
solution += coefficients[i] * phi1(interval_start + i * (interval_end - interval_start) / (num_nodes - 1));
}
//шығару
cout << "Solution: " << solution << endl;
return 0;
}
Есептеу нәтижесі экранға шығарылады.