.....................................................

نقل قول::

کد فرترن محاسبه دترمینان ماتریس n*n.....کدی که نوشته شده است اول از شما n که همان تعداد سطر و ستون ماتریس مربعی n*n است را میگیرد.بعد خط به خط درایه های ماتریس را از شما گرفته و برنامه در هر مرحله نگاه میکند میبیند که روی قطر اصلی عدد 0 وجود دارد یا نه.اگر باشد ان سطر را با سطری که درایه های آن 0 ندارد عوض میکند.این کار را برای همه سطرها در هر مرحله انجام میدهد تا همه آنها درست باشند.حالا شروع میکند به تبدیل این ماتریس به یک ماتریس بالا مثلثی.میدانید که یکی از روش های محاسبه دترمینان همین است. بعد از اینکه بالا مثلثی شد دترمینان برابر حاصلضرب درایه های قطر اصلی ماتریس می شود. این حاصلضرب را محاسبه میکند و دترمینان را نشان میدهد...…

program determinan
implicit none
INTEGER::i,j,g,f,r,j1,k,s,n
REAL::landa,t
REAL,ALLOCATABLE::a(:,:),b(:)
PRINT*
PRINT*," for matrix(n,n) input n : "
READ*,n
ALLOCATE (a(n,n),b(n))
do g=1,n
PRINT*," input line",g,":"
READ*,a(g,:)
end do
k=0
do i=1,n
if (a(i,i)==0) then
do j=i+1,n
if (a(j,i)/=0) then
k=k+1
b=a(j,:)
a(j,:)=a(i,:)
a(i,:)=b
end if
end do
end if
do j1=i+1,n
landa=-a(j1,i)/a(i,i)
a(j1,:)=(landa*a(i,:))+a(j1,:)
end do
end do
PRINT*
PRINT*
t=1
do s=1,n
t=t*a(s,s)
end do
if (MOD(k,2)==0) then
PRINT*,"   Determinant = ",t
else
PRINT*,"   Determinant = ",-t
end if
PRINT*
PRINT*,"---------"
end