sum(p)*dw*dtheta at t = 0.01 = 1.000000
sum(p)*dw*dtheta at t = 0.01 = 1.000000