program MetodoEuler
  implicit none

  real(8) :: x0, y0, xn, h
  real(8) :: x, y, yp, ypp, yppp
  integer :: N, i

  ! Valores iniciales
  x0 = 0.0d0
  y0 = 1.0d0

  ! Valor final de x
  xn = 1.0d0

  ! Tamaño del paso
  h = 0.05d0

  ! Número de subintervalos
  N = int((xn - x0) / h)

  ! Inicialización
  x = x0
  y = y0
  yp = 0.0d0
  ypp = 0.0d0
  yppp = 0.0d0

  ! Iteración
  do i = 1, N
    yp = yp + h * (x - y)  ! y'
    ypp = ypp + h * yp    ! y''
    yppp = yppp + h * ypp  ! y'''
    y = y + h * yp       ! Método de Euler para y
    x = x + h
  end do

  ! Resultado
  print *, "El valor de y en x =", xn, "es", y

end program MetodoEuler