clc
clear

x = (-5, 5, 1000);
y = 2*x.*sin(x).*cos(x) + x.^3;

plot(x, y, 'LineWidth', 2);
title('y = 2x sin(x) cos(x) + x^3');
xlabel('x');
ylabel('y');
grid on;