《刀具设计手册》上几个问题的探讨

一、hG计算时如何编程（第883页），能在EXCEL中实现吗？如果不能有什么好的方法吗？怎么实现这个迭代运算。

第883页的这个公式好像就更不好解了

第885页的hG＝4.538可能计算有误，我计算的结果为4.859

这个迭代运算我卡这里好久了怎么没人会呀

建议:您可以使用VB平台编程计算,保证可以!

现在这个剃前滚刀的计算就卡在这个迭代运算上，
谁能帮帮我呢
谢谢啦

10-14式子用C语言编程如下，但其中还有点缺陷，望各位多多指教，谢谢！
#include <math.h>
main()
{int i,j;
float r,rG2,hG[5000],alf,rb,A1,a1[5000],a2[5000],a3,a[5000],b[5000],t,m;
scanf("%f,%f,%f,%f,%f",&r,&rG2,&alf,&rb,&A1);
for(i=0;i<=1999;i++)
   {hG<i>=3+i*0.001;
   a1<i>=sqrt(rG2*rG2-(r-hG<i>)*(r-hG<i>))/(r-hG<i>);
   a2<i>=sqrt(rG2*rG2-(r-hG<i>)*(r-hG<i>))/r;
   a3=tan(acos(rb/rG2))-acos(rb/rG2);
   a<i>=(r*(atan(a1<i>)-a2<i>+tan(alf)-alf-a3)+A1/cos(alf))/tan(alf);
   b<i>=fabs(a<i>-hG<i>);
   }
for(j=1;j<=999;j++)
   for(i=1000;i<=2000-j;i++)
      if(b<i>>b[i+1])
    {t=b<i>;b<i>=b[i+1];b[i+1]=t;
     m=hG<i>;hG<i>=hG[i+1];hG[i+1]=m;
     }
printf("%fn",hG[1000]);
}

式子10-16用C语言编程如下，但其中也有点缺陷，望各位多多指教，谢谢！
#include <math.h>
main()
{int i,j;
float rG1[5000],rb,r,hG,lmd1,alf,H,x[5000],y[5000],z[5000],t,m,a1[5000],a2[5000],a;
scanf("%f,%f,%f,%f,%f",&r,&hG,&rb,&alf,&H);a=56;
for(i=0;i<=1999;i++)
   {rG1<i>=a+i*0.001;lmd1=(hG*sin(alf)-H)/r/cos(alf)-tan(alf)+alf;
   a1<i>=sqrt(rG1<i>*rG1<i>-(r-hG)*(r-hG))/(r-hG);
   a2<i>=sqrt(rG1<i>*rG1<i>-(r-hG)*(r-hG))/r;
   x<i>=atan(a1<i>)-a2<i>-lmd1;
   y<i>=acos(rb/rG1<i>);
   z<i>=fabs(tan(y<i>)-y<i>-x<i>);
   }
for(j=1;j<=999;j++)
   for(i=1000;i<=2000-j;i++)
      if(z<i>>z[i+1])
    {t=z<i>;z<i>=z[i+1];z[i+1]=t;
     m=rG1<i>;rG1<i>=rG1[i+1];rG1[i+1]=m;
     }
printf("%fn",rG1[1000]);
}

用上面的C语言程序计算出的P885页上的hG＝4.859
rG1＝57.388
αF＝16°58′6″，书上αF＝15°19′49″计算得不准确，但最后取αF＝14°是肯定满足要求的。

这个公式只是一个近似公式，只在A1趋于0时才精确，主要是过渡曲线并不是刀头一点的长幅渐开线形成，而是刀头圆心的长幅渐开线的等距线形成的。