function [Macu,M] = bihacutest(m,q)
for j=1:length(m),
    [U,LapU] = TestBih(m(j));
    [u,lap]=biharmlst(U,q);
    uacu = biharmonic2d(U,LapU,q);
    M(j,1) = norm(u(:)-U(:),1)/norm(U(:),1);
    M(j,2) = norm(u(:)-U(:))/norm(U(:));
    M(j,3) = norm(u(:)-U(:),inf);
    Macu(j,1) = norm(uacu(:)-U(:),1)/norm(U(:),1);
    Macu(j,2) = norm(uacu(:)-U(:))/norm(U(:));
    Macu(j,3) = norm(uacu(:)-U(:),inf);
end

     
