拉格朗日插值c语言

A. 拉格朗日插值算法

全区间拉格朗日插值

功用
本程序用拉格朗日插值公式对一元不等距观测数据进行程组插值 。

方法概要
对给定的n个插值节点x1,x2,…,xn及对应的函数值y1,y2,…,yn,计算给定点x的函数值y(x)。

本程序可以在插值区间内对给定的NJ个插值点进行插值。

程序说明
(1)、程序名

LGRCZ1.FOR (此文件为免费软件，你可以从本站点 下载。)

(2)、子程序语句

SUBROUTINE LGRCZ1（Z,N,CZ,NJ）

(3)、形参说明

Z 输入参数，两维实数组，容量为N×2，存放给定的插值节点数列及对应的观测值。
N输入参数，整变量，插值节点数。
NJ 输入参数，整变量，插值节点数。
CZ 输入输出参数，二维实数组，容量NJ×2，第一列输入插值点数列，第二列输出插值

结果。

计算实例
（1）、例题

设已知下表观测数据

x
1
2
3
4
5
6
7

Y1
1
4
9
16
25
36
49

要在给出2,2.5,3,3.5,4,4.5,5,5.5,6,6.5上进行插值。这里N＝7，NJ＝10。

（2）、主程序设计

N、NJ用参数说明语句定义，观测数据及插值点数用DATA语句输入，屏幕显示插值结果。

（3）、计算结果


程序清单