function [bv] =bvect(a,b,ip,iq,g,h) 
      
for i=1:ip+iq 
      	bv(i)=-g(i);
      	for j=1:ip+iq 
      		if j <= ip 
			t=a(j);
		else
			t=b(j-ip);
		end
		bv(i)=bv(i)+h(i,j)*t;
	end
end