Hickerson series of almost integers: Difference between revisions

m
→‎{{header|Fortran}}: even if there is a real(8) kind, "q" is not standard conformant: replaced this
(Ada version)
m (→‎{{header|Fortran}}: even if there is a real(8) kind, "q" is not standard conformant: replaced this)
Line 576:
 
=={{header|Fortran}}==
<lang fortran>Programprogram Hickersonhickerson
implicit none
! 3 February 2014
integer, parameter :: q = selected_real_kind(30)
! not all Fortran compilers provide REAL*16 and INTEGER*8
integer, parameter :: l = selected_int_kind(15)
implicit none
real(kind=kind(1q0)q) :: s, l2
integer(kind=kind(1_8)) :: i, n,f,is k
 
l2 = log(2.0_q)
do n = 1, 17
do sn = 0.5q0 /1, log(2q0)17
do i s = 1,n0.5_q / l2
do si = (s * i) /1, log(2q0)n
end do s = (s * i) / l2
is = send do
 
k = floor((s - floor(s, l)) * 10)
if (fk == 0 .or. fk == 9) then
write(* print 1,10) n, s,'' ""
else
write(* print 1,10) n, s,' " NOT'"
endif
end do
10 1 format('h(',i2I2,') = ',F23.3,' is',A,' an almost-integer')
end program Hickerson</lang>
 
is = s
f = (s-is)*10 !first digit after decimal point
if (f == 0 .or. f == 9) then
write(*,10)n,s,''
else
write(*,10)n,s,' NOT'
endif
end do
10 format('h(',i2,') = ',F23.3,' is',A,' an almost-integer')
end program Hickerson</lang>
{{out}}
 
<pre>
h( 1) = 1.041 is an almost-integer
Line 619 ⟶ 621:
h(17) = 130370767029135900.458 is NOT an almost-integer
</pre>
 
=={{header|FreeBASIC}}==
{{libheader|GMP}}
1,336

edits