include(`test/CalculusFunctions.txt'):

let(pi,4*atan(1)):

let(r,<<`H/D ratio',`Vol. ratio'>>):

for(<i,0,1,.2>,{
  block({
    let(depth,(i)):
    let(temp,.5-(depth)):
    let(theta,newtonsolve(
      lambda(
        {((pi)/2)*(temp)-(1/4)*sin(2*[x])-(1/2)*[x]}
      ),.5,6)
    ):
    let(ratio,(1-sin(theta))/2):
    let(r,(r)~push~<(ratio),(depth)>):
    print_():
  })
}):

print_(r)
