C1xx = Mxx.coeff(phi_xx)C2xx = Mxx.coeff(theta_xx)C3xx = Mxx.coeff(lambda_xx)D1xx = Mxx.coeff(phi_yy)D2xx D2Xm * D2Yn - T11xx * D2Xm * Yn) * Xm * Yn, (x, 0, a), (y, 0, b))A72 = integrate2((C2xx * D4Xm * Yn + (D2xx D2Xm * D2Yn - T11xx * Xm * D2Yn) * Xm * Yn, (x, 0, a), (y, 0, b))A75 = integrate2((C2xx * Xm * D4Yn + (D2xx