program prog20 implicit none double precision T(10) integer i intrinsic dsin call tab(dsin, 0.d0, 1.d0, 10, T) do i=1,10 print *, i, T(i) end do end