function he=qiuhe(x,p,y,q,z)
% he    x^p*y^q 從1－>m的和
% x,y 行向量 長度相同
% p,q 為x,y系數
%x=[2 3]; y=[ 3 4];
%p=2; q= 2;



m=length(x);
if (nargin<4 )&(m~=length(y))            %輸入量至少為四，x,y行向量長度必需一樣
    error('error check check!');
end
if nargin==4                            %沒有 z , 默認為單位行向量
    z=ones(m,1);
end

he=0;
for i=1:m
    he=he+x(i)^p * y(i)^q*z(i);     % 1-->m 求和
end