function Iter(n)

global Node r

run Node_Coor;

ND(:,:,1)=Node;

for r=1:n
    run Boundarys;
    run Arrays;
    run H_bilr;
    run F;
    run Vel;
    ND(:,:,r+1)=Node;
end
