module gauss
    contains
    subroutine solve(A,b,x,N)
        implicit real*8(a-z)
        integer::i,k,N
        real*8::A(N,N),b(N),x(N)
        real*8::Aup(N,N),bup(N)
        real*8::Ab(N,N+1)
        Ab(1:N,1:N)=A
        Ab(:,N+1)=b
        !+++++++++++++++++++++++++++++++++++++++++
        £¡¸ß˹ÏûÈ¥·¨µÄºËÐIJ¿·Ö
        do k=1,N-1
            do i=k+1,N
                temp=Ab(i,k)/Ab(k,k)
                Ab(i,:)=Ab(i,:)-temp*Ab(k,:)
            end do
        end do
        Aup(:,:)=Ab(1:N,1:N)
        bup(:)=Ab(:,N+1)
        !+++++++++++++++++++++++++++++++++++++++++
        £¡Çó½âÉÏÈý½Ç¾ØÕóÐ͵ÄÏÈÐÔ·½³Ì×é
        call uptri(Aup,bup,x,N)   
        end subroutine solve
     subroutine uptri(A,b,x,N)        
        implicit real*8(a-z)         
        integer::i,j,N
        real*8::A(N,N),b(N),x(N)      
        x(N)=b(N)/A(N,N)              
        !»Ø´ú²¿·Ö
        do i=N-1,1,-1               
            x(i)=b(i)
            do j=i+1,N               
                x(i)=x(i)-a(i,j)*x(j)
            end do
            x(i)=x(i)/A(i,i)
        end do
    end subroutine uptri
end module gauss
program main
    use gauss
    implicit real*(a-z)
    integer,parameter::N=4
    integer::i,j
    real*8::A(N,N),b(N),x(N)
    open(unit=11,file='fin.txt')
    open(unit=12,file='fout.txt')
    read(11,*) ((A(i,j),j=1,N),i=1,N),b
    call solve(A,b,x,N)   !xÔÚÕâÀïÃæµÄ×÷ÓÃ
    write(12,101) x
    101 format(T5,'¸ß˹ÏûÈ¥·¨¼ÆËã½á¹û',/,T4,4(f12.8))
end program main

3Â¥: Originally posted by snoopyzhao at 2014-02-09 05:22:20
ÖÁÉÙÓà gfortran ±àÒëÄãÌù³öÀ´µÄ´úÂë²»Ö¹ÄãÉÏÃæÌáµ½µÄ´íÎ󣬵«°´ÕÕ³ö´íÐÅϢȥÐÞ¸ÄÓÖ²»»áÓÐʲôÎÊÌâ

1£©ÄãÔÚ³ÌÐòÖÐʹÓÃÖÐÎĵĸÐ̾ºÅÈ¥×¢ÊÍ´úÂëÊDz»ÐеÄ
2£©Ö÷³ÌÐòÖеÄÒþʽÉùÃ÷²»È«£¬ËùÒÔ´úÂëÊÇÓÐÎÊÌâµÄ

Ïû³ýÁË ...

4Â¥: Originally posted by physliebe at 2014-02-09 12:28:40
лл£¬Ö÷ÒªÊÇÎÒÓÃÁËcodeblocks13.12×÷Ϊ¼¯³É¿ª·¢»·¾³£¬Ã»ÓÐÅäÖúÃToolchain excutables¡£±ê׼д·¨ÊDz»ÊÇreal,dimension(8)ÄØ£¿...