program phi_angle integer N real phi_delta,angle print*, 'enter half-angle' read(5,*) angle print*, 'enter # of points' read(5,*) N phi_delta= 2.*angle/(N-1) do i=1,N print*, -angle+ phi_delta*(i-1) enddo stop end