Como A é triangular superior, podemos tratar as componentes em cascata: primeiro as equações desacopladas (y3,y4) e depois usar essas soluções em y2 e, por fim, em y1.
Passo 2
Para dot y3=3y3 com y3(0)=1 e dot y4=4y4 com y4(0)=1 a solução é direta, pois dot y=lambda y dá