A numerical scheme for a partial differential system motivated by light-triggered drug delivery

•Study of a system of quasilinear PDEs that can be used to model drug release from a polymeric structure enhanced by light.•A numerical scheme to compute the numerical solution of the IBVP based on the system of PDEs is proposed.•The theoretical support (stability and convergence) under weaker assumptions than those usually used is developed.•Numerical results illustrating the theoretical results established are provided.•A comparation between experimental data and simulation results is also presented.