Par7 <- c(1.124, 14.86, 49.43)
phi1 <- seq(0, pi, len=2000)
g1 <- DSGE(P=Par7, phi=phi1)
dev.new()
plot(phi1, g1, type="l", col=4, cex.lab=1.5, cex.axis=1.5,
xlab=expression(italic(varphi)),
ylab=expression(paste(italic(g), "(", italic(varphi), ")", sep="")))
graphics.off()
Run the code above in your browser using DataLab