add functions for double porosity model with 3 parameters