C++, метод Эйлера с усреднением

[lana_ar]

#include <iostream>
#include <math.h>
using namespace std;

double f(double x, double y, double z)
{
	return 3 * pow(x, 2) - x - 1 + 6 * y - z;
}


int main()
{
	double a, b, n, h, xn = 0, yn = -0.9, zn = 3.2;
	const int size = 6;
	cin >> a;
	cin >> b;
	cin >> n;

	double x[size];
	double y[size];
	double z[size];
	x[0] = xn;
	y[0] = yn;
	z[0] = zn;
	h = (b - a) / n;
	for (int i = 1; i < size; i++)
	{
		x[i] = x[i - 1] + h;
		y[i] = y[i - 1] + h * (z[i - 1] + h * f(x[i - 1], y[i - 1], z[i - 1]));
		z[i] = z[i - 1] + h * ((z[i - 1] + h * f(x[i - 1], y[i - 1], z[i - 1])) + f(x[i - 1] + h, y[i - 1] + h * z[i - 1], z[i - 1] + h * f(x[i - 1], y[i - 1], z[i - 1])))/2;

		
			
	}
	for (int i = 0; i < size; i++)
	{
		cout << i << "\t";
		cout << x[i] << "\t";
		cout << y[i] << "\t";
		cout << z[i] << "\t\n";
		
	}
	return 0;
}

Считает неправильные значения, помогите, пожалуйста
Есть предположение, что накосячила в строчках, которые выделила

[Алексей_2012]

Почему один цикл for начинается с i=1, а другой уже с 0?
С какого индекса массив заполняется с первого или с нудевого, Вы в курсе?

[lana_ar]

с 0, мой косяк но мне сказали, что проблема именно в функциях y[i] и z[i] различны

[BDA]

Во-первых, я считаю, что заполнение именно с единицы, так как для нулевой точки данные даны. Во-вторых, если size явно задан, то нет смысла запрашивать n. Или стоит динамически выделять память под столько точек, сколько задаст пользователь. В-третьих, я не очень дружу с математикой, а на картинках нет формулы для расчета z, так что могу судить только по формулам, приведенным в статье на википедии о модифицированном методе Эйлера с пересчетом, так что просто подставил значение zn.

x[i] = x[i - 1] + h;
x_pr = x[i];
y_pr = y[i - 1] + h * f(x[i - 1], y[i - 1], zn);
y[i] = y[i - 1] + h * (f(x[i - 1], y[i - 1], zn) + f(x_pr, y_pr, zn)) / 2;

[lana_ar]

спасибо огромное! распуталась!