[cdev] Метод RKF45

Grigoriy A. Sitkarev sitkarev на komitex.ru
Вт Апр 13 11:27:17 MSK 2010


Здравствуй, Андрей.

Размеры шагов конечно же ты можешь сам задавать, это именно то что 
пользователь задаёт для управления точностью метода. У тебя в исходниках 
заданы очень большие размеры минимального шага, это 0.1 всего и точность 
  epsilon тоже не большая. Попробуй уменьшить минимальный шаг до 0,025 
хотя бы и epsilon на несколько порядков. Это для начала. Если не пойдёт, 
я посмотрю ещё раз где я мог ошибиться. И перепроверю твою программу. 
Кстати, неплохо "выглядишь", но я потом ряд замечаний выскажу, давай 
доведём сначала метод до решения.

Потом сделаем обвес, чтобы удобно можно было завернуть всё для решения 
задач на практике.

--
Г.А.

Андрей Ануфриев пишет:
> Здравствуйте. Я посмотрел код и нашёл у себя две ошибки, но всё равно результат вычисления не стал точным. Я вынес задачу для вычисления методом РК 4-го порядка в отдельный файл (RK.c), там вычисления идут как в простом методе РК 4-го порядка, который мы писали раньше, но результаты отличаются почти на 0,1. Так-же переделал RKF45. Можно ли в нём самим задавать максимальный или минимальный шаг,или это должен делать пользователь?
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> cdev mailing list
> cdev на wiki.syktsu.ru
> http://wiki.syktsu.ru/cgi-bin/mailman/listinfo/cdev





Подробная информация о списке рассылки cdev