Hey, I have defined
pu = zeros(nx,ny,N);
pv = zeros(nx,ny,N);
pu(:,ny,:) = 1;
and written the loop:
line 38-40
for i = 2:(nx-1);
for i = 2:(ny-1);
ps(i,j) = p(i,j,n) - a1*(pu(i+1,j,n) - pu(i,j,n)) -a2*(pv(i,j+1,n) - pv(i,j,n));
and I get the error
Could someone help me...