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

工程分析程序设计 上机作业（四）结构化编程上机目的：练习使用选择结构和循环结构编制程序。1、 求y值，x由键盘输入。PROGRAM MAIN IMPLICIT NONE REAL:X,Y WRITE(*,*)请输入X READ(*,*)X IF(X=0.AND.X=10.AND.X=20.AND.X30)THEN Y=X*3+X*2+1 ELSE WRITE(*,*)ERROR END IF WRITE(*,*)YEND PROGRAM2、 输入4个数A，B，C，D，按由大到小的顺序打印出来。PROGRAM MAIN IMPLICIT NONE REAL:A,B,C,D,T WRITE(*,*)请输入A,B,C,D READ(*,*)A,B,C,D IF(AB)THEN T=AA=BB=T END IF IF(BC)THEN T=CC=BB=T END IF IF(CD)THEN T=DD=CC=T END IF IF(AB)THEN T=AA=BB=T END IF IF(BC)THEN T=CC=BB=T END IF IF(AB)THEN T=AA=BB=T END IFWRITE(*,*)A,B,C,DEND PROGRAM3、 编写一个子程序，判断一个整数是否素数。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:Argumentsdatatypesareincompatiblewithintrinsicprocedure,assumeEXTERNAL.DEXP4、 利用上一编程结果，对1000以内的所有偶数验证哥德巴赫猜想。即：对于大于2的任一偶数，先分解为两个奇数之和，然后验证第一个奇数是否素数，如果是，再验证第二个奇数是否素数，如果两个奇数都是素数，则输出结果。（要求给出所有把偶数分解成两个素数之和的等式）。并统计，每个偶数可以分解的等式个数（两个素数互换位置的不算一个独立的结果），研究随着偶数增大，分解的等式个数的变化规律。PROGRAM MIAN IMPLICIT NONE REAL:X,A,B X=6 DO WHILE(X=A) IF(SS(A)=1.AND.SS(B)=1)THEN WRITE(*,*)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一个多，一个少5、 输入一个自然数，进行因子分解并输出结果，例如24=12223。（输出格式不限）。PROGRAM MAIN IMPLICIT NONE REAL:X,Y,I WRITE(*,*)请输入一个自然数 READ(*,*)X I=2 WRITE(*,*)1.000000 DO WHILE(I=X/2) IF(SS(I)=1.AND.MOD(X,I)=0)THEN WRITE(*,*)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