当前位置:首页
建筑知识
阅读以下FORTRAN程序:READ(*.*)MIF(M.LT. 100.AND.M.GT999) STOPK=M/100L=(M-K*100)10N=M-K*100-L*10MM=K+NWRITE(
阅读以下FORTRAN程序:READ(*.*)MIF(M.LT. 100.AND.M.GT999) STOPK=M/100L=(M-K*100)10N=M-K*100-L*10MM=K+NWRITE(*,(2X,13)) MMEND键盘输入600,则程序的输出结果是:
A 、8
B 、3
C 、7
D 、6
【正确答案:D】
K=6, L= (0-10)0 =0, N= (60100-0*1)1) =0.MN=K+N-6+0-6。
PROGRAM MAIN !主程序
DOUBLE PRECISION A,AA,XX,XX1,XX2,TYY0,TYX0,YY0,YX0!定义双精度变量
DIMENSION A(101),AA(101)!定义数组
OPEN(10,FILE='OUTPUT.TXT',STATUS='UNKNOWN')!建立输出文件通道
WRITE(*,*)'此程序用于求解一元N次方程最接近给定的X0的一个解。'!屏幕提示
WRITE(*,*)'请输入N(≤100,0=EXIT):'!屏幕提示
2 READ(*,*)N !从屏幕读取N
IF(N.EQ.0) THEN !判定如果N=0则程序转去标识符号10000处
GOTO 10000
END IF
IF((N.Lt.0).OR.(N.GT.100)) THEN !判定如果N小于0或N大于100 则屏幕提示重新输入N
WRITE(*,*)'请重新输入N(≤100,0=EXIT):'
GOTO 2 !转到 标识2句处执行
ENDIF
WRITE(*,*)'请依次输入各项系数AN~A0(高次幂~低次幂)。'屏幕提示
WRITE(*,*)'每个数据后回车:'屏幕提示
DO 5 I=1,N+1 输入数组A(I)
READ(*,*) A(I)
5 CONTINUE
IF(MOD(N,2).EQ.0) THEN
WRITE(*,*)' '
WRITE(*,*)'对N为偶数的情况,请仔细校核。'
WRITE(*,*)' '
END IF
WRITE(*,*)'请输入X0:'
READ(*,*) X0
DO 10 I=1,N 输入数组AA(I)
AA(I)=(N-I+1)*A(I)
10 CONTINUE
ICOUNT=0 给ICOUNT赋值
XX=X0给XX赋值
DO 30 I=1,1000000 对I循环从1到10000000
ICOUNT=ICOUNT+1每循环一次计数加一
YY0=0.给YY0赋值
YX0=0.给YX0赋值
DO 40 J=1,N
TYY0=A(J)*(XX**(N-J+1)) 计算TYY0,YY0
YY0=YY0+TYY0
40 CONTINUE
YY0=YY0+A(N+1) 计算YY0
DO 50 K=1,N-1套另一个循环 对K循环从1到N-1
TYX0=AA(K)*(XX**(N-K))
YX0=YX0+TYX0
50 CONTINUE
YX0=YX0+AA(N) 计算YX0
IF(YX0.EQ.0) THEN 判别 YX0如果=0 则
XX=XX+1XX加1
ICOUNT=ICOUNT-1 计数减1
GOTO 30 结束I这一步大循环开始I的下一步循环
END IF
IF(MOD(I,2).EQ.1)THEN判别
XX1=XX-YY0/YX0
XX=XX1
ELSE
XX2=XX-YY0/YX0
XX=XX2
END IF
XX=.5*((XX1+XX2)-(-1)**I*(XX1-XX2)) 计算赋值
IF(ICOUNT.GE.1000) THEN 判别 如果计数大于等于1000则结束程序
WRITE(*,*)' ' 屏幕输出提示
WRITE(*,*)'此方程无实数解。'屏幕输出提示
GOTO 10000 则结束程序
END IF
IF((ABS(XX2-XX1).LE.1.D-15).AND.(ICOUNT.GT.1)) THEN
IF((ABS(XX).GT.1.D-6).AND.(ABS(XX).LT.9.D6)) THEN
GOTO 1000
ELSE
GOTO 1001
END IF
END IF
30 CONTINUE
1000WRITE(*,100)XX
GOTO 1010
1001WRITE(*,101)XX
100 FORMAT(8X,'X=',F25.6)
101 FORMAT(8X,'X=',D25.6)
1010WRITE(*,*)' '
WRITE(*,*)'本程序计算结果仅供参考。'
WRITE(*,*)' '
PAUSE
10000 STOP
END
是write语句的格式问题
character m
write(m,'(i2)') i
此处不是很理解,m是一个字符(character),为何用作文件的编号呢?
用问题再交流