program Gauss_Seidel
implicit none
real(8),dimension(:,:),allocatable::a
real(8),dimension(:),allocatable::b,x,x_k1,x_result
integer::i,j,n
print,"输入维度n"
read(,)n
allocate(a(n,n),b(n),x(n),x_k1(n),x_result(n))
print,"输入nn的数组a"
read(,)((a(i,j),j=1,n),i=1,n)
print,"n的数组b"
read(,),b
print,"初始数组x0"
read(,),x
call Gauss_Seidel_dd(n,a,b,x,x_result)
print,"x的数值为：",x_result
pause
end program Gauss_Seidel
!****************************************************
subroutine Gauss_Seidel_dd(n,a,b,x,x_result) ! 书上迭代法
implicit none
integer::i,j,n,k
real(8),dimension(n,n)::a
real(8),dimension(n)::b,x,x_k1,x_result
real(8)::sum1,sum2
do k=1,100 !迭代次数
do i=1,n
sum1=0
do j=1,i-1
sum1=sum1+a(i,j)x_k1(j)
end do
sum2=0
do j=i+1,n
sum2=sum2+a(i,j)x(j)
end do
x_k1(i)=(b(i)-sum1-sum2)/a(i,i)
enddo
x=x_k1 !完成一次迭代
enddo
x_result=x
end subroutine Gauss_Seidel_dd
Gauss_Seidel
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- program lagrangeimplicit nonereal(8),dimension(:,:),alloc...
- Last updated December 5, 2012. Try searching this page fo...
