Schrijf de functie solve() met als argumenten:

• een NumPy-tabel adie een vierkante matrix voorstelt met $$M$$ rijen en $$M$$ kolommen
• een NumPy-tabel b, die uit $$M$$ rijen en 1 kolom bestaat.

Het resultaat van de functie is een NumPy-tabel met $$M$$ rijen en 1 kolom (een kolomvector) die de oplossing van het bijhorende stelsel voorstelt (met a de coëfficiëntenmatrix en b het rechterlid). Ga hierbij als volgt tewerk:

• Construeer de verhoogde matrix (voeg aan a een extra kolom toe, die bestaat uit de elementen van b)
• Breng deze matrix naar een rij-echelongedaante via de functie echelon()
• Bepaal de oplossing van dit stelsel door achterwaartse substitutie.

Je hoeft de code voor de functie echelon() niet mee in te dienen op het Dodona-platform.

Voorbeeld

a = np.array([[1., 2., 3.], [7., 10., 12.], [9., 21., 3.]])
b = np.array([[10.], [8.], [6.]])
x = solve(a, b)
#[[-14.634146341463415], [5.951219512195123], [4.2439024390243905]]