1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
#ifndef T8_DIFFEQ_H
#define T8_DIFFEQ_H
#include <types.h>
/*
* `eqfs` - Массив функций
* `init` - массив начальных условий f_n(a) = init[n] (значение функций в начале интервала!)
* `files` = массив файлов, в который записывать получившиеся приближения
* `n` - количество уравнений ( и начальных условий ) ( и переменных )
* `a` - начало интервала решения ( по t? )
* `b` - конец интервала решения ( по t? )
* `dt` - шаг по времени
*/
void difeq_solve_euler(eqf_t *eqfs, f64 *init, const char **files, usz n, f64 a, f64 b, f64 dt);
/*
* То же самое, что и прошлая функция, только использует не метод эйлера, а метод рунге-кутты.
* Аргументы все точно такие же
*/
void difeq_solve_RK(eqf_t *eqfs, f64 *init, const char **files, usz n, f64 a, f64 b, f64 dt);
#endif /* T8_DIFFEQ_H */
|
