还是找不出错误，这是用秦九韶算法求多项式解简单程序

freedom

[Fortran] 纯文本查看 复制代码

program mytp16
implicit none
integer n
double precision x,poly
double precision,allocatable::cf(:),npoly(:)
write(6,*)'input the nombre n'
read(5,*)n
write(6,*)'input the x'
read(5,*)x
allocate(cf(0:n),npoly(n))
cf=poly(n,npoly,cf,x)
npoly=poly(n,npoly,cf,x)
write(6,'("npoly(",i3,")=",d15.8)')n,npoly(n)
end program

double precision function poly(n,npoly,cf,x)
implicit none
integer i,n
double precision cf(0:n),npoly(0:n),x
open(15,file='the coefficient of npoly.txt')
do i=0,n
cf(i)=2.d0*i+1
write(15,*) cf(i)
end do
close(15)
open(15,file='the coefficient of npoly.txt')
do i=0,n
read(15,*) cf(i)
write(6,'("cf(",i3,")=",d15.8)')i,cf(i)
end do
close(15)
npoly(0)=cf(n)
npoly(1)=cf(n)*x+cf(n-1)
do i=2,n
npoly(i)=npoly(i-1)*x+cf(n-i)
end do
poly=npoly(n)
end function poly
!  --------------------------------------------------
!  Silverfrost FTN95 for Microsoft Visual Studio
!  Free Format FTN95 Source File
!  --------------------------------------------------

fcode

1.请描述你遇到的问题
2.屏幕输入输出，请直接用 * ，而不要用 5 和 6！！！语法并没有规定 5 和 6 的具体含义。
3.你为什么要写入文件再读进来？

freedom

我要用秦九韶算法计算多项式，然后在没有用函数时，运行的是正常的，但是把代码专程函数后总是不出结果，我的作业老师要求这样的，先读入再用

fcode

我让你描述你的问题，而不是描述你的需求。我已经知道你要应用秦九昭算法。
我需要知道的是，你遇到的什么问题？是编译错误？运行错误？还是结果不对？如果遇到错误，请给出错误的现象描述或错误提示。这些都是求助的基本信息。

关于你的代码，文件读写为什么要放在函数里呢？多次调用就会多次写入文件，再读文件。你不觉得浪费吗？还有为什么 poly 要调用两次？
poly 已经是单变量的结果了，就不需要再用一个数组来接收了。

以下是我的修改版，简化了很多东西。

[Fortran] 纯文本查看 复制代码

Program mytp16
  implicit none
  integer n , i
  Real(Kind=8) :: x , y
  Real(Kind=8) , allocatable::cf(:)
  write(*,*)'input the nombre n'
  read(*,*)n
  write(*,*)'input the x'
  read(*,*)x
  allocate( cf(0:n) )
  open(15,file='the coefficient of npoly.txt')
  Do i=0,n
    cf(i) = 2.d0 * i + 1.d0
    write(15,*) cf(i)
  End Do
  close(15)    
  y = poly ( cf , x )
  write(*,'("y=",es13.6)') y

contains

  Real(Kind=8) Function poly( cf , x ) result( y )
    integer i , n
    Real(Kind=8) :: cf(0:) , x
    n = size( cf ) - 1
    y = cf(n) * x + cf( n - 1 )
    Do i = 2 , n
      y = y * x + cf(n-i)
    End Do
  End Function poly

End Program mytp16

freedom

版主太好了，太感谢了！