umat中已知结点坐标，把某个属性插值到任意坐标，帮我写段Fortran代码

以下是一个简单的Fortran程序来实现基于节点坐标的插值：

PROGRAM interpolation
  INTEGER, PARAMETER :: n = 10 ! 节点数
  REAL :: x(n), y(n), x_new, y_new ! 节点坐标和要插值的坐标
  REAL :: w(n), w_new ! 插值系数
  INTEGER :: i

  ! 假设已知节点坐标为 (x,y)，属性值为 y
  x = [1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0]
  y = [1.0, 4.0, 9.0, 16.0, 25.0, 36.0, 49.0, 64.0, 81.0, 100.0]

  ! 要插值的坐标为 (x_new,y_new)
  x_new = 4.5
  y_new = 0.0

  ! 计算插值系数
  DO i = 1, n
    w(i) = 1.0
    DO j = 1, n
      IF (i .NE. j) THEN
        w(i) = w(i) * (x_new - x(j)) / (x(i) - x(j))
      END IF
    END DO
    y_new = y_new + w(i) * y(i)
  END DO

  ! 输出插值结果
  WRITE(*,*) 'Interpolated value at (',x_new,',',y_new,') is', y_new

END PROGRAM interpolation

这个程序假设节点的坐标和属性值已经知道，然后通过计算插值系数来得到在要插值的坐标处的属性值。在这个例子中，我们假设要插值的坐标为 (4.5,0.0)。