Hur lösa ut vektorer från linjära system med SymPy?

Hej!

Hur löser man ut vektorer från linjära system med SymPy? Jag har sökt men det finns ingen bra metod som fungerar. Jag vill mest bara skriva upp ett uttryck som sedan anropa detta uttryck med funktionen solve () och sedan lägga till ett argument i funktionen som pekar på vad som ska lösas ut. Men det fungerar inte.

Första problemet.

X = Ax + Bu

X är en 1x3 matris 

A är en 3x3 matris

x är en 1x3 matris som ska lösas ut

B är en 1x3 matris 

u är 1x1 matris 

Andra problem:

det (sI-A)

A är en 2x2 matris

I är en 2x2 diagonalmatris 

s är 1x2 matris som ska lösas ut

Det är tillståndsmodeller(Engelska: State-space) som jag håller på med.