程序代码：

#include <stdio.h>
#include <math.h>

#define pi 3.1415926

void main()
{
  double a=5.0755;
  double b=10.12658;
  double fh;
  double q1=pow((b+1)/(b-1),0.5);
  double q2=pow(a*a+(b+1)*(b+1),0.5);
  double q3=pow(a*a+(b-1)*(b-1),0.5);
  double q4=pow((a*a+(b+1)*(b+1))/(a*a+(b-1)*(b-1)),0.5);
  double q5=pow((b-1)/(b+1),0.5);
  fh=(atan(q1)-(a*a+b*b-1)/q2/q3*atan(q4*q5))/pi;
  printf("=%lf\n",fh);
}

[此贴子已经被作者于2017-10-11 20:25编辑过]