function [Q]=B9A4(k)
f=@(x,y) 1./(x+y+1);
%f=@(x,y) 1 + 0*(x+y);
%xt=@(x1,x2,x3)
%yt=


for l=1:7
  h=2^(-l);
  v=[0 0; 2*h 0; h h];
  [X,Y,Wx,Wy]=gaussQuad_correct(k,v);
  %figure(l);
  %plot(X,Y,'o',v(:,1),v(:,2), [h,0], [h,0], 'r');
  Q(l)=Wx'*feval(f,X,Y)*Wy;
  %disp(['Funktionenwert für h=',num2str(h),':']);
  %disp(Q);

  %disp(feval(f,X,Y));
  %disp(X);
  %disp(Y);
end

end