Your approach is wrong for the question. It asks for numeric verification.
So take x = linspace() a range notably wider than [-pi, pi] because you need several periods in order to verify periodic behaviour.
Then evaluate y = sin(x) + 2*cos(x) for that numeric x.
Now, do numeric differentiation,
and compare d2y + dy to 0