x = [x0:0.1:x1];
y = x.^2 - 0.5;
plot(x,y);
print -dppm foo.ppm
quit

