西安交大fortran上机作业之四.docxVIP

四工程分析程序设计上机作业（四）结构化编程上机目的：练习使用选择结构和循环结构编制程序。求y值，x由键盘输入。PROGRAM MAINIMPLICIT NONE REAL X,Y PRINT*,请输入X值 READ*,X IF (X=0.AND.X10) THEN Y=XELSEIF (X=10.AND.X20) THEN Y=X**2+1ELSEIF (X=20.AND.X30) THEN Y=X**3+X**2+1ELSE PRINT*,X不在取值范围内END IFPRINT*,Y值为：,YEND输入4个数A，B，C，D，按由大到小的顺序打印出来。PROGRAM MAINIMPLICIT NONEREAL A,B,C,D,T REAL::E(4) INTEGER I,J PRINT*,请输入四个数： READ*,A,B,C,DE(1)=AE(2)=BE(3)=CE(4)=D DO I=1,3 DO J=I+1,4 IF (E(I)E(J)) THEN T=E(I)E(I)=E(J)E(J)=T END IFEND DO END DO PRINT*,E END编写一个子程序，判断一个整数是否素数。PROGRAM MAINIMPLICIT NONEINTEGER X,SSPRINT*,请输入一个整数：READ*,X CALL SUSHU(X,SS)CONTAINSSUBROUTINE SUSHU(A,SS) IMPLICIT NONE INTEGER I,A,B,SS,K K=0 DO I=2,A-1 B=MOD(A,I) IF (B==0) THEN K=K+1END IF END DO IF (K==0) THEN PRINT*,是素数 ELSE PRINT*,不是素数 END IFEND SUBROUTINEEND利用上一编程结果，对1000以内的所有偶数验证哥德巴赫猜想。即：对于大于2的任一偶数，先分解为两个奇数之和，然后验证第一个奇数是否素数，如果是，再验证第二个奇数是否素数，如果两个奇数都是素数，则输出结果。（要求给出所有把偶数分解成两个素数之和的等式）。并统计，每个偶数可以分解的等式个数（两个素数互换位置的不算一个独立的结果），研究随着偶数增大，分解的等式个数的变化规律。PROGRAM MAINIMPLICIT NONEINTEGER I,x,y,SS,NDO I=2,1000,2N=0 DO X=1,I/2,2 Y=I-X CALL SUSHU(X,SS) IF (SS==0) THEN CALL SUSHU(Y,SS) IF (SS==0) THEN PRINT*,I,X,Y N=N+1 END IF END IF END DO PRINT*,共有N种,NEND DOCONTAINSSUBROUTINE SUSHU(A,SS) INTEGER J,A,B,SS,K K=0 DO J=2,A-1 B=MOD(A,J) IF (B==0) THEN K=K+1END IF END DO IF (K==0) THEN SS=0 ELSE SS=1 END IFEND SUBROUTINEEND输入一个自然数，进行因子分解并输出结果，例如24=1×2×2×2×3。（输出格式不限）。PROGRAM MAINIMPLICIT NONEINTEGER X,I,A,K,B,JPRINT*,请输入一个自然数:READ*,XB=XDO I=1,B DO J=2,X A=MOD(X,J) IF (A==0) THEN X=X/J PRINT*,J EXIT END IFEND DOEND DOEND