alfa=Dinac.alfa.aw1(126)
eps=Dinac.eps.aw1

f=Dinac.fun(eps=eps,alfa=alfa)

n=100
x=c(rep(0,n/2),1,rep(0,n/2-1))

t=120
Gs.postscript("18-aw1-126.ps",4,5)
latox=c(0,n); latoy=c(0,t)
par(bg="yellow")
Gd.grafica(latox,latoy)
Dinac(f,x,t)
dev.off()