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

    !res = valorY + h * (0.0152_8 * aproxraiz * (1._8 - (aproxraiz / 100000._8))) - aproxraiz
    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 * (aproxraiz * (-19._8 / 6250000._8) + 19._8 / 2500._8) - 1._8
    res = h * (0.005_8 * log (100000._8 / aproxraiz) - 0.005_8) - 1._8
end function derFunImp