library("kernel")
x = aseq(-1, 50, 0.04)
y = bqua2(x,0.5)
t = createdisplay(1,1)
show(t,1,1,x~y)

library("kernel")
library("plot")
n = 200
x =(1:(n))/(n)
x=(-x[(n):1])|0|x
y=bqua2(x,1)
y9=bqua2(x,0.9)
y8=bqua2(x,0.8)
y7=bqua2(x,0.7)
y6=bqua2(x,0.6)
y5=bqua2(x,0.5)
y4=bqua2(x,0.4)
y3=bqua2(x,0.3)
y2=bqua2(x,0.2)
y1=bqua2(x,0.1)
y0=bqua2(x,0)
xy=setmask(x~y,"line","thin")
xy9=setmask(paf(x~y9,x<0.9),"line","thin")
xy8=setmask(paf(x~y8,x<0.8),"line","thin")
xy7=setmask(paf(x~y7,x<0.7),"line","thin")
xy6=setmask(paf(x~y6,x<0.6),"line","thin")
xy5=setmask(paf(x~y5,x<0.5),"line","thin")
xy4=setmask(paf(x~y4,x<0.4),"line","thin")
xy3=setmask(paf(x~y3,x<0.3),"line","thin")
xy2=setmask(paf(x~y2,x<0.2),"line","thin")
xy1=setmask(paf(x~y1,x<0.1),"line","thin")
xy0=setmask(paf(x~y0,x<0),"line","thin")
axes=graxes(#(-1.2,1.2)~#(min(y0),max(y0)),"origin",0|0,"xticks",#(-1,-0.5,0.5,1),"yticks",#(min(y0)|max(y0)))
axesoff()
setsize(600,300)
di2 = createdisplay(1, 1)
show(di2, 1, 1, xy, xy9, xy8, xy7, xy6, xy5, xy4, xy3, xy2, xy1, xy0, axes)