拉格朗日插值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）、計算結果


程序清單