Body := proc (R::list, Bo::list, Hi, hf, kf) local i, q1, q2, q3, q4, q5, q6, w1, w2, w3, w4, w5, w6, w7, w11, w8, w9, w10, w12, w13, w14; global Rho, Cf, Sf, X, Y; q3 := -u + theta; q5 := T^2; w2 := 0; for i to n do q1[i] := R(u)[i]; q2[i] := q1[i]^2; q4[i] := q1[i]*T*cos(q3); q6[i] := ((1/hf*Hi - 1)/Pi*kf*So + q2[i] - 2*q4[i] + q5)^(1/2); w2 := w2 + evalf(Int(q6[i]*(T*diff(R(u)[i], u)*sin(q3) + q2[i] - q4[i])/(q2[i] - 2*q4[i] + q5), u = Bo[i].. Bo[i+1], EIco)); od; w1 := hf^2; w3 := 5^(1/2); w4 := w2^2; w5 := -1 + 1/hf*Hi; w6 := (5*w4 - 24*Pi*So*w5*(-1 + kf))^(1/2); w7 := 1/Pi; w11 := T^2; for i to n do w8[i] := R(u)[i]; w9[i] := w8[i]^2; w10[i] := w8[i]*T*cos(q3); w12[i] := (w5*w7*kf*So + w9[i] - 2*w10[i] + w11)^(1/2); w13[i] := (w9[i] - 2*w10[i] + w11)^(1/2); w14[i] := 1/w13[i]; od; Rho := [seq(unapply(1/2*z*(-z + hf)/w1*(-5*w2 + w3*w6)*w7 + w12[i], u, z), i = 1 .. n)]; Cf := [seq(unapply((w8[i]*cos(u) - T*cos(theta))*w14[i], u), i = 1 .. n)]; Sf := [seq(unapply((w8[i]*sin(u) - T*sin(theta))*w14[i], u), i = 1 .. n)]; X := [seq(unapply(Rho(u, z)[i]*Cf(u)[i], u, z), i = 1 .. n)]; Y := [seq(unapply(Rho(u, z)[i]*Sf(u)[i], u, z), i = 1 .. n)]; print(`To display the compressed body use the command:`); if T = 0 then print(`cylinderplot(r(f,z), f=a..b, z=0..h)`); else print(`plot3d([x(f,z), y(f,z), z], f=a..b, z=0..h)`); fi; end;