Index exceeds matrix dimensions.

I have the following portion:

N = size(p,1);
x = 1:N; y = x;
[X,Y] = meshgrid(x,y);
R = sqrt((2.*X-N-1).^2+(2.*Y-N-1).^2)/N;
Theta = atan2((N-1-2.*Y+2),(2.*X-N+1-2));
R = (R<=1).*R;
Rad = radialpoly(R,n,m);    

Product = p(x,y).*Rad.*exp(-1i*m*Theta);

I get the error in the following line:

Product = p(x,y).*Rad.*exp(-1i*m*Theta);

The error that I get is:

Index exceeds matrix dimensions.

I tried to "debug" the program, and found that the dimensions issue may be due to the following line:

    [X,Y] = meshgrid(x,y);

The weird thing is that the program doesn't give me errors on most of the images I'm applying it on.

Do you know why I'm getting such error?

Thanks.

[Zoltán Csáti]

What is p? Does p have the size of NxN?

[med-sweng]

"p" is just a matrix passed to a function

[med-sweng]

Regarding the size, no, it is not of NXN size