编写一个fortran90程序，利用契贝谢夫多项式（取前十项）计算gamma函数的近似值

以下是一个使用Fortran 90编写的程序，利用契贝谢夫多项式来计算gamma函数的近似值（取前十项）：

program gamma_approximation
  implicit none
  integer, parameter :: k = 10 ! 取前十项
  real(8) :: x, gamma_approx
  real(8) :: chebyshev_poly(k)
  
  ! 契贝谢夫多项式系数
  data chebyshev_poly /1.00000000, -0.57721566, 0.42278434, -0.08423902, 0.02604047, -0.01181581, &
                       0.00571722, -0.00290893, 0.00191876, -0.00126249/
  
  ! 输入x值
  write(*, *) "请输入要计算gamma函数的值（x）："
  read(*, *) x
  
  ! 计算契贝谢夫多项式的和
  gamma_approx = chebyshev_poly(1)
  do i = 2, k
    gamma_approx = gamma_approx + chebyshev_poly(i) * cos((i-1) * acos(x))
  end do
  
  ! 计算gamma函数的近似值
  gamma_approx = gamma_approx * (1 - x) + log(x)
  
  ! 输出结果
  write(*, *) "gamma(", x, ") 的近似值为：", gamma_approx
  
end program gamma_approximation

请注意，这只是一个简单的近似计算程序，并不保证精确性。在实际计算中，可能需要更复杂的算法和更多的项数来获得更准确的结果。