solution=resolution(U0, C1_0, C1_L, L, dt, t_max)[/CODE]