function linearsys
    A=[2,-1,1;1,1,-1;1,-1,1];
    %t=0:0.01:2;
    [T,Y]=ode45(@ff,[0,2],[-1,-1,1]);
    plot(T,Y(:,1),'-',T,Y(:,2),'-',T,Y(:,3),'.')
    axis([0,2,-40,40])
    function z=ff(t,y)
        z=A*y+[-2;-2;3]*exp(t);
    end

end

