포트란으로 생성한 행렬을 Octave로 간단히 시각화하기

*본 글은 ChatGPT를 활용해 작성하였습니다.

​

(1) 다음과 같은 포트란 코드를 작성하여, 3X3 행렬을 array.txt 파일로 저장한다.

program test real :: a(3,3) integer :: i, j, n = 3 a = reshape((/ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0 /), shape(a)) open(unit=10, file='array.txt') do i = 1, n write(10,*) (a(i,j), j=1,n) end do close(10) end

(2) array.txt의 출력값은 다음과 같다.

1.000000 4.000000 7.000000 2.000000 5.000000 8.000000 3.000000 6.000000 9.000000

(3) octave 실행시켜 array.txt 파일이 있는 폴더로 이동한 다음 아래와 같이 입력한다. (터미널에서 octave를 실행시켜 입력해도 된다.)

z=load('array.txt'); mesh(z);

(4) octave로 플롯한 그래프는 다음과 같다.

(5) 빨간 박스의 아이콘을 클릭하면 마우스를 움직여서 원하는 방향에서 그래프를 볼 수 있다.

(6) 'mesh(z)' 명령어보다 색깔입혀진 그래프를 보려면 'surf(z)' 명령어를 사용하면 된다.

​