% Putzer-Algorithmus für Aufgabe 26.1 c, by Ute Feldmann
A=[2 0 0;2 1 -2;-1 0 2]
L=[1 2 2]%L=[2 2 1]%L=[2 1 2]
P0=eye(3)
P1=A-L(1)*P0
P2=(A-L(2)*eye(3))*P1
P3=(A-L(3)*eye(3))*P2
syms x real
w1=exp(L(1)*x)
w2=expand(exp(L(2)*x)*int(exp(-L(2)*x)*w1,0,x)) %x*exp(L(2)*x),%w2=(exp(L(1)*x)-exp(L(2)*x))/(L(1)-L(2))
w3=expand(exp(L(3)*x)*int(exp(-L(3)*x)*w2,0,x))
Y=P0*w1+P1*w2+P2*w3




