Buena biblioteca, si deseas puedes agregar lo siguiente
Code:
Anexo un seudocódigo para solucionar por formula un sistema de ecuaciones lineales en dos variables
el siguiente código se ejecuta con PSEINT
http://pseint.sourceforge.net/
Proceso sistema_de_ecuaciones_lineales
//(c) By JaiMeza aka compSystems mayo 12 2015 www.jaimeza.org
//(c) para ejecutarlo use http://pseint.sourceforge.net/
imprimir "sistema de ecuaciones lineales"
imprimir "en dos variables"
imprimir ""
imprimir "a1*x + b1*y + c1 = 0"
imprimir "a2*x + b2*y + c2 = 0"
imprimir ""
imprimir "/!\ a1*b1 <> 0, a2*b2 <> 0,"
definir a1, b1, c1, a2, b2, c2, d, x, y_ como reales
definir multiplo_ como logico
multiplo_ = falso
imprimir ""
imprimir "entre a1:"
repetir
leer a1
hasta que a1 <> 0
imprimir "entre b1:"
repetir
leer b1
hasta que b1 <> 0
imprimir "entre c1:"
leer c1
imprimir ""
imprimir "entre a2:"
repetir
leer a2
hasta que a2 <> 0
imprimir "entre b2:"
repetir
leer b2
hasta que b2 <> 0
imprimir "entre c2:"
leer c2
// imprimir a1
// imprimir b1
// imprimir a2
// imprimir b2
d = (a1*b2) - (a2*b1)
//imprimir d
imprimir ""
Si d=0 Entonces
multiplo_ = ((a1 es multiplo de a2) o (a2 es multiplo de a1)) y ((b1 es multiplo de b2) o (b2 es multiplo de b1) ) y ((c1 es multiplo de c2) o (c2 es multiplo de c1))
Si multiplo_ = verdadero Entonces
imprimir "Sistema Compatible Indeterminado"
imprimir "existen infinitas soluciones o intersecciones"
imprimir "las dos rectas son paralelas y estan juntas"
Sino
imprimir "Sistema Incompatible"
imprimir "no hay solucion o puntos comunes"
imprimir "las dos rectas son paralelas y separadas"
Fin Si
Sino
x = (b2*c1 - b1*c2)/d
y_ = (a1*c2 - a2*c1)/d
imprimir "Sistema Compatible Determinado "
imprimir "existe solucion unica, las dos rectas se cortan en un solo punto"
imprimir ""
imprimir "x=",x
imprimir "y=",y_
imprimir ""
imprimir "la coordenada de interseccion es (", x, ", ", y_, ")"
Fin Si
FinProceso
// ejemplos
// 1:
// 9*x-1 + -6*y = 3
// 4*x-1 + 6*y = 36
// a1 = 9, b1=-6, c1=3
// a2 = 4, b2=6, c2=36
// x = 3, y=4
// 2:
// 2*x-1 = 4*y
// 5+y = 10*x
// a1 = 2, b1=-4, c1=1
// a2 = 10, b2=-1, c2=5
// x = 0.5, y=0
// 3:
// 2.5*x + 4.6*y = 5
// 1*x - 2.9*y = 7
// a1 = 2.5, b1=4.6, c1=5
// a2 = 1, b2=-2.9, c2=7
// x=3.94092827 y=-1.0548523207
// 4:
// 1.2*x + 2.2*y = 1
// 6*x+11*y = 5
// a1 = 1.2, b1=2.2, c1=1
// a2 = 6, b2=11, c2=5
// no hay solución