Delta = 0.7;

[t,x,y] = sim('L17_interpolation_sim', 4);
vin = y(:, 1) / Delta;
y   = y(:, 2:3);
 
figure(1);
plot(vin, y(:,1), vin, y(:,2), vin, vin*0, ':');
axis([ 0 3 -0.6 0.6 ]);
xlabel('Vin / \Delta');
ylabel('Preamp Output');
legend('A_1', 'A_2', -1);
grid on;
fixfig;

figure(2);
subplot(2, 1, 1);
plot(vin, y(:,1), vin, -y(:,1), vin, y(:,2), vin, -y(:,2), vin, vin*0, 'k:');
axis([ 0 3 -0.6 0.6 ]);
ylabel('Preamp Output');
legend('A_1', '-A_1', 'A_2', '-A_2', 1);
fixfig;
subplot(2, 1, 2);
plot(vin, y(:,1)+y(:,2), vin, vin*0, 'k:');
axis([ 0 3 -0.6 0.6 ]);
xlabel('Vin / \Delta');
ylabel('A_1+A_2');
legend('A_1+A_2', 1);
fixfig;

Delta = 1.3;

[t,x,y] = sim('L17_interpolation_sim', 6);
vin = y(:, 1) / Delta;
y   = y(:, 2:3);
 
figure(3);
subplot(2, 1, 1);
plot(vin, y(:,1), vin, -y(:,1), vin, y(:,2), vin, -y(:,2), vin, vin*0, 'k:');
axis([ 0 3 -0.6 0.6 ]);
ylabel('Preamp Output');
legend('A_1', '-A_1', 'A_2', '-A_2', 1);
fixfig;
subplot(2, 1, 2);
plot(vin, y(:,1)+y(:,2), vin, vin*0, 'k:');
axis([ 0 3 -0.6 0.6 ]);
xlabel('Vin / \Delta');
ylabel('A_1+A_2');
legend('A_1+A_2', 1);
fixfig;