def METHOD_NAME(zetaC, ZetaPanel, vortex_radius, gamma=1.0):
    """
    Induced velocity over point ZetaC of a panel of vertices coordinates
    ZetaPanel and circulaiton gamma, where:
        ZetaPanel.shape=(4,3)=[vertex local number, (x,y,z) component]
    """
    q = np.zeros((3,))
    for ss, aa, bb in zip(svec, avec, bvec):
        q += biot_segment(zetaC, ZetaPanel[aa, :], ZetaPanel[bb, :],
                          vortex_radius, gamma)
    return q