function B=rightmatrix(x,p,y,q,z)
% A*a=B 
% B  為一個列向量 長為p*q
% x y z 為點的坐標 
%p q 為x y指數

if nargin~=5 
    error('error check check! rightmatrix')
end

B=zeros(p*q,1);
for i=1 : p*q
    x_z=quotient(i-1,q);
    y_z=mod(i-1,q);
    B(i,1)=qiuhe(x,x_z,y,y_z,z);
end