拉格朗日插值 适合任意个插值点 解多元一次方程组fortran程序.doc_第1页
拉格朗日插值 适合任意个插值点 解多元一次方程组fortran程序.doc_第2页
拉格朗日插值 适合任意个插值点 解多元一次方程组fortran程序.doc_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

!拉格朗日插值 适合任意个插值点的程序.implicit nonereal,dimension(:),allocatable:xreal,dimension(:),allocatable:yreal,dimension(:),allocatable:areal f,x0integer i,j,h,nprint*,请输入插值点的个数read*,nallocate(x(n)allocate(y(n)allocate(a(n)print*,请输入x(i)的值do i=1,nread*,x(i)enddoprint*,请输入y(i)的值do i=1,nread*,y(i)enddoprint*,请输入x的值read*,x0do i=1,n-1a(i)=1h=i+1do j=h,na(i)=a(i)*(x0-x(j)/(x(i)-x(j)enddoenddoa(n)=1do i=2,nh=i-1do j=1,ha(i)=a(i)*(x0-x(j)/(x(i)-x(j)enddoenddo f=0do i=1,nf=y(i)*a(i)+fenddoprint*,f=,fdeallocate(x)deallocate(y)deallocate(a)end!*!解多元一次方程组的程序implicit noneinteger i,j,k,g,e,h,f,nreal,dimension(:,:),allocatable:wreal,dimension(:),allocatable:mreal x,sumprint*,请输入未知数的个数read*,nallocate(w(n,n+1)allocate(m(n+1)print*,请输入方程的增广矩阵do i=1,n do j=1,n+1 read*,w(i,j) enddoenddoprint*,你输入的增广矩阵是do i=1,n print 10,w(i,:) 10 format(2x,6f5.1)enddo do j=1,n-1 e=j+1 f=e do g=e,n if(w(j,j)w(g,j) then do h=j,n+1 m(h)=w(j,h) w(j,h)=w(g,h) w(g,h)=m(h) enddo endif enddo do i=e,n x=w(i,j)/w(j,j) do k=f,n+1 w(i,k)=w(i,k)-x*w(j,k) enddo w(i,j)=0 enddo enddo m(n)=w(n,n+1)/w(n,n) do i=n-1,1,-1 e=i+1 sum=0 do j=n,e,-1 sum=sum+m(j)*w(i,j) enddo m(i)=(w(i,n+1)-sum)/w(i,i) enddo print*,答案是x(i)=. do i=1,n print 10000,m
人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论