function vI_t = vI_steady(vI_0, I_Y, Ac, Bc, v, A_bar)

	options = optimset('Display', 'off', 'TolX', 1e-10, 'TolFun', 1e-10);

	vI_t = fsolve(@(vI_t) vI_t*(I_Y^(-1)) - (sqrt(v^2-vI_t^2)+vI_t)*(1-2*sqrt(Ac*Bc)) - A_bar*v^2 - 2*Bc,...
		vI_0,options);

	% Improved residual from eq. (19)