programmainimplicitnoneinteger 33603360 I real * 8:3360 x (9)，y )9)，t，z！ //x和y的数据来自matlab的interp1帮助文档，t是插值的节点。 帮助文档变量xqx=[0. d0，0.78539816397448 d0，1.57079632679490 d0，2.35619449019235 d0，3.3参见3.92699081698724 d0，4.71238888888888888 6.28318530717959 d0 ] y=[ 0.d0，0.707106781186548 d0，1.dd0-0.707106781186548 d0，-1.d0，- 0.70710678111186548 d0 file ) ) interp_m.dat ) 102，*tz=0.D0callinterp1(x，y，size(x ) x，t，z ) write(102，* ) t， zenddoclose ) 101 ) close ) 102 ) endprogrammamanage tmp ) implicit none integer : i，j，k1，k2 integer 33603: n real * n //如果插值节点与插值节点相同，则插值节点的值为原始节点的值do i=1，nif(ABS(t-x ) I ) ) 1d-15 ) then tmp=y(i ) return end if end do！ //查找插值的节点所属的区间[X_i，x_I1]I=1doif(tx(I ) ) then exit else i=i 1 end if end do k2=i； K1=K2-1k=(y ) K2-y ) K1 )/x ) K2 )-x ) K1 ) ) tmp=y ) K*(t-x ) k1 ) ) End subroutine interp1