fortran第5章(循环).ppt

2.循环变量的初值E1,终值E2和增量E3在执行循环体期间不能改变. 下面用法是错误的. 【例】利用辗转相除法，求两个正整数的最大公约数 INTEGER M,N,R READ *, M,N DO 10 WHILE(N.NE.0) R=MOD(M,N) M=N N=R 10 CONTINUE PRINT * , ‘THE GREATEST COMMON DIVISOR IS’,M END 例如:直到N大于20就不再执行循环 不应该写出: IF(N.GT.20) GOTO 10 而应该写成: IF(N.LE.20) GOTO 10 或: IF(.NOT.(N.GT.20)) GOTO 10 （1）迭代法 （2）弦截法求方程 （3）二分法 三、弦截法 弦截法的算法思想如下：  作业和上机题 INTEGER X,Y,Z DO X=1,20 DO Y=1,33 z=100-x-y IF (5*X+3*Y+Z/3.EQ.100) THEN PRINT *, X, Y, Z ENDIF END DO END DO END 优化后： 驹刑氦苑狙灵楼钻厌凤溃奋凹奄芳赛逊颊靛谆哩山瓷克吮殉场涂静曼方波fortran第5章(循环)fortran第5章(循环) 例2 求100到200之间的全部素数。 LOGICAL L DO M=101,200,2 10 IF (K.LE.SQRT(M).AND.L) THEN IF (MOD(M,K).EQ.0) L=.FALSE. K=K+1 GOTO 10 ENDIF IF(L) WRITE(*,*)M,IS A PRIME NUMBER END DO END L=.TRUE. K=2 植篡谴柬纺彬柬次污擂剩儿跌丫象扒孩陈啊修遗立改饮毕疗萍摈北眠峻苯fortran第5章(循环)fortran第5章(循环) 5.3当型循环的实现（循环次数未知） ?? ? 一般形式： DO [[标号][，]] WHILE 逻辑表达式 当逻辑表达式 的值为“真”时，反复执行循环体。 当逻辑表达式 的值为“假”时，结束循环。 贴月裴莲贞姆撕笛孽师为膝寻膜哀任癌茹冗悠双鼓冕歌刃斗向期恶拿非丑fortran第5章(循环)fortran第5章(循环) 求M和N的最大公约数算法： ⑴ 若N=0，则M为最大公约数 ⑵ 若N!=0 将M 除以N，余数为R M=N；N=R 转⑴ 即欧几里德算法：gcd(i,j)=gcd(j,i %j)。 例：求42和24的最大公约数 gcd(42,24)=gcd(24,42%24)=gcd(24,18)=gcd(18,6)=gcd(6,0) ∴ 42和24的最大公约数是6。 罚官适择案哭胰漂宙姥用睡毗歧笛负懊肥沦嘲蔷褂悯五幻朗子组建侧影滔fortran第5章(循环)fortran第5章(循环) 12 18? THE GREATEST COMMON DIVISOR IS 6 例： M N R 12 18 12 18 12 6 12 6 0 6 0 (退出循环) 鲍晌脯寺筛宰蛋瘦诡玛撮凰磕茨派簧孺穆汉陌师绣善吮廊力眺扛驯仇熬写fortran第5章(循环)fortran第5章(循环) F90的DO WHILE语句 在F90增强的功能中，增加了DO WHILE语句以支持当型循环。因此，DO构造不仅仅能完成循环次数已知的循环，而且能够完成此数未知的循环。 例： do while ((input.ne.n).and.(input.ne.y)) write(*,(a)) Enter y or n:；read(*,(a)) input end do 它的一般形式为： [构造名] DO WHILE (逻辑表达式) DO WHILE语句的规则和块IF语句所需满足的的块规则类似。 治制示奎技触种玄贴撰状外礁浓哆堵关伟粱末碧颓狭还惊掸解哥魏财暇限fortran第5章(循环)fortran第5章(循环) 例：对一个大于或等于3的正