In this work we present a large signal non quasi static (LSNQS) compact model based on the discretization of the current continuity equation using a spline collocation approach. The underlying charge/current model is based on the VRH theory. The LSNQS model takes in consideration the presence of parasitic regions and is particularly suited for simulations of printed organic thin film transistors. Some results obtained by a Verilog-A implementation are presented.