fortran语言解n元线性方程组
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
PROGRAM MAIN
IMPLICIT NONE
INTEGER I,J,K,N
REAL:: SUM=0
REAL,ALLOCATABLE::A(:,:)
REAL,ALLOCATABLE::B(:)
REAL,ALLOCATABLE::X(:)
PRINT*,'请输入方程组未知数个数'
READ*,N
ALLOCATE(A(N,N+1))
ALLOCATE(B(N+1))
ALLOCATE(X(N))
PRINT*,'请按列输入方程组的系数和列向量b(即输入增广矩阵A)'
READ*,A
DO K=1,N-1
DO WHILE(A(K,K)==0)
B(:)=A(K,:)
DO I=K,N-1
A(I,:)=A(I+1,:)
END DO
A(N,:)=B(:)
END DO
DO I=K+1,N
A(I,:)=A(I,:)-A(I,K)/A(K,K)*A(K,:)
END DO
END DO
X(N)=A(N,N+1)/A(N,N)
DO I=N-1,1,-1
DO J=I+1,N
SUM=SUM+A(I,J)*X(J)
END DO
X(I)=(A(I,N+1)-SUM)/A(I,I)
SUM=0
END DO
PRINT*ຫໍສະໝຸດ BaiduX
END
用其他语言编写的话，算法是一样的，只是要改一下其中的循环结构
IMPLICIT NONE
INTEGER I,J,K,N
REAL:: SUM=0
REAL,ALLOCATABLE::A(:,:)
REAL,ALLOCATABLE::B(:)
REAL,ALLOCATABLE::X(:)
PRINT*,'请输入方程组未知数个数'
READ*,N
ALLOCATE(A(N,N+1))
ALLOCATE(B(N+1))
ALLOCATE(X(N))
PRINT*,'请按列输入方程组的系数和列向量b(即输入增广矩阵A)'
READ*,A
DO K=1,N-1
DO WHILE(A(K,K)==0)
B(:)=A(K,:)
DO I=K,N-1
A(I,:)=A(I+1,:)
END DO
A(N,:)=B(:)
END DO
DO I=K+1,N
A(I,:)=A(I,:)-A(I,K)/A(K,K)*A(K,:)
END DO
END DO
X(N)=A(N,N+1)/A(N,N)
DO I=N-1,1,-1
DO J=I+1,N
SUM=SUM+A(I,J)*X(J)
END DO
X(I)=(A(I,N+1)-SUM)/A(I,I)
SUM=0
END DO
PRINT*ຫໍສະໝຸດ BaiduX
END
用其他语言编写的话，算法是一样的，只是要改一下其中的循环结构