function funImpl (valorX, valorY, aproxraiz, h) result (res)
    implicit none
    real (kind = 8) :: res, valorX, valorY, aproxraiz, h, funcion

    res = valorY + h * (0.005_8 * log ((100000._8) / aproxraiz) * aproxraiz) - aproxraiz
end function funImpl

function derFunImp (valorX, valorY, aproxraiz, h) result (res)
    implicit none
    real (kind = 8) :: res, valorX, valorY, aproxraiz, h, k, k2, funcion

    res = h * (0.005_8 * log (100000._8 / aproxraiz) - 0.005_8) - 1._8
end function derFunImp