## BEGIN OF FUNCTION DEFFINITION ##

RamSpecEval<- function(X0, Y0, FF, NC) {
Yi<- lowess(X0, Y0,f=FF)$y
nx<-length(X0)
for(i in 1:NC) {
for(j in 1:nx) { Yi[j]<- min(Y0[j],Yi[j])
	              }	             
	           Yi<- lowess(X0, Yi,f=FF)$y
                }
Yi<- Y0-Yi
cat('.')  
Yi
}
## END OF FUNCTION DEFFINITION ##

LIST<-dir(patt='R040114')
cols<-c('black','red','green', 'blue', 'purple', 'orange', 'brown', 'yellow')
n<-length(LIST)
quartz(width=11, height=5)
plot(c(0,1800),c(0,1200), typ='n', main="Standards from the RRUF Database")
for(i in 1:n){XX<-read.table(LIST[i], sep=',')
	          attach(XX)
	          V2<- RamSpecEval(V1, V2, 0.1, 10 )
	          mv2<-max(V2)
	          lines(V1,25*i+1000*V2/mv2,col=cols[i])
	          text(700,1100-70*i, LIST[i], col=cols[i])
	          rm(V2)
	          detach(XX)
	          }
rm(LIST, XX, mv2)
	          
LIST<-dir(patt='R050307')
m<-n+2
n<-length(LIST)
for(i in 1:n){XX<-read.table(LIST[i], sep=',')
	          attach(XX)
	          V2<- RamSpecEval(V1, V2, 0.1, 10 )
	          mv2<-max(V2)
	          lines(V1,25*(m+i)+1000*V2/mv2,col=cols[i], lwd=2)
	          text(1500,1100-70*i, LIST[i], col=cols[i], cex=1.2)
	          rm(V2)
	          detach(XX)
	          }