%C_function
function y = C_function(kappa, t)
global L0
y = L0/(L0 + (kappa(2)/kappa(1))).*(1-exp(-(kappa(1).*L0+kappa(2)).*t));