function f_L=temp0(x,T_w,L_ins,m,T_inf)
    f_L=(T_w+L_ins*m*tanh(m*x)*T_inf)/(1+L_ins*m*tanh(m*x))-273-100;
