FORTRAN数值方法技巧及其在物理学中应用.ppt

数值方法及其在物理学中应用;课 程 提 纲 ;第三章 物理学中定积分的数值计算方法 §3-1 定积分基本数值算法及其应用 §3-2 龙贝格法及其应用 §3-3 高斯求积法?;第六章 物理学中非线性方程的求根问题 §6-1 根的有哪些信誉好的足球投注网站和二分法 §6-2 函数迭代法 §6-3 牛顿迭代法?;课 程 特 点 ;课时：46学时＋24机时; 第一章 算法语言简介与误差分析初步;一、FORTRAN语言的三点说明 ; 在FORTRAN 77中有字符型变量，字符型常数;2. 语句书写的格式;3. 与C语言的不同;二、FORTRAN基本语句;⑥call s(d1,d2,…,dn) ;②算数条件语句: if(e) k1, k2, k3 (e0) (e=0) (e0);算数条件语句的计算编程; ;正确表示: if(X.le.2.1) Y=0.5*X+0.95 if(X.gt.2.1) Y=0.7*X+0.53 write(*,*)X,Y;cha1-3.f: read(*,20)x,y,z 20 format(3F10.4) big=x if(y.gt.big)big=y if(z.gt.big)big=z write(*,*)’big=’,big end;④ 循环do语句; 例4：编程求解0.0, 0.1, 0.2, 0.3的平方根 do 10 i=0, 0.3, 0.1; ⑥call s(d1,d2,…,dn) return (在end前);(3) 输入、输出语句;2. 非执行语句;(3)data语句(数据初值语句) data v1/d1/,v2/d2/,…vn/dn/ 或 data v1,v2,…,vn/d1,d2,…,dn/;函数子程序 function举例:; integer p,r write(*,*)n=,r=? read(*,*)n,r call fac(n,m) m1=m call fac(n-r,m) m2=m p=m1/m2 write(*,*)n,r,p end ? subroutine fac(i,m) m=1 if(i.le.1) goto 77 do 10 k=2,i 10 m=m*k 77 return end ; open(1,file=’N!.dat’) write(*,*)’input N=?’ read(*,*)N M=1 I=2 5??? M=M*I I=I+1 if(I.gt.N) goto 10 goto 5 10 write(1,*) M end;三、源程序语句排列顺序;Aw;/;四、FORTRAN常用内部函数（库函数）; 错误 正确 a(-b) a*-b a*(-b) (ab)3 a*b**3 (a*b)**3 sin6t sin6t sin(6*t) aex a*e**x a*exp(x) 6.8log26.5 6.8*log26.5 6.8*alog10(26.5) ; 1. “/”及“*”不能省略；一律用小括号。 2???同类型算术量间才能运算。 ;五、有关循环语句; implicit real*8(A-H,O-Z) open(1,file=’sin.dat’) write(*,*)’input w,N’ read(*,*)w,N pi=3.1415926