clear all; close all;
N  = 100;
s0 = linspace(-1,0,N);
s1 = 1/2;

perr = zeros(1,N);
for i = 1:N
    perr(i) = p3(s0(i),s1);
end
plot(s0,perr);