I'm a student that recently started taking a course on cognitive robotics. The book I use is Probabilistic Robotics by Thrun Burgard and Fox.
In the EKF algorithm, we linearized the action model in the following way
$$g(u_t,\:x_{t-1})\:=\:g(u_t,\:\mu _{t-1})+G_t\cdot (x_{t-1}-\mu \:_{t-1})$$
$g(u_t,\:x_{t-1})$ is the action model and $G_t$ is its Jacobian matrix with respect to the state $x_{t-1}$.
I don't see how this guarantees linearity because $g$ could be nonlinear in $u_t$. The authors don't mention anything about why this is the case.
In other words, I imagined that the multivariate taylor expansion for this where we get a linear function in both $u_t$ and $x_{t-1}$