Вставить в одну координатную плоскость несколько графиков функций

clc; close all; clear all;
a=0; % задаем начало интервала.
b=2; % задаем конец интервала.
h=0.5; % задаем шаг в интервале.
x = (a:h:b); % задаем значение x
y = sin(3*x)+cos(x).^2; % Значения у
N = length(x); % определяем количество точек.
xx = linspace(min(x),max(x),100); % формируем линейный массив размера 1x100 от минимального к максимальному.
% Создаем матрицу степеней х
for i = 1:N 
    for j = 1:N
        A(i,j) = x(i)^(N-j);
    end
end
B = y'; % вектор столбец правой части.
P = A\B; % считаем коэффициенты многочлена.
yy = polyval(P,xx); % вычисляем функцию в новых точках.
plot(x,y,'ob', x,y,'b--',xx,yy,'r') % вывод узлов интерполяции, исходной функции и графика полинома (интерполянты). 
%syms x
%ezplot(sin(3*x)+(cos(x))^2,[0 2])
axis tight % устанавливаем диапазоны координат по осям в соответствии с диапазонами изменения данных;
legend('Узлы интерполяции', 'Исходная функция','Интерполянта')
xlabel('x') % подписываем ось x.
ylabel('y') % подписываем ось y.
title('Полиномиальная интерполяция')

clc; close all; clear all;
syms xxx
 
a=0; % задаем начало интервала.
b=2; % задаем конец интервала.
h=0.5; % задаем шаг в интервале.
x = (a:h:b); % задаем значение x
y = sin(3*x)+cos(x).^2; % Значения у
N = length(x); % определяем количество точек.
xx = linspace(min(x),max(x),100); % формируем линейный массив размера 1x100 от минимального к максимальному.
% Создаем матрицу степеней х
for i = 1:N 
for j = 1:N
A(i,j) = x(i)^(N-j);
end
end
B = y'; % вектор столбец правой части.
P = A\B; % считаем коэффициенты многочлена.
yy = polyval(P,xx); % вычисляем функцию в новых точках.
plot(x,y,'ob',xx,yy,'r') % вывод узлов интерполяции, исходной функции и графика полинома (интерполянты). 
axis tight % устанавливаем диапазоны координат по осям в соответствии с диапазонами изменения данных;
legend('Узлы интерполяции', 'Интерполянта')
xlabel('x') % подписываем ось x.
ylabel('y') % подписываем ось y.
title('Полиномиальная интерполяция')
hold on
ezplot(sin(3*xxx)+(cos(xxx))^2,[a b]);