tu as raison et sinon ce que tu peux faire avec ton système :
si F = a1(t).E1(t) + ... + an(t).En(t) : tu veux résoudre donc F' = PF + B
où F = (f,f',...f^(n-1)) la solution que tu cherches et Ei = (ei,...ei^(n-1)) avec les ei la base de solution (Ei' = PEi), P est la matrice associée à ton système et B = (g,0,0,..0)
tu remarques que F' = a1'E1 + ... + an'.En + (a1.E1'(t) + ... + an.En'(t))
= a1'E1 + ... + an'.En + (a1.PE1 + ... + an.PEn)
= a1'E1 + ... + an'.En + P.F
Donc a1'E1 + ... + an'.En = (g,0,0,...,0)
ce qui te donne avec les formules de cramer les valeurs de ai'(t) donc ai(t) par intégration