График (парабола + касательная)

clc; clear all; close all;
% зададим значение коэффициента а
a=5;
% зададим значение констатнты b
b=50;
% определим диапазон изменения переменной х
x=-2*b:0.1:2*b;
%зададим функцию кривой (параболы)
y=a.*sym('x').^2;
% вычисляем значение функции в точке b
Yb=subs(y,'x',b);
% касательной к функции в точке является производной функции
Ycasat=diff(a.*sym('x').^2);
% нахождение коэффициента качательной
K=subs(Ycasat,'x',b);
% Определим уравнение касательной
Yt=sym('y0+Ycasat*(x-x0)');
% подстановка коэффициентов
Yt=subs(Yt,'Ycasat',K);
Yt=subs(Yt,'x0',b);
Yt=subs(Yt,'y0',Yb);
%выводим график обеих функций
hold on
xlabel('X')
ylabel('Y')
grid on
ezplot(char(y),[-2*b 2*b])
ezplot(Yt, [-2*b 2*b])
plot(b,Yb,'or')
title('y=a*x^2 и касательная в точке b')
legend('a*x^2','касательная')
hold off

a = 2;
b = 3;
f = @(x)a*x.^2; % уравнение параболы
x0 = b;
y0 = f(x0);
g = @(x) 2*a*x0*(x - x0) + y0; % уравнение касательной к точке параболы с координатами (x0,y0)
diapason = [-2*x0 2*x0];
fplot(f, diapason)
Ylim = ylim;
hold on
fplot(g, diapason, 'g')
plot(x0,y0,'.r')
ylim(Ylim);
xlabel('x')
ylabel('y')
legend('{\itf(x)} = a * \itx^2','касательная в точке {\itx_0}','точка ({\itx_0},{\itf}({\itx_0}))')