MODULE functions CONTAINS FUNCTION my_exp(x) result(y) real(kind=8), intent(in) :: x real(kind=8) :: y y = exp(x) END FUNCTION my_exp FUNCTION my_sin(x) result(y) real(kind=8), intent(in) :: x real(kind=8) :: y y = sin(x) END FUNCTION my_sin FUNCTION my_poly(x) result(y) real(kind=8), intent(in) :: x real(kind=8) :: y !! for the task some polynomial of order 10 should be tested... !! so to have a simple equation I chose Chebyshev polynomial y = cos(10 * acos(x)) END FUNCTION my_poly END MODULE functions