clear all
close all


load('HwPart3varyKpKdMSE');

p0 = linspace(1,10,10);
p1 = linspace(1,10,10);

figure (77)
surf(p0,p1,MeanSquare')
grid on;
% axis square;
title('Kp 1:1:10  Kd 5:5:50')
xlabel('Order of fit');
ylabel('Order of fit');
zlabel('Mean square error');


load('HwPart3varyKpKdMin');

p0 = linspace(5,50,10);
p1 = linspace(1,10,10);


figure (88)
surf(p1,p0,MinValueMatrix')
grid on;
% axis square;
title('Kp 1:1:10  Kd 5:5:50')
xlabel('Kp');
ylabel('Kd');
zlabel('Order of fit');