h = [0.1 0.05 0.025 0.0125 0.00625 0.003125]
eL2 = [0.00082538 0.00020652 5.164e-5 1.2911e-5 3.2277e-6 8.0692e-7]
eH1 = [0.051577 0.025796 0.012899 0.0064494 0.0032247 0.0016124]
degree = 1;
figure(10); clf
            h1=loglog(h,eL2);
            hold on
            h2=loglog(h,eH1,'--');
            %
            aa=h(4);
            e1=eL2(end-1);
            bb=aa*6;
            e2=e1*exp(2*log(6));
            loglog([aa bb],[e1,e2],'LineWidth',2)
            loglog([aa bb],[e1,e1],'LineWidth',1)
            loglog([bb bb],[e1,e2],'LineWidth',1)
            %
            aa=h(5);
            e1=eH1(end);
            bb=aa*6;
            e2=e1*exp(1*log(6));
            loglog([aa bb],[e1,e2],'LineWidth',2)
            loglog([aa bb],[e1,e1],'LineWidth',1)
            loglog([bb bb],[e1,e2],'LineWidth',1)

            

qL2 = log(eL2(2:6)./eL2(1:5)) ./ log(h(2:6)./h(1:5))
qH1 = log(eH1(2:6)./eH1(1:5)) ./ log(h(2:6)./h(1:5))