In a recent article I read, which referred to the ODE model of the photopolymerisation reaction, the relationship between the photoinitiator concentration and the time during the initiation phase of the reaction needs to be calculated, but this equation also includes the reaction rate kd, which is related to both the photoinitiator concentration, and the depth z. How should I calculate this equation? I tried it in MATLAB and got results, but wish I could model it in COMSOL.