We present a method to obtain the frequency response of step index (SI) plastic optical fibers (POFs) based on the power flow equation generalized to incorporate the temporal dimension where the fibre diffusion and attenuation are functions of the propagation angle. To solve this equation we propose a fast implementation of the finite-difference method in matrix form. Our method is validated by comparing model predictions to experimental data. In addition, the model provides the space-time evolution of the angular power distribution when it is transmitted throughout the fibre which gives a detailed picture of the POFs capabilities for information transmission. Model predictions show that angular diffusion has a strong impact on temporal pulse widening with propagation.