Numerische-Methoden/MonteCarlo-Kotlin
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Multiple Runs

Browse Source
This commit is contained in:
Parric007
2024-05-16 11:56:38 +02:00
parent b5ed399b0c
commit 042ce5d141
1 changed files with 14 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
src/Main.kt
View File

@@ -96,18 +96,20 @@ fun main() {
val format = DecimalFormat("#,###.##") val format = DecimalFormat("#,###.##")
var pointsInside = 0.0 var pointsInside = 0.0
val samples = 20_000_000 val samples = 10_000_000
val bounds = Bounds(LinearBounds(lowerX.toFloat(), upperX.toFloat()), LinearBounds(getLowerLimit(), getUpperLimit())) val bounds = Bounds(LinearBounds(lowerX.toFloat(), upperX.toFloat()), LinearBounds(getLowerLimit(), getUpperLimit()))
var integral = 0.0 var integral = 0.0
var error = 0.0 //var error = 0.0
var result = 0.0 //var result = 0.0
val listOfIntegrals = mutableListOf<Double>();
val listOfErrors = mutableListOf<Double>();
val realValue = getRealIntegral(bounds.getX()) val realValue = getRealIntegral(bounds.getX())
var toTestPoint : Point var toTestPoint : Point
val elapsed = measureTimeMillis { val elapsed = measureTimeMillis {
for(i in 0..50){ for(i in 0..20){
for(ii in 0..samples) { for(ii in 0..samples) {
toTestPoint = bounds.getRandomPoint() toTestPoint = bounds.getRandomPoint()
if(getIsInside(toTestPoint)){ if(getIsInside(toTestPoint)){
@@ -117,15 +119,17 @@ fun main() {
integral = (pointsInside/samples) * bounds.area() integral = (pointsInside/samples) * bounds.area()
result += integral listOfErrors.add(abs(realValue-integral))
listOfIntegrals.add(integral)
integral = 0.0 integral = 0.0
pointsInside = 0.0 pointsInside = 0.0
} }
result /= 50.0 //result /= 50.0
error = abs(realValue-result) //error = abs(realValue-result)
} }
println("The approximated Integral of the function is: " + format.format(result)) println("The approximated Integral of the function is: " + format.format(listOfIntegrals.first()))
println("The real Integral of the function is: " + format.format(realValue)) println("The real Integral of the function is: " + format.format(realValue))
println("That's an error of " + format.format(error) + " or " + format.format((error/realValue)*100) + "%") println("That's an error of " + format.format(listOfErrors.first()) + " or " + format.format((listOfErrors.first()/realValue)*100) + "%")
println("And the whole thing took " + format.format(elapsed/1000.0) + " Seconds for 50 times" + format.format(samples) + " Samples") println("And the whole thing took " + format.format(elapsed/1000.0) + " Seconds for 20 times " + format.format(samples) + " Samples")
println(listOfErrors); println(listOfIntegrals)
} }