!Главная программа program test use eyler implicit none real (8) :: x0 !Начальное X real (8) :: x1 !Конечное X real (8) :: y0 !Начальное Y real (8) :: y1 !Конечное Y integer :: n !Число шагов интегрирования integer :: error !Обшибка метода Ньютона (зацикливание) !Начальные условия x0 = 0.0 x1 = 3.0 y0 = 1.0 n = 500 call eyler_n (x0, x1, n, y0, y1, error) print '(a, 1p, e15.7)', "Решение = ", y1 print '(a, i3)', "Ошибка = ", error end program test !Функция диф. уравнения pure real (8) function f (x, y) real (8), intent (in) :: x, y f = 50 * (exp(-x)-y) end function f