FORTRAN数值方法及其在物理学中应用课件51708.pptVIP

数值方法及其在物理学中应用 一、FORTRAN语言的三点说明1. 常数与变量的说明 2. 语句书写的格式 FORTRAN77扩展名为: *.f或*.for FORTRAN90无严格限制（以.f90为后缀） 二、FORTRAN基本语句 1、可执行语句（2）流程控制语句 ⑥call sd1,d2,…,dn②算数条件语句: ife k1, k2, k3e0 e0 e0 算数条件语句的计算编程 cha1-1.f:read*,40 X40formatF8.2ifx10,20,3010Y-1.57079goto 10020Y0goto 10030Y1.57079goto 100100write*,50x,Y50format1X,2HX,F10.6,4H, YF10.6end ③逻辑条件语句 cha1-3.f:read*,20x,y,z 20format3F10.4bigxify.gt.bigbigyifz.gt.bigbigzwrite*,*’big’,bigend 例4：编程求解0.0, 0.1, 0.2, 0.3的平方根do 10 i0, 0.3, 0.1 2. 非执行语句 3data语句数据初值语句data v1/d1/,v2/d2/,…vn/dn/或 data v1,v2,…,vn/d1,d2,…,dn/ 函数子程序 function举例:integer fac,p,rwrite*,*n,r?read*,*n,rpfacn/facn-rwrite*,*n,r,pendinteger function facnfac1ifn.le.1 goto 77do 10 k2,n 10facfac*k 77returnendinteger p,rwrite*,*n,r?read*,*n,rcall facn,mm1mcall facn-r,mm2mpm1/m2write*,*n,r,pendsubroutine faci,mm1ifi.le.1 goto 77do 10 k2,i 10 mm*k 77 returnendopen1,file’N!.dat’write*,*’input N?’read*,*NM1I2 5MM*III+1ifI.gt.N goto 10goto 5 10write1,* Mend 三、源程序语句排列顺序四、FORTRAN常用内部函数（库函数）错误正确 a-ba*-ba*-bab3a*b**3a*b**3 sin6tsin6tsin6*taexa*e**xa*expx 6.8log26.5 6.8*log26.56.8*alog1026.5 1. “/”及“*”不能省略；一律用小括号。 2．同类型算术量间才能运算。 五、有关循环语句 （1）循环变量最好不要第二次赋值。 （2）可以从循环体内转向循环体外，不允许循环体外转向循环体内。 （3）注意多重循环问题。implicit real*8A-H,O-Zopen1,file’sin.dat’write*,*’input w,N’read*,*w,Npi3.1415926do 10 I1,Nt2.*pi/wtt*floatI/Nxsinw*twrite*,*t, x 10write1,*t, xend （指数增加每次加1，可以用i代替）a1do 10 i1, 63 10aa+2**iwrite*,* aendwrite*,*’input x, N’read*,*x,Ny0do 20 M1,Nt1sxdo 100 J1,Mtt*floatJ 100ss*xyy+s/t 20write*,* M,yend§1.2 MATLAB 语言简介 观看教学录像片read*,* tdo 10 K1,5dt10.**1-Kv19.*t+4.5*dt-6.*t*t-6.*t*dt-2.*dt*dta19.-12.*t-6.*dt 10write*,*dt,v1,a1v9.*t-6.*t*ta9.-12.*twrite*,*v,aend 二、运动方程问题弦截法 弦截法 弦截法 弦截法read*,*a,b,v,efaa+alog1.+a-vfbb+alog1.+b-viffa*fb.gt.0 goto 1000ifabsfa.gt.e goto 5cagoto 100 5ifabsfb.gt.e goto 10cbgoto 100 10cb-a-b/fa-fb*fbfcc+alog1.+c-vifabsfc.le.e goto 100iffa*fc.gt.0 goto 15fbfcbc!c与b同号goto 10 15fafcac!c与a同号goto 10 100write*,* ’c’, cgoto 111 1