We introduce a new state variable v, the velocity, so y′=v. Then we have a pair of equations y′=v and v′=-βv-γy+u. We put these together in a matrix equation
so the differential equation is