Modeling gene regulatory networks with piecewise linear differential equations