subroutine fun(x, f, f8, g)
real(8) x8,f8
real(16) x16,f16

!f  = x*(sqrt(x +1) -sqrt(x))
!x8 =dble(x)
!f8 = x8*(dsqrt(x8 +1) -dsqrt(x8))
!!f16(x16) =x16*(qsqrt(x16 +1) -qsqrt(x16))
!g  = x/(sqrt(x +1) +sqrt(x))


f  = x*x*(sqrt(x +10) -sqrt(x-5))
x8 =dble(x)
f8 = x8*x8*(dsqrt(x8 +10) -dsqrt(x8-5))
g  = 15*x*x/(sqrt(x +10) +sqrt(x-5))

end

