From 40320dc5b69818887df867bacf51c34ebac6e1cd Mon Sep 17 00:00:00 2001 From: Julian Schacher Date: Fri, 3 Aug 2018 00:07:36 +0200 Subject: [PATCH 1/3] Print pi estimate. --- contents/monte_carlo_integration/code/julia/monte_carlo.jl | 1 + 1 file changed, 1 insertion(+) diff --git a/contents/monte_carlo_integration/code/julia/monte_carlo.jl b/contents/monte_carlo_integration/code/julia/monte_carlo.jl index fee58cdcd..65289ccbc 100644 --- a/contents/monte_carlo_integration/code/julia/monte_carlo.jl +++ b/contents/monte_carlo_integration/code/julia/monte_carlo.jl @@ -24,6 +24,7 @@ function monte_carlo(n::Int64) # are only using the upper quadrant and the unit circle, so we can use # 4*pi_count/n instead pi_estimate = 4*pi_count/n + println("The pi estimate is: ", pi_estimate) println("Percent error is: ", signif(100*(pi - pi_estimate)/pi, 3), " %") end From 87537080b4f38651aeb196709dcbd35effb8605e Mon Sep 17 00:00:00 2001 From: Julian Schacher Date: Fri, 3 Aug 2018 00:07:56 +0200 Subject: [PATCH 2/3] Fix percent error calculation. --- contents/monte_carlo_integration/code/julia/monte_carlo.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/contents/monte_carlo_integration/code/julia/monte_carlo.jl b/contents/monte_carlo_integration/code/julia/monte_carlo.jl index 65289ccbc..181976634 100644 --- a/contents/monte_carlo_integration/code/julia/monte_carlo.jl +++ b/contents/monte_carlo_integration/code/julia/monte_carlo.jl @@ -25,7 +25,7 @@ function monte_carlo(n::Int64) # 4*pi_count/n instead pi_estimate = 4*pi_count/n println("The pi estimate is: ", pi_estimate) - println("Percent error is: ", signif(100*(pi - pi_estimate)/pi, 3), " %") + println("Percent error is: ", signif(100*(pi_estimate - pi)/pi, 3), " %") end monte_carlo(10000000) From e5fe3ad42d37f645a67ad85ad78399a300ef3b71 Mon Sep 17 00:00:00 2001 From: Julian Schacher Date: Fri, 3 Aug 2018 00:37:40 +0200 Subject: [PATCH 3/3] Make the percent error always be a positive value, by using the absolute value of a critical calculation step. Also space out the calculation of the error a bit. --- contents/monte_carlo_integration/code/julia/monte_carlo.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/contents/monte_carlo_integration/code/julia/monte_carlo.jl b/contents/monte_carlo_integration/code/julia/monte_carlo.jl index 181976634..51b9335a8 100644 --- a/contents/monte_carlo_integration/code/julia/monte_carlo.jl +++ b/contents/monte_carlo_integration/code/julia/monte_carlo.jl @@ -25,7 +25,7 @@ function monte_carlo(n::Int64) # 4*pi_count/n instead pi_estimate = 4*pi_count/n println("The pi estimate is: ", pi_estimate) - println("Percent error is: ", signif(100*(pi_estimate - pi)/pi, 3), " %") + println("Percent error is: ", signif(100 * abs(pi_estimate - pi) / pi, 3), " %") end monte_carlo(10000000)