Вычисление производной и построение графика

x = sym('x'); %объявление символьной переменной x
a = sym('a'); %объявление символьной переменной a
y = a.*x.^2;
dy = diff (y); %вычисление производной в символьном виде
disp(dy); %выведение производной на экран
x = (-15:0.1:15); %задаем диапазон значений x
a = (-15:0.1:15); %задаем диапазон значений a
plot(x, subs(dy),'k'); %график проивзодной 
hold on %объединение графиков
plot(x, subs(y)); %график функции 
axis( [-15, 15,-15, 15] ) %установка диапазонов координат по осям х и у для графика
grid on %включение сетки 
xlabel('ось Ох'); %подпись оси Ox
ylabel('ось Оy'); %подпись оси Oy
title('график функции y = a*x^2 и график ее производной'); 
legend('график производной','график функции'); %создание легенды

% a = (-15:0.1:15); %задаем диапазон значений a
a=1; % или какое-нибудь другое число