全国县城天气预报:求 N阶行列式的程序~~~~

才10分就想要个递归调用的程序啊。
加个0我还考虑下

方法很多啊，你需要什么方法的？（LU，QR,SVD）什么语言编写的？(C,ForTran?)

给你一个ForTran编写的Schmidt_QR方法：

cccc --------------- 主程序 ---------------- cccc
ccc --------------- 主程序 ---------------- ccc
program test
real temp
real s,t
allocatable::A(:,:),B(:,:),R(:,:)
call preenter(n,'matrix_A.dat')
allocate(A(n,n),B(n,n),R(n,n))
call enter(n,A,'matrix_A.dat')

do m=1,n
call prebeita(A,B,R,n,m)
call beita_f2(B,n,m,s)
call beita(B,n,m,s)
call R_array(A,B,R,n,m,s)
end do
call output(B,R,n)
end program
cccc --------------主程序结束---------------- cccc
ccc ---------------------------------------- ccc
cc -------------以下为子程序--------------- cc
ccc ************* cc *************** ccc
subroutine preenter(n,filename)
character(*) filename

open(10,file=filename,status='old')
read(10,*) n
close(10)
end subroutine

cccc ************ cc **************** cccc

subroutine enter(n,A,filename)
real A(n,n)
character(*) filename
open(10,file=filename,status='old')
read(10,*)
read(10,*) ((A(i,j),j=1,n),i=1,n)
close(10)
end subroutine

c *********************************************
subroutine prebeita(A,B,R,n,m)
real A(n,n),B(n,n),R(n,n)
do i=1,n
b(i,m)=0.0
do k=1,m-1
b(i,m)=b(i,m)+r(k,m)*b(i,k)
end do
b(i,m)=a(i,m)-b(i,m)
end do
end subroutine

c ***********************************************
subroutine beita_f2(B,n,m,s)
real B(n,n),s
s=0.0
do i=1,n
s=s+b(i,m)*b(i,m)
end do
s=sqrt(s)
end subroutine

c ************************************************
subroutine beita(B,n,m,s)
real B(n,n),s
do i=1,n
b(i,m)=b(i,m)/s
end do
end subroutine

c *************************************************
subroutine inner_multiply(A,B,n,m,t,p)
real A(n,n),B(n,n),t
integer p
t=0.0
do i=1,n
t=t+a(i,p)*b(i,m)
end do
end subroutine

c **********************************************
subroutine R_array(A,B,R,n,m,s)
real A(n,n),B(n,n),R(n,n),t
r(m,m)=s
do i=m+1,n
call inner_multiply(A,B,n,m,t,i)
r(m,i)=t
end do
end subroutine
c ***********************************************

subroutine output(B,R,n)
real B(n,n),R(n,n)
open(10,file='vector_X.dat',status='unknown')
write(10,*) n
do i=1,n
write(10,*) (B(i,j),j=1,n)
end do
write(10,*) '**********************************'
do i=1,n
write(10,*) (R(i,j),j=1,n)
end do
end subroutine