5.1
根据题目条件和函数方程，就可以调用已经写好的 N-R 法函数。
先画出y=u0*N^2*pi*a^2*((l^2+a^2)^(1/2)-a)/l^2-L的图像，可
以看出图像和直线y=0的交点就是函数的解，大概在0.04左右。函数图像如下
图所示：
用 N-R 法算出的方程的解的结果如下图所示：
5.6
根据题目条件和函数方程，就可以调用已经写好的 N-R 法函数。
Z 的变化范围是 0.001~0.1cm，分辨率为 0.01cm，可以得到横轴的变化范围
和变化区间，z 的取值为 0.001，0.011，0.021，0.031，0.041，0.051，0.061，
0.071，0.081，0.091，每取一个 z 值，计算一次 rcru，在图中有十个点，连接的
曲线就是临界距离 rcru 随 z 的变化曲线，如下图所示：
