행렬계산에 관해서

[흰둥이-ㅡv]

(dn1/dt) = (5 6)(n1)

(dn2/dt) = (7 8)(n2) 참고로 행렬을 쓴겁니다.

이것을 포트란으로 짜볼려고 하는데요, 저의 의도는 n1의 값이 시간에 따라서 변하고 그 값이 다시 저 행렬식에 들어가서

계산이 되었으면 합니다.

이번에 포트란을 배운적이 없어서 너무 힘이듭니다. 조언좀 부탁드릴게요...모가 문제인지도 모르겠습니다.

제가 한거는요

program phone2 
 
 implicit none

 integer :: A(2,2)
 integer :: B(2), C(2)
 integer n1, n2, dt, i, dn1, dn2

  
 write(*,*) "n1, n2의 값을 입력하시오."
 write(*,*) 'n1='
 read(*,*) n1
 write(*,*) 'n2='
 read(*,*) n2
 write(*,*) 'dt='
 read(*,*) dt
 
 A=reshape((/5,6,8,4/),shape(A))
 B=reshape((/dn1,dn2/),shape(B))

 write(*,*) 'A=    '
 write(*,"(2i2)")  A !writes 5 6
                     !       8 4
 write(*,*) 'B='
 write(*,"(2i1)")  B

 do i=1,1000
  
  dn1= (5*n1+6*n2)*dt
  dn2= (8*n1+4*n2)*dt
  n1(i)= n1 + dn1
  n2(i)= n2 + dn2
  C(i) = matmul(B,A)

 end do
 

 write(*,*) 'C =    A     *     B'
 write(*,"(2i1)")  C

 stop

 end program phone2