syms z h r s1 b a
z0=h+s1*tan(r)-b;
x=a/b*sqrt(b^2-z^2);
y=(z0-z)/tan(r);
w=2*x*y;
v=int(w,z,-b,z0)
运行结果:
Warning: Explicit integral could not be found.
> In sym.int at 58
   In test at 11
v =
int(2*a/b*(b^2-z^2)^(1/2)*(h+s1*tan(r)-b-z)/tan(r),z = -b .. h+s1*tan(r)-b)