x=[];
for t = linspace(0, 30, 100)
  x = [x, springs(t,[1;1;1;2;3])];
endfor

plot(x')

% ;;; Local Variables: ***
% ;;; mode: octave ***
% ;;; End: ***
