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
        !+++++++++++++++++++++++++++++++++++++++++
        ！高斯消去法的核心部分
        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）你在程序中使用中文的感叹号去注释代码是不行的
2）主程序中的隐式声明不全，所以代码是有问题的

消除了 ...

4楼: Originally posted by physliebe at 2014-02-09 12:28:40
谢谢，主要是我用了codeblocks13.12作为集成开发环境，没有配置好Toolchain excutables。标准写法是不是real,dimension(8)呢？...