clear; clc;
Aimax=0.1;
Ao=0.001;M=100;g=9.8;L=1;rho=1000;
Ai=linspace(3*Ao,Aimax,100);
vo=((M+L*rho.*Ai)*g./(2*rho*Ao*(1-2*Ao./Ai))).^0.5;
plot(Ai,vo);
xlabel('Ai');ylabel('vo');
