Как сохранить 3 графика в Matlab-е?

clc; close all; clear all; 
   %syms X             
a=-10; % задаем начало интервала.
b=2; % задаем конец интервала.
h=1; % задаем шаг в интервале.
x=(a:h:b); % задаем значение x
y=x.^4+2*x.^3+7*x.^2-12; % Значения у
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); % вычисляем функцию в новых точках.
 
%% Выводим график исходной функции и узлов интерполяции
figure 
   %Y=ezplot(X.^4+2*X.^3+7*X.^2-12,[a b]); hold on
plot(x,y,'b-',x,y,'or');
axis tight
legend('Исходная функция','Узлы интреполяции')
xlabel('x') % подписываем ось x.
ylabel('y') % подписываем ось y.
 
%% Выводим график интерполянты и узлов интерполяции
figure 
plot(xx,yy,'g',x,y,'or'); % вывод узлов интерполяции, исходной функции и графика полинома (интерполянты). 
axis tight % устанавливаем диапазоны координат по осям в соответствии с диапазонами изменения данных.
legend('Интерполянта','Узлы интерполяции')
xlabel('x') % подписываем ось x.
ylabel('y') % подписываем ось y.
 
%% Выведем погрешность в виде графика для наибольшей наглядности.
yz=polyval(P,x);
Z=y-yz; % расчитываем погрешность.
figure
plot(x,Z,'b-');  % выводим в графическом виде.
axis tight
legend('Погрешность',-1)
m=max(Z);
disp(m);

cla

cla (handles.axes1,'reset') % стираем нарисованные графики. 
cla (handles.axes2,'reset') 
cla (handles.axes3,'reset')