The attached file has a 3d curve fit of pressure data on a disc. The curve fitting routine could only do so in x,y,z coordinates and the this is best handled in r,z theta. I replaced x with r*cos(theta) and y with r*sin(theta). The integral will not work. I believe that it is having a problem because it is evaluating one variable at a time. Does anyone have a method to integrate this surface.