[t, x] = ode45('test', [0 2*pi], 3);
plot(t, x, 'linewidth', 3)
title('The solution to x'' = cos(t)/(2x-2), with x(0) = 3.')
xlabel('t')
ylabel('x')
grid
