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

x8 =x
x16 =x

! Cod pentru f; f8 sau f16; g;
f  = 0.1*x*(sqrt(x +100) -sqrt(x))
f8 = 0.1*x8*(dsqrt(x8 +100) -dsqrt(x8))

f16 =0.1*x16*(qsqrt(x16 +100) -qsqrt(x16))
g  = 10*x/(sqrt(x +100) +sqrt(x))


end

