fortran里的stop问题

Mr.Yes

麻烦各位大神能不能帮我看看最后一行的stop应该怎么改才可以让主程序的print*，b1，b2，b4输出？ps.如果可以，再麻烦可不可以教教我为什么主程序的N=14，子程序qn里的m=13才可以求逆？而如果两个都是13就求不了。。。我的数据就是13个数啊？谢谢各位了！附程序： 新建文本文档.txt (2.96 KB, 下载次数: 1)

2013-12-15 21:55 上传

点击文件名下载附件

1. program ex5external qninteger,parameter::n=14real x1(1,n),x2(1,n),x4(1,n),y(1,n),xx1(n,1),xx2(n,1),xx4(n,1),x11(n,n),x22(n,n),x44(n,n),b(n,n),a(n,n),b11(n,n),b22(n,n),b44(n,n),xx11(n,1),xx22(n,1),xx44(n,1),b1,b2,b4!读数据data x1/7,1,11,11,7,11,3,1,2,21,1,11,10/data x2/26,29,56,31,52,55,71,31,54,47,40,66,68/data x4/60,52,20,47,33,22,6,44,22,26,34,12,12/data y/78.5,74.3,104.3,87.6,95.9,109.2,102.7,72.5,93.1,115.9,83.8,113.3,109.4/!求转置矩阵do i=1,nxx1(i,1)=x1(1,i)xx2(i,1)=x2(1,i)xx4(i,1)=x4(1,i)end do !X与X'相乘do i=1,ndo j=1,nx11(i,j)=xx1(i,1)*x1(1,j)x22(i,j)=xx2(i,1)*x2(1,j)x44(i,j)=xx4(i,1)*x4(1,j)end doend do!求X11,X22,X44的逆矩阵,b矩阵是相应的逆矩阵call qn(x11)b11=x11call qn(x22)b22=x22call qn(x33)b44=x44!逆矩阵与X'相乘do i=1,ndo j=1,nxx11(i,1)=xx11(i,1)+b11(i,j)*xx1(i,1)xx22(i,1)=xx22(i,1)+b22(i,j)*xx2(i,1)xx44(i,1)=xx44(i,1)+b44(i,j)*xx4(i,1)end doend do!XX11,XX22,XX44与y相乘，求bdo i=1,nb1=xx11(i,1)*y(1,i)b2=xx22(i,1)*y(1,i)b4=xx44(i,1)*y(1,i)end doprint*,b1,b2,b4endmodule LinearAlgebra  implicit nonecontains! 求逆矩阵subroutine inverse(A,IA)  implicit none  real    :: A(:,:), IA(:,:)  real, allocatable :: B(:,:)  integer :: i,j,N  N = size(A,1)    allocate(B(N,N))  ! 先把IA设定成单位矩阵  forall(i=1:N,j=1:N,i==j) IA(i,j)=1.0  forall(i=1:N,j=1:N,i/=j) IA(i,j)=0.0  ! 保存原先的矩阵A, 使用B来计算  B=A   ! 把B化成对角线矩阵(除了对角线外,都为0)  call Upper(B,IA,N) ! 先把B化成上三角矩阵  call Lower(B,IA,N) ! 再把B化成下三角矩阵  ! 求解  forall(i=1:N) IA(i,:)=IA(i,:)/B(i,i)   returnend subroutine! 输出矩阵的子程序subroutine output(matrix)  implicit none  real    :: matrix(:,:)  integer :: m,n,i  character(len=20) :: for='(??(1x,f7.3))'  m = size(matrix,1)  n = size(matrix,2)  ! 用字符串来设定输出格式  write( FOR(2:3), '(I2)' ) N  do i=1,N        write( *, FMT=FOR ) matrix(i,:)  end do  returnend subroutine output! 求上三角矩阵的子程序subroutine Upper(M,S,N)  implicit none  integer :: N  real    :: M(N,N)  real    :: S(N,N)  integer :: I,J  real :: E  do I=1,N-1    do J=I+1,N                    E=M(J,I)/M(I,I)      M(J,I:N)=M(J,I:N)-M(I,I:N)*E      S(J,:)=S(J,:)-S(I,:)*E    end do  end do  returnend subroutine Upper! 求下三角矩阵的子程序subroutine Lower(M,S,N)  implicit none  integer :: N  real    :: M(N,N)  real    :: S(N,N)  integer :: I,J  real :: E  do I=N,2,-1    do J=I-1,1,-1                 E=M(J,I)/M(I,I)      M(J,1:N)=M(J,1:N)-M(I,1:N)*E      S(J,:)=S(J,:)-S(I,:)*E    end do  end do  returnend subroutine Lowerend module! 求解联立式subroutine qn(a)  use LinearAlgebra  implicit none  integer, parameter :: m=13 ! Size of Matrix  real :: A(m,m)  real :: IA(m,m)  integer :: i  write(*,*) "原矩阵"  call output(A)  call inverse(A,IA)  write(*,*) "逆矩阵"  call output(IA)stopend subroutine qn

program ex5external qninteger,parameter::n=14real x1(1,n),x2(1,n),x4(1,n),y(1,n),xx1(n,1),xx2(n,1),xx4(n,1),x11(n,n),x22(n,n),x44(n,n),b(n,n),a(n,n),b11(n,n),b22(n,n),b44(n,n),xx11(n,1),xx22(n,1),xx44(n,1),b1,b2,b4!读数据data x1/7,1,11,11,7,11,3,1,2,21,1,11,10/data x2/26,29,56,31,52,55,71,31,54,47,40,66,68/data x4/60,52,20,47,33,22,6,44,22,26,34,12,12/data y/78.5,74.3,104.3,87.6,95.9,109.2,102.7,72.5,93.1,115.9,83.8,113.3,109.4/!求转置矩阵do i=1,nxx1(i,1)=x1(1,i)xx2(i,1)=x2(1,i)xx4(i,1)=x4(1,i)end do !X与X'相乘do i=1,ndo j=1,nx11(i,j)=xx1(i,1)*x1(1,j)x22(i,j)=xx2(i,1)*x2(1,j)x44(i,j)=xx4(i,1)*x4(1,j)end doend do!求X11,X22,X44的逆矩阵,b矩阵是相应的逆矩阵call qn(x11)b11=x11call qn(x22)b22=x22call qn(x33)b44=x44!逆矩阵与X'相乘do i=1,ndo j=1,nxx11(i,1)=xx11(i,1)+b11(i,j)*xx1(i,1)xx22(i,1)=xx22(i,1)+b22(i,j)*xx2(i,1)xx44(i,1)=xx44(i,1)+b44(i,j)*xx4(i,1)end doend do!XX11,XX22,XX44与y相乘，求bdo i=1,nb1=xx11(i,1)*y(1,i)b2=xx22(i,1)*y(1,i)b4=xx44(i,1)*y(1,i)end doprint*,b1,b2,b4endmodule LinearAlgebra  implicit nonecontains! 求逆矩阵subroutine inverse(A,IA)  implicit none  real    :: A(:,:), IA(:,:)  real, allocatable :: B(:,:)  integer :: i,j,N  N = size(A,1)    allocate(B(N,N))  ! 先把IA设定成单位矩阵  forall(i=1:N,j=1:N,i==j) IA(i,j)=1.0  forall(i=1:N,j=1:N,i/=j) IA(i,j)=0.0  ! 保存原先的矩阵A, 使用B来计算  B=A   ! 把B化成对角线矩阵(除了对角线外,都为0)  call Upper(B,IA,N) ! 先把B化成上三角矩阵  call Lower(B,IA,N) ! 再把B化成下三角矩阵  ! 求解  forall(i=1:N) IA(i,:)=IA(i,:)/B(i,i)   returnend subroutine! 输出矩阵的子程序subroutine output(matrix)  implicit none  real    :: matrix(:,:)  integer :: m,n,i  character(len=20) :: for='(??(1x,f7.3))'  m = size(matrix,1)  n = size(matrix,2)  ! 用字符串来设定输出格式  write( FOR(2:3), '(I2)' ) N  do i=1,N        write( *, FMT=FOR ) matrix(i,:)  end do  returnend subroutine output! 求上三角矩阵的子程序subroutine Upper(M,S,N)  implicit none  integer :: N  real    :: M(N,N)  real    :: S(N,N)  integer :: I,J  real :: E  do I=1,N-1    do J=I+1,N                    E=M(J,I)/M(I,I)      M(J,I:N)=M(J,I:N)-M(I,I:N)*E      S(J,:)=S(J,:)-S(I,:)*E    end do  end do  returnend subroutine Upper! 求下三角矩阵的子程序subroutine Lower(M,S,N)  implicit none  integer :: N  real    :: M(N,N)  real    :: S(N,N)  integer :: I,J  real :: E  do I=N,2,-1    do J=I-1,1,-1                 E=M(J,I)/M(I,I)      M(J,1:N)=M(J,1:N)-M(I,1:N)*E      S(J,:)=S(J,:)-S(I,:)*E    end do  end do  returnend subroutine Lowerend module! 求解联立式subroutine qn(a)  use LinearAlgebra  implicit none  integer, parameter :: m=13 ! Size of Matrix  real :: A(m,m)  real :: IA(m,m)  integer :: i  write(*,*) "原矩阵"  call output(A)  call inverse(A,IA)  write(*,*) "逆矩阵"  call output(IA)stopend subroutine qn

topmad

你这代码贴的有问题啊...编辑一下   没法看啊

言深深

代码不仅仅贴乱了，还不完全，直接粘贴到里面吧，不要使用贴代码的工具贴···

Leslie

{:soso_e163:}{:soso_e163:}

lysx

完全没法看！