Thanks for your reply. I tried “ erru3 = sqrt( Integrate(errDD,mesh,order = 2*order, definedon=[2]) )” , but it doesn't seem to be able to use it that way.
TypeError: Integrate(): incompatible function arguments. The following argument types are supported: