说明：  采用有限差分法计算旋转对称电场分布,北京理工大学光电学院小学期程序。编了一个多月
(The finite difference method rotationally symmetric field distribution, Beijing Institute of Opto College School of procedures. Made up more than a month)