Fortran Coder

查看: 148|回复: 1

[求助] 无报错运行错误-拉格朗日多项式插值算法

[复制链接]

1

帖子

1

主题

0

精华

新人

F 币
18 元
贡献
9 点
发表于 2021-10-8 19:15:14 | 显示全部楼层 |阅读模式
1F 币
以下为代码,捣鼓一下午了
[Fortran] 纯文本查看 复制代码
Program lagrange
Implicit none
integer::t,n,i,j,k
real::tmp
real,dimension(1)::x,y
real,dimension(2)::p(17,2)
print*,"请输入气压值"
read*,t
open(10,file='d:\my fortran\air_tmp.txt')
read(10,*)((p(k,j),j=1,2),k=1,17)!17行2列
close(10)
do i=1,17
x(i)=p(i,1)
y(i)=p(i,2)
enddo
print*,x,y
tmp=fun(t)
print*,tmp

contains
function fun(t) result(tmp)
implicit none
integer::t,n,i,j
real::tmp
real,dimension(1)::x,y,l
do i=1,n
    do j=1,n
        if(i/=j) then
        l(i)=l(i)*(t-x(j))/(x(i)-x(j))
    else
        return
    endif
    enddo
tmp=l(i)*y(i)+tmp
enddo
endfunction fun

end program

回复

使用道具 举报

1637

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
2047 元
贡献
1292 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2021-10-9 08:56:58 | 显示全部楼层
你的主要问题是:
1. 累加和累乘,没用赋初值。累加应该赋值0,累乘应该赋值1
2. contains下面的子程序中使用主程序的变量,不要重新定义,否则会覆盖主程序中变量的作用。
3. 数组定义混乱。当前的Fortran,基本上很少用到 dimension。(除非同时定义多个大小相同的数组)

[Fortran] 纯文本查看 复制代码
Program lagrange
  Implicit none
  Integer , parameter :: S = 17
  real :: tmp , t , x(S) , y(S)
  integer :: i
  write(*,*) "请输入气压值"
  read(*,*) t
  open(10,file='a.txt')
  do i = 1 , size(x)
    read(10,*) x(i),y(i)
    write(*,*) x(i),y(i)
  end do
  close(10)
  tmp=fun(t)
  write(*,*) tmp

contains

  function fun(t) result(tmp)
    integer :: i , j
    real :: tmp , t , l
    tmp  = 0.0
    do i = 1 , size(x)
      l = 1.0
      do j = 1 , size(x)
        if( i /= j ) then
          l = l * (t-x(j))/(x(i)-x(j))
        end if
      end do
      tmp= l * y(i) + tmp
    end do
  endfunction fun

end program lagrange
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2021-10-24 08:58

Powered by Tencent X3.4

© 2014-2021 Tencent

快速回复 返回顶部 返回列表