р, 1979, 312с. br/>
Додаток
Текст програми
# include
# include
# include
# include
# include
# include
# define EPSILON 0.00001
# define MAXSTEP 1
# define VERSION 1.43
// ---------------------------------------- -------------------------------// f (double x, double y); do_step (double h, double x_cur, double y_cur); title (void); main (void);
// ---------------------------------------- -------------------------------// f (double x, double y)
{(x, y) (pow (2.718, x) * y);
}
// ---------------------------------------- -------------------------------// main (void)
{i; metka; flag = 0; metka1, metka2; err = 0; x0, y0; big2_step_res, super_step_res; k = 1; zoom = 1; big_step_res, small_step_res; a, b; temp ; x_cur, y_cur; h; f_max = 0, f_min = 0; norma = 0; int c = 8; * myfile; gdriver = DETECT, gmode, errorcode; (& gdriver, & gmode, ""); = graphresult (); (errorcode! = grOk)
{("Помилка ініціалізації графікі:% s n", grapherrormsg (errorcode)); ();
} (0); (0); (); ("y '= f (x, y), y (x0) = y (a) = y0, [a, b] - відрізок інтеграції n ");: printf (" na = "); ("% lg ", & a); (" b = "); ("% lg ", & b); (a> b)
{= a; = b; = temp;
} (a == b)
{("Початок відрізання інтеграції співпадає з его кінцем, повторити введення! n"); label1;
} ("y (% lg) =", a); ("% lg", & y0); (); ("[% lg,% lg] - Межі інтеграції, у (% lg) =% lg - початкова Умова. n ", а, b, а, y0); = fabs (b - a)/10; (h> 0.1) h = 0.1; _cur = a; _cur = y0; _max = y_cur; _min = y_cur; = fopen ("rk4.txt", "w"); (myfile, "Program: Ilya RK4 Version% g n", VERSION); ( myfile, "Method: Runge-Kutta n"); (myfile, "The order of method: 4 n"); (myfile, "Automatic integration step select: Enabled n"); (myfile, " ; [a, b] = [% lg,% lg], y (% lg) =% lg n ", a, b, a, y0); (x_cur <= b)
{(flag> 1) break; _step_res = do_step (h, x_cur, y_cur); = do_step (h/2, x_cur, y_cur); _step_res = do_step (h/2, x_cur + h/2, temp); = fabs (big_step_res - small_step_res); (err> EPSILON)
{= h/2;;
} _step_res = do_step (h, x_cur + h, big_step_res); _step_res = do_step (2 * h, x_cur, y_cur); (fabs (big2_step_res - super_step_res)
{* = 2;;
} (h> MAXSTEP) h = MAXSTEP; (h
{(...
