We discuss, both theoretically and numerically, the propagation of self-confined beams in uniaxial nematic liquid crystals. Our model accounts for optical anisotropy and goes beyond the perturbative regime, that is, it allows addressing the unique properties of nematicons, including their tunable nonlinearity, high nonlocality, and ability to spatial self-steer.