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]
    """
    vortex_radius_sq = vortex_radius*vortex_radius
    Cfact = cfact_biot * gamma
    q = np.zeros((3,))
    R_list = zetaC - ZetaPanel
    Runit_list = [R_list[ii] / algebra.norm3d(R_list[ii]) for ii in svec]
    for aa, bb in LoopPanel:
        RAB = ZetaPanel[bb, :] - ZetaPanel[aa, :]  # segment vector
        Vcr = algebra.cross3(R_list[aa], R_list[bb])
        vcr2 = np.dot(Vcr, Vcr)
        if vcr2 < (vortex_radius * algebra.normsq3d(RAB)):
            continue
        q += ((Cfact / vcr2) * np.dot(RAB, Runit_list[aa] - Runit_list[bb])) * Vcr
    return q