关于三角形若干算法的作业

li913

算法没问题，可以改进。

[Fortran] 纯文本查看 复制代码

PROGRAM sistemas
IMPLICIT NONE
REAL a1,b1,c1,a2,b2,c2,a3,b3,c3,m12,m13,m23,xa,ya,xb,yb,xc,yc!,F1,F2,F
integer i

OPEN (11,FILE='dato5.txt')
OPEN (12, FILE='result.txt')
do
READ (11,*,iostat=i) a1,b1,c1,a2,b2,c2,a3,b3,c3
if(i/=0) exit
m12= a1*b2 - a2*b1   ! determinante de la matriz m12
IF (m12 == 0) STOP ' el sistema no tiene solucion unica'
xa = (-c1*b2 + c2*b1) / m12
ya = (-c2*a1 + c1*a2) / m12 

m13= a1*b3 - a3*b1   ! determinante de la matriz m13
IF (m13 == 0) STOP ' el sistema no tiene solucion unica'
xb = (-c1*b3 + c3*b1) / m13
yb = (-c3*a1 + c1*a3) / m13

m23= a2*b3 - a3*b2   ! determinante de la matriz m23
IF (m23 == 0) STOP ' el sistema no tiene solucion unica'
xc = (-c2*b3 + c3*b2) / m23                                  
yc = (-c3*a2 + c2*a3) / m23 

WRITE (12,9000) xa,ya,xb,yb,xc,yc
9000 FORMAT (6(3X,F12.4))
write(*,*) a1*xa+b1*ya+c1,a2*xa+b2*ya+c2
write(*,*) a1*xb+b1*yb+c1,a3*xb+b3*yb+c3
write(*,*) a3*xc+b3*yc+c3,a2*xc+b2*yc+c2
end do
pause
end program sistemas