Проблема с динамическим выделением памяти

[Bugrimov]

Добрый день!
Даже не знаю как описать проблему....
Программа выполняет решение дифференциального уравнения методом Рунге-Кутта.
Проблема в динамическом выделении памяти. На сколько я могу судить неверно выделяется динамическая память... Скажу даже происходит выход за пределы массива x[i + 1]. Решение проблемы не по глазам... К тому же компилятор возмущается...
warning C4244: функция: преобразование 'double' в 'size_t', возможна потеря данных естественно в строке где выделяется память под массивы. Сможет кто-нибудь намякнуть, как разрешить проблему в коде..???

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

/* ------- Прототипы функций ------- */
double myFunc(double, double);

int main()
{
	double  a = 0.0;	// Начальное значение
	double  b = 1.0;	// Конечное значение
	int		i;			// Счетчик
	double	n;
	double  h = 0.05;	// Шаг
	double	*x = NULL; 
	double	*y = NULL;
	double	k1, k2, k3, k4;
 
	n = (b - a) / h; // Количество шагов 
 
	x = (double*)malloc(n * sizeof(double)); // Выделение памяти под массив X
	y = (double*)malloc(n * sizeof(double)); // Выделение памяти под массив Y
	
	if(x == NULL || y == NULL)
	{
		printf("Error memory!!!");
		exit(1);
	}
 
	x[0] = a;
	y[0] = 0;
 
	for(i = 0; i <= n; i++)
	{
		k1 = myFunc(x[i], y[i]);
		k2 = myFunc(x[i] + h * 0.5, y[i] + (h * 0.5 * k1));
		k3 = myFunc(x[i] + h * 0.5, y[i] + (h * 0.5 * k2));
		k4 = myFunc(x[i] + h, y[i] + h * k3);
		x[i + 1] = x[i] + h;
		y[i + 1] = y[i] + (h / 6) * (k1 + (2 * k2) + (2 * k3) + k4);
	}
 
	for(i = 0; i <= n; i++)
		printf("x = %.3lf - y = %.3lf\n", x[i], y[i]);
	 
	free(x); // Освобождение памяти
	free(y);
	 
	return 0;
}

/* ------- Дифференциальное уравнение ------- */
double myFunc(double x, double y)
{
	return cos(7 * x + y) + 7 * (x - y);
}
/* ------------------------------------------ */

[Аватар]

циклы то до n-1 нужно делать. Тем более присвоение в x[i + 1] и y[i + 1] есть, даже так будет выход за пределы

[Bugrimov]

Сообщение от Аватар

циклы то до n-1 нужно делать. Тем более присвоение в x[i + 1] и y[i + 1] есть, даже так будет выход за пределы

Если посмотреть то конечно понятно, алгоритм который мне скинули противоречит... Может я чего-то не до понял... Глаз замылился..
http://www.simumath.net/library/book...f_Ur_method_RK
Может есть какой-нибудь вариант???

[Аватар]

на вскидку сделайте 1-ый цикл до n-2, 2-ой до n-1

[Bugrimov]

Как вариант.... Одной итерации не хватает для расчета x = 1.0. Что можно придумать.
Может увеличить объем выделяемой памяти...???
Например

x = (double*)malloc(n+1 * sizeof(double));

[Аватар]

Память для n+1 и оба цикла до n-1? Можно и так

x = (double*)malloc((n+1) * sizeof(double));

Не?

[Bugrimov]

Работает, но что-то не то... Наверно слишком дотошный
Может сам алгоритм реализации неправильный...?
А что с предупреждением компилятора (пост №1)

[Аватар]

А что с предупреждением компилятора

Не скажу, код понимаю, но не мой язык. Алгоритмически похоже верно

[Bugrimov]

Сообщение от Аватар

Не скажу, код понимаю, но не мой язык

А по поводу алгоритма выполнения??? Может как-то по другому построить...?

[Аватар]

А что смущает? Разве что в лишней точке вычисления. Но она за пределами отрезка справа. А в алгоритме по ссылке описка, должно быть Yi+1=Yi + ...