编写牛顿均差插值多项式的程序用fortran语言
来源:学生作业帮 编辑:搜搜做题作业网作业帮 分类:综合作业 时间:2024/07/10 16:50:59
编写牛顿均差插值多项式的程序用fortran语言
![编写牛顿均差插值多项式的程序用fortran语言](/uploads/image/z/5596958-38-8.jpg?t=%E7%BC%96%E5%86%99%E7%89%9B%E9%A1%BF%E5%9D%87%E5%B7%AE%E6%8F%92%E5%80%BC%E5%A4%9A%E9%A1%B9%E5%BC%8F%E7%9A%84%E7%A8%8B%E5%BA%8F%E7%94%A8fortran%E8%AF%AD%E8%A8%80)
去年学的Fortran语言,许多算法都忘了.给你一个去年写的牛顿插值程序,需要改动的地方自己改,貌似循环计算(多层计算)那块处理得不错.
real*8 x(0:5),y(0:5),c(0:5),u1,u2,ci
data y/0.796,0.773,0.744,0.704,0.656,0.602/
open(1,file='7-2.dat')
write(*,*)'input u1=?,u2=?'
read(*,*)u1,u2
do i=0,5
x(i)=(i+1)*0.125
enddo
c(0)=y(0)
do i=1,5
call chashang(x,5,i,ci)
c(i)=ci
write(*,*)'c(i)=,i=',c(i),i
write(1,*)'c(i)=,i=',c(i),i
enddo
call newchazhi(u1,c,x,5)
call newchazhi(u2,c,x,5)
end
subroutine newchazhi(u,c,x,m)
real*8 c(0:m),x(0:m),u,N
S=c(0)
do i=1,5
N=c(i)
do j=0,i-1
N=N*(u-x(j))
enddo
S=S+N
enddo
write(*,*)u,s
write(1,*)u,s
end
subroutine chashang(x,m,i,ci)
real*8 y(0:5),x(0:m),f(1:i),ci
y(0)=0.796
y(1)=0.773
y(2)=0.744
y(3)=0.704
y(4)=0.656
y(5)=0.602
k=0
10 do j=1,i-k
f(j)=(y(j)-y(j-1))/(x(j+k)-x(j-1))
enddo
if(j==1) then
goto 5
else
do j=1,i-k
y(j-1)=f(j)
enddo
k=k+1
goto 10
endif
5 ci=f(j)
return
end
real*8 x(0:5),y(0:5),c(0:5),u1,u2,ci
data y/0.796,0.773,0.744,0.704,0.656,0.602/
open(1,file='7-2.dat')
write(*,*)'input u1=?,u2=?'
read(*,*)u1,u2
do i=0,5
x(i)=(i+1)*0.125
enddo
c(0)=y(0)
do i=1,5
call chashang(x,5,i,ci)
c(i)=ci
write(*,*)'c(i)=,i=',c(i),i
write(1,*)'c(i)=,i=',c(i),i
enddo
call newchazhi(u1,c,x,5)
call newchazhi(u2,c,x,5)
end
subroutine newchazhi(u,c,x,m)
real*8 c(0:m),x(0:m),u,N
S=c(0)
do i=1,5
N=c(i)
do j=0,i-1
N=N*(u-x(j))
enddo
S=S+N
enddo
write(*,*)u,s
write(1,*)u,s
end
subroutine chashang(x,m,i,ci)
real*8 y(0:5),x(0:m),f(1:i),ci
y(0)=0.796
y(1)=0.773
y(2)=0.744
y(3)=0.704
y(4)=0.656
y(5)=0.602
k=0
10 do j=1,i-k
f(j)=(y(j)-y(j-1))/(x(j+k)-x(j-1))
enddo
if(j==1) then
goto 5
else
do j=1,i-k
y(j-1)=f(j)
enddo
k=k+1
goto 10
endif
5 ci=f(j)
return
end
编写牛顿均差插值多项式的程序用fortran语言
用牛顿切线法解决一个问题(方程)用fortran语言编写程序
请用fortran语言帮忙编写一下程序
用Fortran语言编写两个4*4矩阵乘的程序
求教 fortran 高手! 怎么用fortran语言编写程序 用来构建一个三维网格
fortran程序 第一次接触fortran语言,试着编写了一个简单的程序:program tsimplicit non
c语言:用递归方法编写程序,求n阶勒让德多项式的值
Fortran语言小程序编写,请各路大神帮忙
用c语言编写程序计算5!+10!的值
如何知道自己编写的fortran程序哪里错了
用c语言编写一个程序,计算多项式2x^3-5x^2+6的值,其中x=2.55,怎么编啊
用fortran语言编一个矩阵相乘的通用程序 急用