sistema_lineare.py

Script

# Risolutore di sistemi lineari con il metodo di Cramer

print("Programma per risolvere un sistema lineare in x ed y")
print("⎰ a x + b y = e")
print("⎱ c x + d y = f")
print("Inserisci la prima equazione")
a = float(input("a = "))
b = float(input("b = "))
e = float(input("e = "))
print("Inserisci la seconda equazione")
c = float(input("c = "))
d = float(input("d = "))
f = float(input("f = "))

# calcolo il determinante D
D = a*d-b*c
# calcolo il determinante Dx
Dx = d*e-b*f
# calcolo il determinante Dy
Dy = a*f-c*e
# calcolo le soluzioni
x = Dx/D
y = Dy/D

print("La soluzione del sistema è:")
print("⎰ x =", x)
print("⎱ y =", y)

Shell

Programma per risolvere un sistema lineare in x ed y
⎰ a x + b y = e
⎱ c x + d y = f
Inserisci la prima equazione
a = 1
b = 2
e = 3
Inserisci la seconda equazione
c = 4
d = 5
f = 6
La soluzione del sistema è:
⎰ x = -1.0
⎱ y = 2.0

Script
(versione con il controllo delle situazioni di sistema indeterminato e impossibile)

# Risolutore di sistemi lineari con il metodo di Cramer

print("Programma per risolvere un sistema lineare in x ed y")
print("⎰ a x + b y = e")
print("⎱ c x + d y = f")
print("Inserisci la prima equazione")
a = float(input("a = "))
b = float(input("b = "))
e = float(input("e = "))
print("Inserisci la seconda equazione")
c = float(input("c = "))
d = float(input("d = "))
f = float(input("f = "))

# calcolo il determinante D
D = a*d-b*c

# calcolo il determinante Dx
Dx = d*e-b*f

# calcolo il determinante Dy
Dy = a*f-c*e

if D!=0:   # sistema determinato
    # calcolo le soluzioni
    x = Dx/D
    y = Dy/D

    print("La soluzione del sistema è:")
    print("⎰ x =", x)
    print("⎱ y =", y)
elif Dx==Dy==0:
    print("Il sistema è indeterminato")
else:
    print("Il sistema è impossibile")

Shell

Programma per risolvere un sistema lineare in x ed y
⎰ a x + b y = e
⎱ c x + d y = f
Inserisci la prima equazione
a = 1
b = 2
e = 3
Inserisci la seconda equazione
c = 2
d = 4
f = 5
Il sistema è impossibile