def METHOD_NAME(self, x):
    parameters = self.parameters
    Vcs = x['qcs']/parameters['Cs']
    Vcp = x['qcp']/parameters['Ccp']
    SOC = (parameters['CMax'] - parameters['qMax'] + x['qb'])/parameters['CMax']
    Cb = parameters['Cbp0']*SOC**3 + parameters['Cbp1']*SOC**2 + parameters['Cbp2']*SOC + parameters['Cbp3']
    Vb = x['qb']/Cb
    return self.OutputContainer(np.array([
        np.atleast_1d(x['tb']),            # t
        np.atleast_1d(Vb - Vcp - Vcs)]))   # v