from numpy import zeros, where

def ReactionTitrante (C_ATitrer, C_Titrante, V_ATitrer, V_Titrante):
    n1 = C_ATitrer * V_ATitrer
    n2 = C_Titrante * V_Titrante
    n1p= zeros(n2.shape)
    n2p= zeros(n2.shape)
    AvantEquiv = where (n1 > n2)
    n1p[AvantEquiv] = n1- n2[AvantEquiv]
    ApresEquiv = where (n2 > n1)
    n2p[ApresEquiv] = n2[ApresEquiv]-n1
    return n1p, n2p