ALpMontana
КЯaL´da TaNıMaM КuяaL´da
Arkadaşlar bu program n bilinmeyenli,n denklemli sistemi Gauss Eleminasyon yöntemi ile çözer.n. mertebeden A katsayılar matrisi ve sağ taraf vektörü girilir ve x çözüm vektörünü bulur.
real,allocatable::a
print*,"Denklem sayisini giriniz:";read*,n
allocate(a(n,n),b
print*,"Katsayilar matrisini giriniz:";read*,a
print*,"Katsayilar matrisi:"
do i=1,n
do j=1,n
print '(F7.3\)',a(i,j)
enddo
print*
enddo
print*,"Denklem sisteminin sag taraf vektorunu giriniz:";read*,b
print*,"Sag taraf vektoru:"
print '(F7.3)',b
do k=1,n-1
do i=k+1,n
do j=k+1,n
a(i,j)=a(i,j)-a(i,k)*a(k,j)/a(k,k)
enddo
b(i)=b(i)-a(i,k)*b(k)/a(k,k)
enddo
enddo
print*,"Elemanter satir islemleri sonrasında"
print*," A matrisi:"
do i=1,n
do j=1,n
if(j<i) a(i,j)=0
print '(F7.3\)',a(i,j)
enddo
print*
enddo
print*," B matrisi:"
do i=1,n
print*,b(i)
enddo
do i=n,1,(-1)
x(i)=(b(i)+a(i,i)*x(i)-(dot_product(a(i,: ),x)))/a(i,i)
enddo
print*,"Cozum vektoru:"
print '(F7.3/)',x
end
not 2 nokta parantezlerde bır bosluk war gulen yuz