Fortran循环结构的.pptVIP

第四章 循环结构程序设计;4.1 用DO语句实现循环;DO循环一般格式;P45 例 求5!;P46 例;计算e1、e2、e3的值;1. DO语句只执行一次。如果循环参数表达式e1、e2、e3中含有变量，若在循环体中改变变量的值，但循环参数并不改变。P47例：;2. 不允许在循环体内给循环变量赋值。;3. 退出循环后循环变量的值与最后一次循环时循环变量的值不同，前者比后者多一个步长。;P48 例4.1 求;P48 例4.2 一个整数的因子(不包括该数本身)之和等于它本身，则该数为完数。如6的因子有1,2,3，且1+2+3=6，因此6是完数。输入一个整数，判断它是否完数。;P49 例4.3 Fibonacci数列定义如下： ? F1=1 ? F2=1 ? Fn=Fn-1+Fn-2 (n2) ?求Fibonacci数列的前30项。;P49 例4.4 所谓“水仙花数”是指一个三位整数，其各位数字立方和等于该数本身。例如，153就是一个水仙花数。输出全部“水仙花数”。 ;与循环有关的控制语句;P51 例4.6 求两个整数a与b的最大公约数和最小公倍数。;求a,b最大公约数和最小公倍数,最大公约数*最小公倍数=a*b;CYCLE语句 (continue);P51 例4.7 求1～100之间的全部奇数之和。;4.2 DO WHILE语句（当型循环）;P52 例：输出所输入的全部正数，直到输入负数或零，程序结束。;P53 例4.8 输入一个整数，输出其位数。;integer :: n=1 real :: y=0.0, f do while (y3) f=1.0/(2*n-1) y=y+f n=n+1 end do write(*,*) y=, y-f write(*,*) n=, n-2 end;4.3 几种循环组织方式的比较 ;程序2：用DO WHILE循环实现 IMPLICIT NONE INTEGER M, I, J READ *, M I=2 J=SQRT(REAL(M)) DO WHILE(I=J.AND.MOD(M, I)/=0) I=I+1 END DO IF (IJ.AND.M1) THEN PRINT *,M,’ is a prime number’ ELSE PRINT *,m,’ is not prime number’ END IF END;程序3：用DO循环和逻辑IF语句的嵌套实现。 IMPLICIT NONE INTEGER M, I, J READ *, M I=2 J=SQRT(REAL(M)) DO IF (MOD(M, I)==0 .OR. IJ) EXIT I=I+1 END DO IF ( IJ. AND. M1) THEN PRINT *, M, ’ is a prime number’ ELSE PRINT *, m, ’ is not prime number’ END IF END;4.4 循环的嵌套 ;P59 例4.15 球从100m高度自由落下，落地后弹起。每次弹起的高度都是上次高度的一半。求第10次落地后反弹起的高度和球所经过的路程。 ;INTEGER N REAL H, S, R H=100 S=0 DO N=1, 10 R=H/2 S=S+H+R H=R ENDDO PRINT *, “H=”, H, “ S=”, S END 运行结果如下： H= 9.765625E-02 S= 299.707000;解一元方程; (2) 估计一个根x，给x定一个初值x0，把它代入上式等号的右边，求出x的第一次近似值x1。;例：用迭代法求解x3+2x2+2x+1=0;例4.16 用牛顿迭代法求方程f(x)=2x3-4x2+3x-7=0 在x=2.5附近的实根，直到满足|xn-xn-1|≤10-6 为止。;real :: x=2.5, d=1.0 ! d=f(x)/ f ’(x) do while(abs(d)1.0e-6) d=(((2.0*x-4.0)*x+3.0)*x-7.0)/((6.0*x-8)*x+3.0) x=x-d end do write(*,*) The root is , x end;3.二分法;4.弦截法;real:: x1, x2, f1, f2, f =1.0 do