R=1;
A=1;
L=5;
dH=0.1;
H=zeros(1,11);
VC=zeros(1,11);
VB=zeros(1,11);
V=zeros(1,11);
H(0)=0;
for i=0:1:10
H(i)=(i-1)*dH
VC(i)=L*R^2*(acos(1-H(i)/R)-(1-H(i)/R)*(1-(1-H(i)/R)^2)^(1/2))
VB(i)=1/3*(R^2*A*asin((2*R*H(i)-H(i)^2)^0.5/R)-2*A*(R-H(i))*(2*R*H(i)-H(i)^2)^0.5+A/R*(R-H(i))^3*log((R+(2*R*H(i)-H(i)^2)^0.5)/(R-H(i))))
V(i)=L*R^2*(acos(1-H(i)/R)-(1-H(i)/R)*(1-(1-H(i)/R)^2)^(1/2))+2/3*(R^2*A*asin((2*R*H(i)-H(i)^2)^0.5/R)-2*A*(R-H(i))*(2*R*H(i)-H(i)^2)^0.5+A/R*(R-H(i))^3*log((R+(2*R*H(i)-H(i)^2)^0.5)/(R-H(i))))
end
A=1;
L=5;
dH=0.1;
H=zeros(1,11);
VC=zeros(1,11);
VB=zeros(1,11);
V=zeros(1,11);
H(0)=0;
for i=0:1:10
H(i)=(i-1)*dH
VC(i)=L*R^2*(acos(1-H(i)/R)-(1-H(i)/R)*(1-(1-H(i)/R)^2)^(1/2))
VB(i)=1/3*(R^2*A*asin((2*R*H(i)-H(i)^2)^0.5/R)-2*A*(R-H(i))*(2*R*H(i)-H(i)^2)^0.5+A/R*(R-H(i))^3*log((R+(2*R*H(i)-H(i)^2)^0.5)/(R-H(i))))
V(i)=L*R^2*(acos(1-H(i)/R)-(1-H(i)/R)*(1-(1-H(i)/R)^2)^(1/2))+2/3*(R^2*A*asin((2*R*H(i)-H(i)^2)^0.5/R)-2*A*(R-H(i))*(2*R*H(i)-H(i)^2)^0.5+A/R*(R-H(i))^3*log((R+(2*R*H(i)-H(i)^2)^0.5)/(R-H(i))))
end