У меня есть двадцать уравнений, с которыми я работал от Hodgkin Huxley model до cressman model, и я хотел бы проанализировать модель, которую я получил с помощью числовых методов.
Я математик и плохо разбираюсь в кодировании, поэтому, если бы кто-нибудь мог помочь мне закодировать эти уравнения с помощью scilab, используя runge kutta IV method, я был бы очень благодарен.
В ответ я могу поделиться своей работой с нейронами activiy, если вам интересно.
С Уважением.
Вы можете использовать без программирования что-либо, кроме пользовательской функции ODE, прямо из коробки,
y = ode(y0, t0, times, rtol, atol, odefunc)
где (t0,y0) - начальные точки, times список раз, когда образцы решения сообщаются в y, а function dydt = odefunc(t, y) - функция, реализующая производные.
Поскольку за этим стоит проверенный и надежный алгоритм, вы можете предположить, что фактические числовые ошибки значимым образом связаны с допусками rtol, atol.
Примеры см. Во встроенной документации и интерактивных примерах, таких как Решение ODE в Scilab, Использование Scilab для решения и построения дифференциальных уравнений
odeв scilab? RK4 с фиксированным шагом, хотя и неплохой, почти всегда неоптимален по сравнению с решателями с адаптивным размером шага. - person Lutz Lehmann schedule 25.05.2017