#pdf(file="3.36plot.pdf")
par(mfrow=c(2,1))
par(pty="s")
x <- seq(-1,1, length=500)
mu <- 0
sigma <- 1
fx <- (1/sigma)*(63/4)*(((x-mu)/sigma)^6 - ((x-mu)/sigma)^8)
plot(x, fx, ylab="f(x)", type="l")
title("f(x) on -1 < x < 1")

x <- seq(-1,5,length=500)
#(a)
mu <- 0
sigma <- 1
fx <- (1/sigma)*(63/4)*(((x-mu)/sigma)^6 - ((x-mu)/sigma)^8)
plot(x, fx, ylab="f(x)", ylim=c(0,2), type="l")
#(b)
mu <- 3
sigma <- 1
fx <- (1/sigma)*(63/4)*(((x-mu)/sigma)^6 - ((x-mu)/sigma)^8)
lines(x, fx, lty=3)
#(c) 
mu <- 3
sigma <- 2
fx <- (1/sigma)*(63/4)*(((x-mu)/sigma)^6 - ((x-mu)/sigma)^8)
lines(x, fx, lty=4)
title("(a) solid, (b) dots, (c) dot dash, all on same axes, -1 < x < 5")
#dev.off()