Simulation of PK models I