def svm_model(r3d_kpc, n0, r_c, beta, r_s=1.0, gamma=3.0, epsilon=0.0, alpha=0.0):
    """
    Compute a Simplified Vikhlinin model

    Parameters
    ----------
    - r3d_kpc: array of radius in kpc
    - r_c : core radius parameter
    - n_0  : normalization
    - beta : slope of the profile
    - r_s  : characteristic radius parameter (kpc)
    - gamma : slope of the profile
    - epsilon : slope of the profile

    Outputs
    --------
    - SVM model profile as a function of the input radius vector

    """

    term1 = n0 * (1 + (r3d_kpc / r_c)**2)**(-3.0*beta/2.0)
    term2 = (r3d_kpc / r_c)**(-alpha/2.0)
    term3 = (1 + (r3d_kpc / r_s)**gamma)**(-epsilon/2.0/gamma)
    
    return term1 * term2 * term3