基于几何光学的光线折射和反射原理提出了变折射率介质中的光线追迹算法,给出了在光线追迹过程中进行递推的反射光线和折射光线的方向余弦方程;解决了空间中法线的确定和光线回转等问题。最后通过对多种不同的存在解析解的特殊折射率分布的数值解值和解析解值的比较,验证了算法的正确性,并且讨论了应用算法时应该注意的问题和算法中还有待解决的问题。该算法不但具有广泛的适用性,基本不受折射率分布的制约,而且其数值解精度达10-5。该算法在变折射率介质的光学设计和空间成像补偿等方面有着广阔的应用前景。