myprojects-hg-reborn

Fork

rm(list=ls())
lambda<-read.table("lambda-value",header=FALSE)
lambda<-lambda[ ,1]
mass<-read.table("mass--values",header=FALSE)
mass<-mass[ ,1]
lambda<- -lambda #NB: the quantity I saved as lambda is as a matter of fact 
                 # the energy -Lambda; that is why I have to change its sign
pdf("mass-vs-lambda.pdf")
plot(lambda,mass,"l",lwd=2.5,col=300,xlab="Lambda",ylab="Mass")
dev.off()
energycentr<-read.table("energy-centr",header=FALSE )
 energycentr<-energycentr[ ,1]
 kincentr<-read.table("central-kine",header=FALSE )
 kincentr<-kincentr[ ,1]
 potcentr<-read.table("central-pote",header=FALSE )
 potcentr<-potcentr[ ,1]
 pdf("site-energy-vs-lambda.pdf")
 plot(lambda,energycentr,"l",lwd=2.5,col=300,xlab="Lambda",ylab="Energy Central Well",
 ylim=range(c(-6,6)))
 lines(lambda,potcentr,lwd=2.5,col="black")
 lines(lambda,kincentr,lwd=2.5,col="orange")
 dev.off()

tkin<-read.table("total-kineti",header=FALSE )
 tkin<-tkin[ ,1]
 
 tpot<-read.table("total-potent",header=FALSE )
 tpot<-tpot[ ,1]
 
 pdf("energy-vs-mass.pdf")
 # since I now want to plot the energy, I need to change again the sign of lambda
 # - (- lambda)=lambda, but what I originally saved as lambda was the energy.
plot(mass,-lambda,"l",lwd=2.5,col=300,xlab="Mass",ylab="Energy",ylim=range(c(-6,5)))
lines(mass,tkin,"l",lwd=2.5,col="black")
lines(mass,tpot,"l",lwd=2.5,col="orange")
dev.off()