# Spread out flux in 2D gaussian

I want to generate a 2D Gaussian model with flux spread out for a given radius. From theory, I found out that the total integrated flux of a Gaussian model. Now I want to image it such that total flux is spread out within the Gaussian model of radius r. If I sum the flux within radius r, it should be same as the total integrated flux. Below I have given my code-

import numpy as np

(xx,yy) = (512,512) # size of an array
r=100 # r is the 2D gaussian radius
A = 100 * (1/r)**2 # A is the peak flux; 100 is the total integrated flux
xc=xx/2.; yc=yy/2. # centre of 2D Gaussian
circle=zeros((xx,yy)) # Empty array
for i in range(xx):
for j in range(yy):
circle[j,i]=  A*exp(-((i-xc)/r)**2.0 - ((j-yc)/r)**2.0)


Is this correct to divide the peak flux by solid angle of the Gaussian to spread out the peak flux within Gaussian?