subroutine size(sigmas) real sigmas vector x(640),y(640) real theta do i=1,640 theta = i/100. x(i) = 1.7*sigmas*cos(theta) y(i) = 0.7*sigmas*sin(theta) enddo end