工程分析程序设计 上机作业(四).docx

工程分析程序设计上机作业（四）结构化编程上机目的：练习使用选择结构和循环结构编制程序。求y值，x由键盘输入。PROGRAM MAIN IMPLICIT NONE REAL::X,Y WRITE(*,*)请输入XREAD(*,*)X IF(X=0.AND.X10)THEN Y=X ELSE IF(X=10.AND.X20)THEN Y=X**2+1 ELSE IF(X=20.AND.X30)THEN Y=X**3+X**2+1 ELSEWRITE(*,*)ERROR END IF WRITE(*,*)YEND PROGRAM输入4个数A，B，C，D，按由大到小的顺序打印出来。PROGRAM MAIN IMPLICIT NONE REAL::A,B,C,D,T WRITE(*,*)请输入A,B,C,DREAD(*,*)A,B,C,DIF(AB)THEN T=AA=BB=T END IFIF(BC)THEN T=CC=BB=T END IF IF(CD)THEN T=DD=CC=T END IFIF(AB)THEN T=AA=BB=T END IFIF(BC)THEN T=CC=BB=T END IFIF(AB)THEN T=AA=BB=T END IFWRITE(*,*)A,B,C,DEND PROGRAM编写一个子程序，判断一个整数是否素数。PROGRAM MAIN IMPLICIT NONE REAL::X,I REAL::Y I=2 WRITE(*,*)请输入一个整数READ(*,*)X Y=SQRT(X) DO WHILE(I=Y)IF(MOD(X,I)==0)THEN WRITE(*,*)该整数不是素数 STOP END IFI=I+1 END DO WRITE(*,*)该整数是素数 END PROGRAM注：在此题中，若定义X为整数则无法运行，Warning:?Arguments?data?types?are?incompatible?with?intrinsic?procedure,?assume?EXTERNAL.???[DEXP]利用上一编程结果，对1000以内的所有偶数验证哥德巴赫猜想。即：对于大于2的任一偶数，先分解为两个奇数之和，然后验证第一个奇数是否素数，如果是，再验证第二个奇数是否素数，如果两个奇数都是素数，则输出结果。（要求给出所有把偶数分解成两个素数之和的等式）。并统计，每个偶数可以分解的等式个数（两个素数互换位置的不算一个独立的结果），研究随着偶数增大，分解的等式个数的变化规律。PROGRAM MIAN IMPLICIT NONE REAL::X,A,B X=6 DO WHILE(X=1000)WRITE(*,*)X,= A=3 B=X-ADO WHILE(B=A)IF(SS(A)==1.AND.SS(B)==1)THENWRITE(*,*)A,+,B END IF A=A+2 B=B-2END DOX=X+2 END DO CONTAINS FUNCTION SS(X) REAL X,I,Y REAL SS I=2 Y=SQRT(X) DO WHILE(I=Y)IF(MOD(X,I)==0)THEN SS=0 GO TO 100 END IF I=I+1 END DO SS=1 100END FUNCTIONEND PROGRAM一个多，一个少输入一个自然数，进行因子分解并输出结果，例如24=1×2×2×2×3。（输出格式不限）。PROGRAM MAIN IMPLICIT NONE REAL::X,Y,I WRITE(*,*)请输入一个自然数READ(*,*)X I=2WRITE(*,*)1.000000 DO WHILE(I=X/2)IF(SS(I)==1.AND.MOD(X,I)==0)THENWRITE(*,*)I Y=X/I DO WHILE(MOD(Y,I)==0)WRITE(*,*)I Y=Y/I END DOEND IF I=I+1 END DO CONTAINS FUNCTION SS(X) REAL X,I,Y REAL SS I=2 Y=SQRT(X) DO WHILE(I=Y)IF(MOD(X,I)==0)THEN SS=0 GO TO 100 END IF I=I+1 END DO SS=1 100END FUNCTIONEND