function volume( theta, height)
    return math.pi * h^2 * math.tan ( theta /2 ) ^ 2
end