In [2]:

```
from __future__ import division
import math
#Initializing the variables
c = 1250;
Dt = 0.02;
Dv = 0.5;
rho = 1000;
v =0.5;
#Calculations
cDt = c*Dt;
Dp = rho*c*Dv;
DOv_DOt = Dv/Dt;
vDOv_DOt = v*Dv/cDt;
DOp_DOt = Dp/Dt;
vDOp_DOx = v*Dp/cDt;
Error = [vDOv_DOt*100/DOv_DOt ,vDOp_DOx*100/DOp_DOt];
print "The percentage errors are :",Error,"\nThese are very small hence can be neglected"
```

In [3]:

```
from __future__ import division
import math
#Initializing the variables
f = 0;
Atunnel = 1.227;
Ashaft = 12.57;
Q =2;
L = 200;
g = 9.81;
#Calculations
Zmax = (Q/Ashaft)*(Ashaft*L/(Atunnel*g))**0.5;
T = 2*math.pi*(Ashaft*L/(Atunnel*g))**0.5;
print "Peak water level (m) :",round(Zmax,3)
print "Mass Oscillation Period (s) :",round(T,1)
```