Schrick-Noah_Project-2.R

@@ -111,6 +111,7 @@ fun.2c <- function(x) {3-2*exp(-x)}
 fun.2c.fx <- function(x) {x+2*exp(-x)-3}
 fun.2c.string <- "x+2*exp(-x)-3"
 fun.2c.estimate <- findZeroRelax(fun.2c,-.5)
+fun.2c.estimate[1]
 # Find the other root with bisection
 fun.2c.bisect <- findZeroBisect(fun.2c,-1,0,1e-6)
 # Plot 2c
@@ -211,10 +212,7 @@ ggplot(data.frame(x=seq(1,80,.1)), aes(x)) +
   geom_hline(aes(yintercept=0, col = "y=0"), show.legend=TRUE)+
   geom_vline(aes(xintercept=func.five.estimate[1], 
                  col = "bisection estimate"), show.legend=TRUE)+
-  ggtitle("zeroth-order dimensional 
-          perturbation theory approximation
-          for the chemical potential of the
-          Gross-Pitaevskii equation in d dimensions") +
+  ggtitle("zeroth-order dimensional perturbation theory approximation for the chemical potential of the Gross-Pitaevskii equation in d dimensions") +
   xlab("x") + ylab("y") +
   theme(text = element_text(size=20), plot.title = element_text(hjust = 0.5))
 
@@ -239,9 +237,6 @@ ggplot(data.frame(x=seq(5,80,.1)), aes(x)) +
   geom_hline(aes(yintercept=0, col = "y=0"), show.legend=TRUE)+
   geom_vline(aes(xintercept=func.five.estimate[1], 
                  col = "bisection estimate"), show.legend=TRUE)+
-  ggtitle("zeroth-order dimensional 
-          perturbation theory approximation
-          for the chemical potential of the
-          Gross-Pitaevskii equation in d dimensions") +
+  ggtitle("zeroth-order dimensional perturbation theory approximation for the chemical potential of the Gross-Pitaevskii equation in d dimensions") +
   xlab("x") + ylab("y") +
   theme(text = element_text(size=20), plot.title = element_text(hjust = 0.5))