拉格朗和牛顿插值的fortran程序

wangyu0824 · 发表于 2016-4-9 22:47:59

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

拉格朗和牛顿插值的fortran程序，自己写的总是不对，求解，还出现NaN。
function Lagrange(n,x,y,xx)
real(8):: x(5),y(5), xx
real :: sum = 0, l = 1, ch, i, j
do i = 1, n + 1
do j = 1, n + 1
if (j /= i) then
l = l * (xx - x(j)) / (x(i) - x(j))
end if
end do
ch = l * y(i)
sum = sum + ch
end do
Lagrange = sum

end function Lagrange

wangyu0824 · 发表于 2016-4-9 22:48:35

求解，非常感谢

wangyu0824 · 发表于 2016-4-9 22:49:44

牛顿插值自己也是很晕，求一个完整的程序

wangyu0824 · 发表于 2016-4-9 22:50:54

俺是新手上路

聆风 · 发表于 2016-4-10 09:32:41

这是我们上课时的程序，希望对你有用

wangyu0824 · 发表于 2016-4-10 13:29:12

好的，我看看，谢谢

wangyu0824 · 发表于 2016-4-10 13:30:41

这个程序我也有，我就是想知道自己错在哪里了？

wodedipan · 发表于 2016-4-10 13:44:22

没看懂啥意思啊

robindct · 发表于 2016-4-11 09:11:21

你的程序只是子函数，我不太确定，你定义的x和y容量都是5，倒是i和j循环都是到n+1，你就要确定一下n<5的问题了~~

hab · 发表于 2016-10-9 16:09:18

聆风发表于 2016-4-10 09:32
这是我们上课时的程序，希望对你有用

多谢，刚开始写，好蒙阿。。。

		自动登录	找回密码
密码			立即注册

[求助] 拉格朗和牛顿插值的fortran程序

登录后查看更多精彩内容~