_jh_B_1_47b.wxmx

kill(all);

ro:1000;g:9.81;R:2;B:1;

P_x

P_x:integrate(ro*g*(R-h)*B,h,0,R);

m:0.5;P_x*m;

P_z

W:(R**2-((%pi*(R**2))/4))*B;

P_z:W*ro*g,numer;

P:sqrt(P_x**2+P_z**2);

P/2;

fi_rad:atan(P_z/P_x);

fi:(180/%pi)*fi_rad,numer;

z_C

fb(h):=(R-sqrt(R**2-h**2));

fz(h):=(R-(fb(h)/2));

S1:integrate(fb(h),h,0,R);

mS:integrate((fb(h)*fz(h)),h,0,R);

z_c:mS/S1,numer;