..................................................

نقل قول:::

کد فرترن برنامه ی روش تکرار گاوس سایدل...کدی که نوشته شده است ابتدا n که تعداد معادلاتی هست که میخواید وارد کنید رو گرفته و ماتریس n در n+1 رو تشکییل میدهد.سپس مقادیر ماتریس را سطر به سطر میگیرد.سپس نرخ همگرایی رو میگیرد و با توجه به نرخ همگرایی شروع به محاسبه به روش گاوس سایدل میکند.خطاها هم خطای نسبی هست و نرخ همگرایی نسبت به خطای نسبی مرتبه بینهایت سنجیده میشود.در آخر هم جواب ها را به همراه تعداد دفعات تکرار نمایش میدهد.توضیحاتی که در برنامه به کاربر میدهد به فینگیلیش نوشته شده است.اگه انگلیسی میخواید خودتون عوضش کنید.…

program gauss_sidel
implicit none
INTEGER::i,j,i1,j1,k,t1,n,j2,k1
REAL::s,e
REAL,ALLOCATABLE::a(:,:),x(:),y(:),t(:)
PRINT*
PRINT*
PRINT*,"             1            2       . . .      n       n+1"
PRINT*,"       _     -            -         -        -        -       _ "
PRINT*,"    1 |  X(1)a(1,1) + X(2)a(1,2) + ... + X(n)a(1,n) = a(1,n+1) |"
PRINT*,"    2 |  X(1)a(2,1) + X(2)a(2,2) + ... + X(n)a(2,n) = a(2,n+1) |"
PRINT*,"    . |    .   .        .  .        .      .  .     .   . . .  |"
PRINT*,"    . |    .   .        .  .        .      .  .     .   . . .  |"
PRINT*,"    . |    .   .        .  .        .      .  .     .   . . .  |"
PRINT*,"    . |    .   .        .  .        .      .  .     .   . . .  |"
PRINT*,"    n |_ X(1)a(n,1) + X(2)a(n,2) + ... + X(n)a(n,n) = a(n,n+1)_|"
PRINT*,"                                                              n*n+1"
PRINT*
PRINT*,"baraye matris n*n+1 bala lotfan n (tedad moadalat) ra vared konid :"
READ*,n
ALLOCATE (a(n,n+1),x(n),y(n),t(n))
do i1=1,n
PRINT*,"khate",i1,"ra vared konid ( az a(",i1,", 1)  ta  a(",i1,",",n+1,")) : "
READ*,a(i1,:)
end do
PRINT*,"nerkh hamgarayi ra vared konid : "
READ*,e
do j2=1,n
x(j2)=0
y(j2)=0
end do
k=1
do
do i=1,n
s=0
do j=1,n
if (j.ne.i) then
s=s+x(j)*a(i,j)
end if
end do
x(i)=(a(i,n+1)-s)/a(i,i)
end do
do t1=1,n
t(t1)=(abs(x(t1)-y(t1)))/ABS(x(t1))
end do
if (MAXVAL(t) < e) then
exit
else
k=k+1
y(1:n)=x(1:n)
end if
end do
PRINT*,"javab ha dar tekrar ",k," ba hadse avaliye X(1:n)=0  ::"
PRINT*
do k1=1,n
PRINT*,"  X(",k1,") = ",x(k1)
end do
PRINT*
end