function p3
    p3_do_plot(-5e-8, 5e-8, 1000)
    pause
end

function do_plot(a, b, n)
    x = linspace(a, b, n);
    y = 1-cos(x);
    plot(x,1-cos(x));
end

