roots([1 5 -5 1])
m=ans
etothem = exp(m)
A = [etothem';  m' .* etothem'; (m').^2 .* etothem']
b = [0; 0; 2]
Bbbb = [0 0 2]
b = [0 0 2]'
c=A\b
t=1:.01:3
t=1:.01:3;
c
plot(t,c(1)*exp(m(1)*t)+c(2)*exp(m(2)*t)+c(3)*exp(m(3)*t))
t=3; c(1)*exp(m(1)*t)+c(2)*exp(m(2)*t)+c(3)*exp(m(3)*t)
c * exp(m*t)
c * exp(m'*t)
c' * exp(m*t)