//M N
function z=mod(z1)
z=sqrt(z1(1)^2+z1(2)^2)
endfunction
function phi=arg(z1)
if z1(2)<0 then phi=-acos(z1(1)/mod(z1))
else phi=acos(z1(1)/mod(z1))
end
endfunction