使用时输入Re和X_d,f为需要求的量
令X_r=34.2/Re^0.875
如果Re<2*10^3,f= 64/Re
  画出图像
如果Re>4*10^3
{
如果 Re>3500/X_d
f=〖[2 log⁡(0.5Xd^(-1) )+1.74]〗^(-2)    画出该公式的函数图像（f-Xd）
如果Re<3500/X_d
   {
    如果X_r>X_d    f^(-(1/2))=2log⁡〖(Ref^(-(1/2)))〗-0.8  画出该公式的函数图像（f-Re）
   如果X_r<X_d   f=〖（2.51/Re-1）〗^2 〖(2 log⁡〖Xd/3.71〗)〗^(-2)    画出这个公式的函数图像（f-Re的函数，也就是Xd当作已知量）
}

}

画图时f和Re坐标都用对数坐标，Re的范围是10^2~10^8
f的范围是0.008~0.1

目前我写的是这样的，未完待续……公式看不清，加个附件吧……