clear;m=1;c=1;k=10;
x=roots([m,c,k]);
t=0:0.1:20;
y=1/(x(1)-x(2))*(exp(x(1)*t)-exp(x(2)*t));
plot(t,y);grid('on');