Sistema De Ecuaciones Lineales Y Determinantes

program calculadora sist ec det

real::w1,w2,w3,w4,x1,x2,x3,x4,y1,y2,y3,y4,z1,z2,z3,z4,r1,r2,r3,r4

real::D,Dw,Dx,Dy,Dz,w,x,y,z

integer::a,b,c,f

print*,'Solucion de determinantes y sistemas lineales'

100 print*,'¿Que desea hacer?'

print*,'(1) Determinantes de una matriz'

print*,'(2) Sistema de ecuaciones lineales'

read(*,*)a

if (a==1) then

print*,'Determinantes de una matriz'

print*,'(1) 2x2'

print*,'(2) 3x3'

print*,'(3) 4x4'

read(*,*)c

if(c==1) then

print*,'Determinantes de 2x2'

print*,'Ingrese x1,x2'

read(*,*)x1,x2

print*,'Ingrese y1,y2'

read(*,*)y1,y2

D=x1*y2-y1*x2

print*,'El resultado de la Determinante=',D

else if (c==2) then

print*,'Determinantes de 3x3'

print*,'Ingrese x1,x2,x3'

read(*,*)x1,x2,x3

print*,'Ingrese y1,y2,y3'

read(*,*)y1,y2,y3

print*,'Ingrese z1,z2,z3'

read(*,*)z1,z2,z3

D=x1*y2*z3+x2*y3*z1+y1*z2*x3-(x3*y2*z1+x1*z2*y3+x2*y1*z3)

print*,'El resultado de la Determinante=',D

else if(b==3) then

print*,'Determinante de 4x4'

print*,'Ingrese w1,w2,w3,w4'

read(*,*)w1,w2,w3,w4

print*,'Ingrese x1,x2,x3,x4'

read(*,*)x1,x2,x3,x4

print*,'Ingrese y1,y2,y3,y4'

read(*,*)y1,y2,y3,y4

print*,'Ingrese z1,z2,z3,z4'

read(*,*)z1,z2,z3,z4

D=w1*(x2*y3*z4+x3*y4*z2+y2*z3*x4-(x4*y3*z2+x3*y2*z4+x2*y4*z3))

& -w2*(x1*y3*z4+x3*y4*z1+x4*y1*z3-(x4*y3*z1+x3*y1*z4+x1*y4*z3))

& +w3*(x1*y2*z4+x2*y4*z1+x4*y1*z2-(x4*y2*z1+x2*y1*z4+x1*y4*z2))

& -w4*(x1*y2*z3+x2*y3*z1+x3*y1*z2-(x3*y2*z1+x2*y1*z3+x1*y3*z2))

print*,'El resultado de la Determinante=',D

end if

else if (a==2) then

print*,'Sistema de ecuaciones lineales'

print*,'(1) 2x2'

print*,'(2) 3x3'

print*,'(3) 4x4'

read(*,*)b

if(b==1) then

print*,'Sistema de ecuaciones de 2x2'

print*,'Ingrese x1,x2'

read(*,*)x1,x2

print*,'Ingrese y1,y2'

read(*,*)y1,y2

print*,'Ingrese r1,r2'

read(*,*)r1,r2

D=x1*y2-y1*x2

Dx=r1*y2-y1*r2

Dy=x1*r2-r1*x2

x=Dx/D

y=Dy/D

print*,'Resultados'

print*,'x=',x

print*,'y=',y

else if(b==2) then

print*,'Sistema de ecuaciones de 3x3'

print*,'Ingrese x1,x2,x3'

read(*,*)x1,x2,x3

...