[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