real a(100,100),b(100) open(13,file='data.d') read (13,*) isize print *, 'The size is ',isize,' x ',isize read (13,*) ((a(i,j), j=1,isize), i=1,isize) do 30 i=1,isize print 20, (a(i,j), j=1,isize) 30 continue 20 format(100(f6.1, 1x)) read (13,*) (b(i), i=1,isize) print *, 'B is =' print 20, (b(i), i=1,isize) close(13) end