PROGRAM GaussElimi
PARAMETER(N=3)
DIMENSION A(N,N),B(N),X(N)
 
DATA A/1,2,3,1,3,2,-1,5,3/
 
DATA B/2,-3,6/
PRINT 20,((A(I,J),J=1,3),B(I),I=1,3)
20   FORMAT(3X,3F8.3,3X,F8.3)
DO K=1,N-1
P=A(K,K)
DO I=K+1,N
F=A(I,K)/P
DO J=K+1,N
A(I,J)=A(I,J)-F*A(K,J)
ENDDO
B(I)=B(I)-F*A(K,J)
ENDDO
ENDDO
PRINT 20,((A(I,J),J=1,3),B(I),I=1,3)
X(N)=B(N)/A(N,N)
DO K=N-1,1,-1
SUM=0.0
DO J=K+1,N
SUM=SUM+A(K,J)*X(J)
ENDDO
X(K)=(B(K)-SUM)/A(K,K)
ENDDO
WRITE(6,*)"SOLUTIONS"
PRINT*,"X=",X(1)
PRINT*,"Y=",X(2)
PRINT*,"Z=",X(3)
END PROGRAM

UFJPR1JBTSBHYXVzc0VsaW1pClBBUkFNRVRFUihOPTMpCkRJTUVOU0lPTiBBKE4sTiksQihOKSxYKE4pCgpEQVRBIEEvMSwyLDMsMSwzLDIsLTEsNSwzLwoKREFUQSBCLzIsLTMsNi8KUFJJTlQgMjAsKChBKEksSiksSj0xLDMpLEIoSSksST0xLDMpCjIwICAgRk9STUFUKDNYLDNGOC4zLDNYLEY4LjMpCkRPIEs9MSxOLTEKUD1BKEssSykKRE8gST1LKzEsTgpGPUEoSSxLKS9QCkRPIEo9SysxLE4KQShJLEopPUEoSSxKKS1GKkEoSyxKKQpFTkRETwpCKEkpPUIoSSktRipBKEssSikKRU5ERE8KRU5ERE8KUFJJTlQgMjAsKChBKEksSiksSj0xLDMpLEIoSSksST0xLDMpClgoTik9QihOKS9BKE4sTikKRE8gSz1OLTEsMSwtMQpTVU09MC4wCkRPIEo9SysxLE4KU1VNPVNVTStBKEssSikqWChKKQpFTkRETwpYKEspPShCKEspLVNVTSkvQShLLEspCkVORERPCldSSVRFKDYsKikiU09MVVRJT05TIgpQUklOVCosIlg9IixYKDEpClBSSU5UKiwiWT0iLFgoMikKUFJJTlQqLCJaPSIsWCgzKQpFTkQgUFJPR1JBTQ==