function wmp = win_mod_poly(N, t, poly, w)
R = size(poly,1);
ply_val = zeros(N,R);
w_mat = repmat (w, 1, R);
for k = 1:R
  ply_val(:,k) = polyval(poly(k,:),t).';
end
wmp = sum (w_mat.*exp(ply_val),1);
end