新手求教The evaluation period for this product ends

zhongczlb

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	real a(70,70),y(70)

02	dimension s(70),nx(70),ny(70)

03	CHARACTER*15 FILE1,file2

04	write(*,*)'(yang pin)N,(bian liang)M=?'

05	read(*,*)n,m

06	print*,'file1.old=?'

07	read(*,'(a)') file1

08	open(3,file=file1,status='old')

09	c       read(3,*)((a(i,j),j=1,m),i=1,n)

10	read(3,*)((a(i,j),i=1,n),j=1,m)

11	read(3,*)(s(i),i=1,m)

12

close(3)

13	print*,'file2.new=?'

14	read(*,'(a)') file2

15	open(5,file=file2,status='new')

16	write(5,*)' YUAN SHI SHU JU '

17	c       CALL PIR(A,N,M,70,70)

18	call lxd(a,s,y,n,m)

19	do 20 i=1,m

20	20      nx(i)=i

21	do 25 i=1,n

22	25      ny(i)=i

23	do 30 i=1,m-1

24	do 40 j=i+1,m

25	if(s(j).gt.s(i)) then

26

t=s(i)

27

s(i)=s(j)

28

s(j)=t

29

nt=nx(i)

30	nx(i)=nx(j)

31

nx(j)=nt

32

end if

33	40      continue

34	30      continue

35	do 50 i=1,n-1

36	do 60 j=i+1,n

37	if(y(j).gt.y(i)) then

38

t=y(i)

39

y(i)=y(j)

40

y(j)=t

41

nt=ny(i)

42	ny(i)=ny(j)

43

ny(j)=nt

44

end if

45	60      continue

46	50      continue

47	s(m+1)=0.0

48	y(n+1)=0.0

49	write(5,*)'* * quan xi shu * *'

50	write(5,*)'        s(i)    s(i)-s(i+1)'

51	do 70 i=1,m

52	70      write(5,100) i,nx(i),s(i),s(i)-s(i+1)

53	write(5,*)'* * dei xiang quan * *'

54	write(5,*)'        y(i)    y(i)-y(i+1)'

55	do 80 i=1,n

56	80      write(5,100) i,ny(i),y(i),y(i)-y(i+1)

57	100     format('(',i3,')','(',i3,')',2f8.4)

58	write(5,200) n,m

59	200     format('n=',i4,'m=',i4)

60

end

61

c

62	subroutine lxd(x,d,y,n,m)

63	real x(70,70),y(70),d(70)

64	do 1 i=1,n

65

y(i)=0.0

66	do 1 j=1,m

67	y(i)=y(i)+x(i,j)*d(j)

68	1    continue

69

return

70

end

71

c

72	subroutine pir(r,n,m,N1,M1)

73	real r(N1,M1)

74	do 100 ib=1,m,10

75

ie=ib+9

76	if(ie-m) 2,2,1

77

1    ie=m

78	2    write(5,4) (i,i=ib,ie)

79	do 99 j=1,n

80	write(5,5) j,(r(j,k),k=ib,ie)

81	99    continue

82	100    continue

83	4      format(5x,10i8)

84	5      format(i5,10f8.4)

85

RETURN

86

end

这样一个fortran文件应该如何编译运行呢？尝试了一下codeblocks和intel visual fortran+vs2012，都得不到运行结果。求大神教。

zhongczlb

已被 @竹叶知风飘 大神解决！thank you！