| if (par1 == '0') bw <- 'nrd0'if (par1 != '0') bw <- as.numeric(par1)
 bitmap(file='density1.png')
 mydensity1<-density(x,bw=bw,kernel='gaussian',na.rm=TRUE)
 plot(mydensity1,main='Gaussian Kernel',xlab=xlab,ylab=ylab)
 grid()
 dev.off()
 mydensity1
 bitmap(file='density2.png')
 mydensity2<-density(x,bw=bw,kernel='epanechnikov',na.rm=TRUE)
 plot(mydensity2,main='Epanechnikov Kernel',xlab=xlab,ylab=ylab)
 grid()
 dev.off()
 bitmap(file='density3.png')
 mydensity3<-density(x,bw=bw,kernel='rectangular',na.rm=TRUE)
 plot(mydensity3,main='Rectangular Kernel',xlab=xlab,ylab=ylab)
 grid()
 dev.off()
 bitmap(file='density4.png')
 mydensity4<-density(x,bw=bw,kernel='triangular',na.rm=TRUE)
 plot(mydensity4,main='Triangular Kernel',xlab=xlab,ylab=ylab)
 grid()
 dev.off()
 bitmap(file='density5.png')
 mydensity5<-density(x,bw=bw,kernel='biweight',na.rm=TRUE)
 plot(mydensity5,main='Biweight Kernel',xlab=xlab,ylab=ylab)
 grid()
 dev.off()
 bitmap(file='density6.png')
 mydensity6<-density(x,bw=bw,kernel='cosine',na.rm=TRUE)
 plot(mydensity6,main='Cosine Kernel',xlab=xlab,ylab=ylab)
 grid()
 dev.off()
 bitmap(file='density7.png')
 mydensity7<-density(x,bw=bw,kernel='optcosine',na.rm=TRUE)
 plot(mydensity7,main='Optcosine Kernel',xlab=xlab,ylab=ylab)
 grid()
 dev.off()
 load(file='createtable')
 a<-table.start()
 a<-table.row.start(a)
 a<-table.element(a,'Properties of Density Trace',2,TRUE)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'Bandwidth',header=TRUE)
 a<-table.element(a,mydensity1$bw)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'#Observations',header=TRUE)
 a<-table.element(a,mydensity1$n)
 a<-table.row.end(a)
 a<-table.end(a)
 table.save(a,file='mytable.tab')
 
 |