I have the following grid in the beginning of my code:
a = 25,
r = linspace(0,sqrt(2)*2*a,501);
phi = linspace(0,2*pi,501);
[R,Phi] = meshgrid(r,phi);
X = R.*cos(Phi);
Y = R.*sin(Phi);
And I plot all my figures by using the code segment:
h = pcolor(X,Y,abs(u_sum).^2);
The output of this code is:
Now, I want to draw a circle full of ones with radius "a" and place it on the same grid. Thus, I want the circle at the center of the above picture to be filled with ones, and zeros at rest. How can I do that using pcolor and the same grid?
Alternatively, how can I implement a Gaussian shape for the same grid using pcolor?