function [FPperOD, dFPperODdt, degradation, synthesis] = FPProcess(t,OD,gammad,FP)

FPperOD = FP./OD;
dFPperOD = diff(FPperOD);
dt = mean(diff(t));
dFPperODdt = dFPperOD./dt;
degradation = gammad.*FPperOD(2:end,:,:);
synthesis = dFPperODdt + degradation;