function F = fun_sol(s)
global T_r T_l match_full
F = match_full - T_r*T_l*s;