function [c]=multh(a,b);
% function [c]=multh(a,b);
% multiplies 2 harmonics a and b
% output: c harmonic.
a1 = zeros(1,length(a));
a1(1)=a(1);
a1(2:length(a)) = a(2:length(a))/2;
b1 = zeros(1,length(b));
b1(1)=b(1);
b1(2:length(b)) = b(2:length(b))/2;
c=multbpf(a1,b1);                                
c(2:length(c))=c(2:length(c))*2;
end
