function [enorm,sqnorm,renorm,rsqnorm]=webnorm(x,cc)
global n h px qx fx u_exact du_exact weight weightt
[BSPL,DBSPL]=bsplfunc;
sqnorm=0;enorm=0;renorm=0;rsqnorm=0;
xx=linspace(0,1,51);
for i=1:(length(x)-2*n-1)
    uh=zeros(1,51);sigh=zeros(1,51);
    xxx=linspace(x(i+n),x(i+n+1),51);
    ww=weight(xxx);dww=weightt(xxx);
    xk=px(xxx);xb=qx(xxx);
    for j=1:n+1,       
        BX=BSPL{j}(xx)/sqrt(h);
        DBX=DBSPL{j}(xx)/h/sqrt(h);
        uh=uh+ww.*BX*cc(i+j-1);
        sigh=sigh+(dww.*BX+ww.*DBX)*cc(i+j-1);
    end
    sqnorm=sqnorm+trapz(xxx,(uh-u_exact(xxx)).^2);
    enorm=enorm+trapz(xxx,xb.*(uh-u_exact(xxx)).^2+1./xk.*(sigh-du_exact(xxx)).^2);
    rsqnorm=rsqnorm+trapz(xxx,(uh-u_exact(xxx)).^2)/trapz(xxx,u_exact(xxx).^2);
    renorm=renorm+trapz(xxx,xb.*u_exact(xxx).^2+1./xk.*du_exact(xxx).^2);
end
enorm=sqrt(enorm);sqnorm=sqrt(sqnorm);renorm=sqrt(renorm);rsqnorm=sqrt(rsqnorm);