close all;
clear all;

t = linspace(-10,10,1000);
x0 = 1;
xd0 = 0;
a0 = x0;
a(1) = xd0;
a(2) = 0;
a(3) = -a0/(3*2);
x = a0*ones(1,1000);

for n=1:100
   a(n+3) = -a(n)/((n+3)*(n+2)); 
   x = x + a(n)*t.^n;
end

plot(t,x);
axis([-10 10 -3 3]);