function Re_theta = Suzen_et_al(lambda,Tu)
    % Suzen et al 2002

    f = @(R,T,l) (120+150*T^(-2/3))*coth(4*(0.3-100000*l/R^2))-R;
    
    Re_theta = zeros(size(Tu));

    for i = 1:length(Tu)
        T = Tu(i);
        l = 0;
        ff = @(R) f(R,T,l);
        Re_0 = fzero(ff,400);
        l = lambda(i);
        ff = @(R) f(R,T,l);
        Re_theta(i) = fzero(ff,Re_0+50);
    end

end

