In [ ]:
function my_Gale_Shap(p_prefs::Matrix{Int}, r_prefs::Matrix{Int}, p_caps::Vector{Int}, c_caps::Vector{Int})
    m = size(p_prefs, 2)
    n = size(r_prefs, 2)
    
    K = sum(p_caps)
    L = sum(p_caps)
    
    p_matched = zeros(Int64, m)
    r_matched = zeros(Int64, L)
    prop_pool = collect(1:m)