Runge Kutta с использованием массивов в C ++. Как правильный способ объявить мои переменные? ⇐ C++

[Anonymous]

Я пытаюсь имитировать цепь связанных осцилляторов и для этого мне нужно использовать Runge-Kutta Order 4 с массивом (для всех различных генераторов). Моя проблема в том, что я действительно плоха в правильном объявлении переменных, и RK4 не делает то, что должно быть, вероятно, потому что я использую указатель для расчета. (Используя один, потому что он не будет работать без него, почему). Если бы кто -нибудь знал, как я могу решить эту проблему, это было бы потрясающим. Спасибо! < /P>
Мои функции F1 и F2 находятся в файле заголовка. Любая функция, которая не зависит от моего индекса, я работаю совершенно хорошо. (Как правило, если F2 возвращает 1, я получаю, как и ожидалось, 1/2x^2 для второго интеграла). Как только я положил F2, это зависит от того, что я всегда получаю что -то, что расходится. Код ниже: (Удачи) < /p>
Файл заголовка: < /p>
#include
#include
using namespace std;

double f1(double theta1) {
return theta1;
}

double f2(double t, double theta[], double alpha, double beta, int i) {
//return sin(theta)-beta*(theta[i+1]+theta[i-1]-2*theta);
return sin(theta);
}
< /code>
Основной файл: < /p>
#include "func_solition.h"
#include
using namespace std;

int main() {
// Paramètres intégrale
double dt = 0.01;
double t = 0;
double TF = 20;

// CI pendules
int n = 10; // ATTENTION ON A DEUX PENDULES POUR LES BC, PREVOIR +2 PENDULES
double theta[n];
double theta1[n];
for (int i = 0; i < n; i++) {
theta = 0;
}
theta[1] = 3.141592/2;
for (int i = 0; i < n; i++) {
theta1 = 0;
}
theta1[1] = 0;

// Constantes équation du mouvement
double alpha = 0.1;
double beta = 0.01;

// Variables pour rk4
double tmp1=0.0, tmp2=0.0;
double k1, k2, k3, k4;
double l1, l2, l3, l4;

cout

Подробнее здесь: https://stackoverflow.com/questions/794 ... -variables