P0 <- c(11.5, 7.8, 1.1, 5.6)
x <- seq(-11.5/2, 11.5/2, len=2000)
y1 <- NRGE(P=P0, x=x)
y2 <- -NRGE(P=P0, x=x)
dev.new()
plot(x, y1, cex.lab=1.5, cex.axis=1.5, type="l",
col=4, ylim=c(-4, 4), asp=1,
xlab=expression(italic(x)), ylab=expression(italic(y)) )
lines(x, y2, col=2)
graphics.off()
Run the code above in your browser using DataLab