00001 ! 00002 implicit none 00003 ! 00004 real(8) :: pi = 3.14159265358979d+0 00005 real(8) :: x(5), y(5) 00006 integer :: ii, jj 00007 ! 00008 do ii = 1, 5 00009 x(ii) = -1.0 + 0.5*dble(ii-1) 00010 y(ii) = -1.0 + 0.5*dble(ii-1) 00011 end do 00012 do jj = 1, 5 00013 do ii = 1, 5 00014 write(6,'(5es11.3)') x(ii),y(jj),atan2(y(jj),x(ii))/pi, & 00015 & dmod(2.0*pi + atan2(y(jj),x(ii)),2.0*pi)/pi 00016 end do 00017 end do 00018 ! 00019 end