diff --git a/server.r b/server.r
new file mode 100644
index 0000000000000000000000000000000000000000..93e38ecf45ec4de5499d9f66fe9c0f7c47286826
--- /dev/null
+++ b/server.r
@@ -0,0 +1,57 @@
+
+# Define a server for the Shiny app
+function(input, output) {
+  
+  
+
+  # Fill in the spot we created for a plot
+  output$Plot <- renderPlot({
+
+	# mu=seq(0,1,length=100)
+	
+	# mypow1=ifelse(mu==input$mu0,input$myalpha1,
+	# 1-pnorm((sqrt(input$n)/input$sigma)*(input$mu0-mu)+qnorm(1-input$myalpha1)))
+    # plot(mu,mypow1,type="l",lwd=2,ylim=c(0,1),cex.lab=2.2,xlab=expression(mu-mu[0]),
+     # main="",ylab=expression(pi[n](mu)),cex=1.4, cex.axis=2,bty="n",
+	 # xlim=c(0,1),col="black")
+
+disable("mu")
+	 
+myalpha1=0.01
+mu0=1
+sigma=input$sigma
+n=input$n
+mu=seq(0,2,length=100)
+mypow1=ifelse(mu==mu0,myalpha1,1-pnorm((sqrt(n)/sigma)*(mu0-mu)+qnorm(1-myalpha1)))
+par(mar=c(4,5.5, 0.5, 0.5))	
+plot(mu-mu0,mypow1,type="l",lwd=2,ylim=c(0,1),cex.lab=2.2,xlab=expression(mu*-~mu[0]),
+     main="",ylab=expression(pi[n](mu)),cex=1.4, cex.axis=2,bty="n",
+	 xlim=c(-1,1),col="black",xaxt="n")
+axis(1, c(-1,0,1),labels=c(-1,0,1),cex.axis=2, tck=-.02)
+
+	 
+myalpha2=0.05
+mypow2=ifelse(mu==mu0,myalpha2,1-pnorm((sqrt(n)/sigma)*(mu0-mu)+qnorm(1-myalpha2))) 
+lines(mu-mu0,mypow2,cex.lab=2.2,xlab="",ylab="",col=c("red"),cex=1.4, 
+yaxt="n", xaxt="n",bty="n",type="l",lty=2,lwd=2)
+
+myalpha3=0.1
+mypow3=ifelse(mu==mu0,myalpha3,1-pnorm((sqrt(n)/sigma)*(mu0-mu)+qnorm(1-myalpha3))) 
+lines(mu-mu0,mypow3,cex.lab=2.2,xlab="",ylab="",col=c("blue"),cex=1.4, 
+yaxt="n", xaxt="n",bty="n",type="l",lty=3,lwd=3)
+	 
+legend("topleft",lty=c(NA,1,2,3),col=c("black","black","red","blue"),bty="n",
+cex=2,lwd=3,legend=c(
+as.expression(bquote(n*"="*.(n))),
+as.expression(bquote(alpha*"="*.(format(myalpha1, digits = 2)))),
+as.expression(bquote(alpha*"="*.(format(myalpha2, digits = 2)))),
+as.expression(bquote(alpha*"="*.(format(myalpha3, digits = 2))))))
+	 
+	 
+	 
+  })
+}
+
+
+
+