координаты в matlab

clear 
 
Nx=3; % количество элементов по оси x
Ny=3; % количество элементов по оси y
dx=1; % расстояние между элементами по оси x
dy=1; % расстояние между элементами по оси y
dpsix=10.*pi/180; % фазовый сдвиг между элементами по оси x
dpsiy=15.*pi/180; % фазовый сдвиг между элементами по оси y
 
N=Nx*Ny; % количество элементов
x=                             % координата элемента по x        
y=                             % координата элемента по y 
z=0;                           % координата элемента по z 
Am=0;                         % Am данного элемента 
P=                             % фазовый сдвиг между элементами 
 
N=;
x1=;   y1=;    z1=;   Am1=;   P1=; 
x2=;   y2=;    z2=;   Am2=;   P2=;
xn=;   yn=;    zn=;   Amn=;   Pn=;

clear 
 
Nx=3;                   % количество элементов по оси x
Ny=3;                   % количество элементов по оси y
dx=10;                  % расстояние между элементами по оси x
dy=10;                  % расстояние между элементами по оси y
dpsix=30;               % фазовый сдвиг между элементами по оси x
dpsiy=30;               % фазовый сдвиг между элементами по оси y
 
N=Nx*Ny; % количество элементов
q=0:(Nx*Ny)-1
x=0:dx:dx*(Nx-1);
y=0:dy:dy*(Ny-1);
[xe ye]=meshgrid(x,y);
x=reshape(xe',1,[]); 
y=reshape(ye',1,[])
z=(0:1:N-1)*0
Am=(0:1:N-1)*0
Px=0:dpsix:dpsix*(Nx-1);
Py=0:dpsiy:dpsiy*(Ny-1);
[P1 P2]=meshgrid(Px,Py);
Px=reshape(P1',1,[])
Py=reshape(P2',1,[])
% q=q', x=x', y=y', z=z', Am=Am', Px=Px', Py=Py';
A=[q;x;y;z;Am;Px;Py];
B=A'
save('elementsave_system.txt','B','-ascii', '-tabs')

%  save('elementsave_system.txt','B','-ascii', '-tabs')
 dlmwrite('elementsave_system.txt', B)