Schrijf de functie solve_multi() 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 $$K$$ kolommen componenten bestaat.

Het resultaat van de functie is een NumPy-tabel met $$M$$ rijen en $$K$$ kolommen die de oplossing van het bijhorende multi-stelsel voorstelt (met a de coëfficiëntenmatrix en b het rechterlid van de matrixvergelijking $$\mathbf{ax} = \mathbf{b}$$). Ga hierbij als volgt tewerk:

• Construeer de verhoogde matrix (voeg aan a een extra kolommen toe, namelijk de kolommen van b)
• Breng deze verhoogde matrix naar een rij-echelongedaante via de functie echelon()
• Bepaal de oplossing van dit multistelsel 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., 5.], [8., 4.], [6., 2.]])
x = solve_multi(a, b)

#[[-14.634146341463415, -7.268292682926829], [5.951219512195123, 2.9024390243902434],      [4.2439024390243905, 2.154471544715447]]