The cable equation is one of the most fundamental equations for modeling neuronal dynamics. Cable equations with fractional order temporal operators have been introduced to model electrotonic properties of spiny neuronal dendrites. In this paper, we consider the following fractional cable equation involving two fractional temporal derivatives: $\u2202u(x,t)/\u2202t=D0t1\u2212\gamma 1(\kappa (\u22022u(x,t)/\u2202x2))\u2212\mu 02Dt1\u2212\gamma 2u(x,t)+f(x,t),$ where $0<\gamma 1$, $\gamma 2<1$, $\kappa >0$, and $\mu 02$ are constants, and $D0t1\u2212\gamma u(x,t)$ is the Rieman–Liouville fractional partial derivative of order $1\u2212\gamma $. Two new implicit numerical methods with convergence order $O(\tau +h2)$ and $O(\tau 2+h2)$ for the fractional cable equation are proposed, respectively, where $\tau $ and $h$ are the time and space step sizes. The stability and convergence of these methods are investigated using the energy method. Finally, numerical results are given to demonstrate the effectiveness of both implicit numerical methods. These techniques can also be applied to solve other types of anomalous subdiffusion problems.