[問題] 反矩陣的副程式&主程式 (已解決)
副程式的部分是課本上的 應該不會有錯
現在就是在煩惱怎麼寫出主程式~ ~
需要求20*20矩陣的反矩陣
所以我想先用2*2來跑看看 可以執行但是一直觸發中斷點
應該是我主程式寫錯的關係
有人可以幫看一下主程式哪裡有問題嗎ˊ ˋ"
implicit none
real a(2,2),kk(2,2)
data a/1,3,2,4/
call inver(a,kk)
write(*,*)kk
pause
end
subroutine inver(n,ss)
real*8 ss,d,tt
integer n,i,j,k
dimension ss(n,n)
do i=1,n
d=1.d0/ss(i,i)
tt=-d
do j=1,n
ss(i,j)=ss(i,j)*tt
end do
do k=1,i-1
tt=ss(k,i)
do j=1,i-1
ss(k,j)=ss(k,j)+tt*ss(i,j)
end do
do j=i+1,n
ss(k,j)=ss(k,j)+tt*ss(i,j)
end do
ss(k,i)=tt*d
end do
do k=i+1,n
tt=ss(k,i)
do j=1,i-1
ss(k,j)=ss(k,j)+tt*ss(i,j)
end do
do j=i+1,n
ss(k,j)=ss(k,j)+tt*ss(i,j)
end do
ss(k,i)=tt*d
end do
ss(i,i)=d
end do
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 111.252.47.251
→
05/20 20:29, , 1F
05/20 20:29, 1F
→
05/20 22:15, , 2F
05/20 22:15, 2F
→
05/20 22:15, , 3F
05/20 22:15, 3F
→
05/20 22:22, , 4F
05/20 22:22, 4F
→
05/20 22:23, , 5F
05/20 22:23, 5F
→
05/20 22:24, , 6F
05/20 22:24, 6F
推
05/20 22:49, , 7F
05/20 22:49, 7F
→
05/20 22:50, , 8F
05/20 22:50, 8F
→
05/20 22:51, , 9F
05/20 22:51, 9F
→
05/20 22:51, , 10F
05/20 22:51, 10F
→
05/20 22:56, , 11F
05/20 22:56, 11F
→
05/20 22:57, , 12F
05/20 22:57, 12F
→
05/20 23:02, , 13F
05/20 23:02, 13F
→
05/20 23:02, , 14F
05/20 23:02, 14F
推
05/20 23:07, , 15F
05/20 23:07, 15F
把呼叫副程式改成 call inver(2,kk)
應該就會把2傳到 subroutine inver(n,ss) 的n吧..
→
05/20 23:09, , 16F
05/20 23:09, 16F
→
05/20 23:09, , 17F
05/20 23:09, 17F
我把主程式改成
implicit none
real kk(2,2)
data kk/1,3,2,4/
call inver(2,kk)
write(*,*)kk
pause
end
宣告丟過去的kk是矩陣 資料是 1 3 2 4
還是錯在主程式的第四排 call過去的資料..
※ 編輯: joelife 來自: 111.252.47.251 (05/20 23:20)
推
05/20 23:43, , 18F
05/20 23:43, 18F
→
05/20 23:44, , 19F
05/20 23:44, 19F
有.. 我是ftn95
但是執行之後依然跳出錯誤 http://ppt.cc/ahtR <- 執行畫面..
※ 編輯: joelife 來自: 111.252.47.251 (05/20 23:50)
→
05/21 01:19, , 20F
05/21 01:19, 20F
→
05/21 01:20, , 21F
05/21 01:20, 21F
推
05/21 10:33, , 22F
05/21 10:33, 22F
真的已經解決了 感謝大家的幫忙Q_Q
晚點挑戰20*20矩陣 還要CALL TIME ˋ(′_‵||)ˊ
※ 編輯: joelife 來自: 111.252.43.122 (05/21 11:38)
推
05/25 15:00, , 23F
05/25 15:00, 23F
→
05/25 18:48, , 24F
05/25 18:48, 24F
→
05/26 20:35, , 25F
05/26 20:35, 25F
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
