t = linspace(-4,2,1000);
