如果您无法下载资料,请参考说明:
1、部分资料下载需要金币,请确保您的账户上有足够的金币
2、已购买过的文档,再次下载不重复扣费
3、资料包下载后请先用软件解压,在使用对应软件打开
一、FORTRAN语言的三点说明1.常数与变量的说明在FORTRAN77中有字符型变量,字符型常数2.语句书写的格式3.与C语言的不同二、FORTRAN基本语句⑥calls(d1,d2,…,dn)②算数条件语句:if(e)k1,k2,k3(e<0)(e=0)(e>0)算数条件语句的计算编程正确表示:if(X.le.2.1)Y=0.5*X+0.95if(X.gt.2.1)Y=0.7*X+0.53write(*,*)X,Ycha1-3.f:read(*,20)x,y,z20format(3F10.4)big=xif(y.gt.big)big=yif(z.gt.big)big=zwrite(*,*)’big=’,bigend④循环do语句例4:编程求解0.0,0.1,0.2,0.3的平方根do10i=0,0.3,0.1⑥calls(d1,d2,…,dn)return(在end前)(3)输入、输出语句2.非执行语句(3)data语句(数据初值语句)datav1/d1/,v2/d2/,…vn/dn/或datav1,v2,…,vn/d1,d2,…,dn/函数子程序function举例:integerp,rwrite(*,*)'n=,r=?'read(*,*)n,rcallfac(n,m)m1=mcallfac(n-r,m)m2=mp=m1/m2write(*,*)n,r,pendsubroutinefac(i,m)m=1if(i.le.1)goto77do10k=2,i10m=m*k77returnendopen(1,file=’N!.dat’)write(*,*)’inputN=?’read(*,*)NM=1I=25M=M*II=I+1if(I.gt.N)goto10goto510write(1,*)Mend三、源程序语句排列顺序Aw/四、FORTRAN常用内部函数(库函数)错误正确a(-b)a*-ba*(-b)(ab)3a*b**3(a*b)**3sin6tsin6tsin(6*t)aexa*e**xa*exp(x)6.8log26.56.8*log26.56.8*alog10(26.5)1.“/”及“*”不能省略;一律用小括号。2.同类型算术量间才能运算。五、有关循环语句implicitreal*8(A-H,O-Z)open(1,file=’sin.dat’)write(*,*)’inputw,N’read(*,*)w,Npi=3.1415926do10I=1,Nt=2.*pi/wt=t*float(I)/Nx=sin(w*t)write(*,*)t,x10write(1,*)t,xend(指数增加每次加1,可以用i代替)a=1do10i=1,6310a=a+2**iwrite(*,*)aendwrite(*,*)’inputx,N’read(*,*)x,Ny=0do20M=1,Nt=1s=xdo100J=1,Mt=t*float(J)100s=s*xy=y+s/t20write(*,*)M,yend六、FORTRAN语言的优缺点本课程后续章节中,有关例题均采用Fortran编译,如若用其他算法语言编程,Fortran源码仅供参考。EX1-1:编程计算,直至大于为止。