Parameters as a function of time

I apologize but I don’t quite understand precisely what you want to achieve. Are you having issues with determining the correct distribution to implement or is it the translation of the model into code that is the problem?