实验报告6
1.子例行子程序ZZ的功能是将N行M列的整型矩阵A转置。调用该子程序，求一个3行4列的整型矩阵X的转置矩阵Y并打印。请在画线处填空,以完成题目的要求。 PROGRAM EXAM IMPLICIT NONE
INTEGER,DIMENSION(3,4)::X,Y(4,3) INTEGER::I,J
DO I=1,3; READ (*,'(4I4)') (X(I,J),J=1,4); END DO PRINT '(1X,4I4)',((X(I,J),J=1,4),I=1,3) CALL ZZ(X,3，4，Y) PRINT * ,'输出转置矩阵：'
WRITE(*,'(1X,3I4)')((Y(I,J),J=1,3),I=1,4) END PROGRAM EXAM SUBROUTINE ZZ(A,N,M,B) IMPLICIT NONE INTEGER::N,M,I,J
INTEGER,DIMENSION(:,:)::A(N,M),B(M,N) DO I=1,N
DO _J=1，M_____ _ ____ B(J,I)=A(I,J) _______ ____
END DO END DO RETURN
END SUBROUTINE
⑴ 编写一个子例行子程序IAC用来求N！(值为实型)。子程序中应对自变量进行检查，对不合理的变量给出错误信息。然后求下式(R和K的值由键盘输入)的值。 C?R! K!(R?K)!
FUNCTION IA(CN) RESULT(IAC-RES) INTEGER::N,IAC-RES,I IF(N<=0) THEN PRINT*,”ERROR” ELSE
IAC-RES=1 DO I=2，N
IAC-RES=IAC-RES*I END DO END IF
END FUNCTION IAC
INTEGER::R,K
DO WHILE(R<0,OR,K<0,OR,R C=1.0*IAC(R)/(IAC(K)*IAC(R-K)) PRINT*,C END PROGRAM EXAM 若输入8,6. 则运行结果为：28.000000 ⑵ 求s= ?n!=1!+2!+......+10! n?110PROGRAM EXAM IMPLICIT NONE INTEGER::S=0,T,N,J DO N=1，10 T=1 DO J=1,K T=T*J END DO S=S+T END DO PRINT ‘(1X,”S=”,I10)’,S END PROGRAM EXAM 运行的结果为S= 4037913 七数组及应用 1.下面程序是用折半查找法在给定的一维有序（假设是从小到大）整型数组中查找给定的元素K。若找到，则返回该元素在数组中的下标，若数组中没有该元素，则返回0。请在画线处填空,以完成题目的要求。 PROGRAM EXAM IMPLICIT NONE INTEGER, DIMENSION(10)::X INTEGER ::K,LOW,HIGH,MID READ *,X; PRINT *,'请输入要查找的K：' ; READ*,K LOW=1;HIGH=10 DO WHILE(LOW<=HIGH) MID=(LOW+HIGH)/2 IF(______ __) EXIT IF(X(MID)>K) THEN HIGH=MID-1 ELSE ____ ____ END IF END DO IF(________) THEN PRINT *,0 ELSE PRINT *,MID END IF END 2.编写程序: ⑴输入10个整型数据，然后任选一种方法进行排序。 ⑵ 设数组A(3,4)，编程实现下列操作： ①找出最大元素及其所在的行与列； ②找出第K行的的最小元素及其所在的列，并求出第K行元素之和； ③交换第I列与第J列的元素。 ⑶ 编写一个程序，打印杨辉三角形的前10行。 1 1 1 1 2 1 1 3 3 1 1 4 6 4 1 ?? 1,三个空:(1)X(MID)==K (2)LOW=MID+1 (3)LOW>HIGH 2,编写程序 <1>PROGRAM EXAM IMPLICIT NONE INTEGER, DIMENSION(10)::A INTEGER::I,J,K,T PRINT*,”请输入数组 A:” READ*,A DO I=1,9 K=I DO J=I+1,10 IF(A(J) T=A(I);A(I)=A(K);A(K)=T END DO PRINT*,” 排序后的数组为：” PRINT ‘(1X,10I4)’,A END PROGRAM EXAM 请输入数组 A: 3,9,-3,0,15,19,4,8,1,2 排序后的数组为： -3 0 1 2 3 4 8 9 15 19 ⑵ 设数组A(3,4)，编程实现下列操作： ①找出最大元素及其所在的行与列； ②找出第K行的的最小元素及其所在的列，并求出第K行元素之和； ③交换第I列与第J列的元素。 <2>1. PROGRAM EXAM IMPLICIT NONE INTEGER, DIMENSION(3,4)::A INTEGER::I,J,MAX,ROW=1,COLV=1 DO I=1,3 READ*,(A(I,J),J=1,4) END DO MAX=A(1,1) DO I=1,3 DO J=1,4 IF(A(I,J)>MAX) THEN MAX= A(I,J); ROW=I; COLV=J END IF END DO END DO PRINT*,MAX,ROW,COLV END PROGRAM EXAM 2. PROGRAM EXAM IMPLICIT NONE INTEGER, DIMENSION(3,4)::A INTEGER::I,J,SUM=0,K,MIN, COLV=1 DO I=1,3 READ*,(A(I,J),J=1,4) END DO READ*,K MIN=A(K,1) DO J=1,4 SUM= SUM+A(K,J) IF(A(K,J)>MIN) THEN MIN=A(K,J) ; COLV=J END IF END DO PRINT ‘(1X,”最小元素为”,I4,2X,”所在的列号为”,I3)’, MIN, COLV PRINT ‘(1X,”…………….”,I3,”……………………”,I6)’,K,SUM END PROGRAM EXAM 3. PROGRAM EXAM IMPLICIT NONE INTEGER, DIMENSION(10,10)::A INTEGER::I,J DO I=1,10 DO I=1,I IF((J=-1)OR(I==J) THEN !!!!!!!!!!!!!!!!!!!!!!!!! A(I,J)=1 ELSE A(I,J)= A(I-1,J)+ A(I-1,J-1) END IF END DO END DO DO I=1,10 PRINT ‘(10I4,( A(I,J),J=1,I)’ END DO END PROGRAM EXAM