Sleeping Beauty problem: Difference between revisions

{{trans|Python}}

 

‘

Run the Sleeping Beauty Problem experiment `repetitions` times, checking to see

 

V CREDENCE = sleeping_beauty_experiment(1'000'000)

 

{{out}}

{{trans|Wren}}

 

wakings: 0

heads: 0

pc: sleepingBeauty 100000

 

{{out}}

==={{header|BASIC256}}===

{{trans|FreeBASIC}}

iteraciones = 1000000

cara = 0

print "Percentage probability of heads on waking = "; (cara/dormir*100); "%"

end

{{out}}

<pre>

 

==={{header|FreeBASIC}}===

Const iteraciones = 1000000

Randomize Timer

Print using "Percentage probability of heads on waking = ###.######%"; (cara/dormir*100)'; "%"

Sleep

{{out}}

<pre>

==={{header|GW-BASIC}}===

In this simulation, Sleeping Beauty flips a coin of her own.

20 MONDAY = 0 : TUESDAY = 1

30 HEADS = 0 : TAILS = 1

250 IF GUESS = HEADS THEN WHEADS = WHEADS + 1 ELSE CTAILS = CTAILS + 1

260 NEXT DAY

{{out}}<pre>

Sleeping Beauty was put through this experiment 300000 times.

==={{header|Yabasic}}===

{{trans|FreeBASIC}}

iteraciones = 1000000

cara = 0

print "Percentage probability of heads on waking = ", (cara/dormir*100), "%"

end

{{out}}

<pre>

 

=={{header|C++}}==

#include <random>

 

std::cout << "Sleeping Beauty should estimate a credence of: "

<< double(heads) / wakenings << '\n';

 

{{out}}

 

=={{header|CLU}}==

% to use the random number generator.

 

chance: real := experiment$run(N)

stream$putl(po, f_form(chance, 1, 6))

{{out}}

<pre>Chance of waking up with heads: 0.333758</pre>

{{trans|Swift}}

 

var heads = 0

var wakenings = 0

}

print("Wakenings over \(experiments) experiments: \(wakenings)")

 

=={{header|Excel}}==

 

{{Works with | Office 365 betas 2021}}

=LAMBDA(n,

LET(

)

)

{{Out}}

 

 

=={{header|F_Sharp|F#}}==

// Sleeping Beauty: Nigel Galloway. May 16th., 2021

let heads,woken=let n=System.Random() in {1..1000}|>Seq.fold(fun(h,w) g->match n.Next(2) with 0->(h+1,w+1) |_->(h,w+2))(0,0)

printfn "During 1000 tosses Sleeping Beauty woke %d times, %d times the toss was heads. %.0f%% of times heads had been tossed when she awoke" woken heads (100.0*float(heads)/float(woken))

{{out}}

<pre>

=={{header|Factor}}==

{{works with|Factor|0.99 2021-02-05}}

 

: sleeping ( n -- heads wakenings )

"Wakenings over 1,000,000 experiments: " write

1e6 sleeping dup . /f

{{out}}

<pre>

{{trans|Wren}}

{{libheader|Go-rcu}}

 

import (

pc := sleepingBeauty(1e6)

fmt.Printf("Percentage probability of heads on waking = %f%%\n", pc)

 

{{out}}

 

=={{header|Haskell}}==

import System.Random (randomIO)

import Control.Monad (replicateM)

let (Sum w, Sum h) = foldMap anExperiment tosses

let ratio = fromIntegral h / fromIntegral w

 

<pre>*Main> main

=={{header|J}}==

Simulation code:

monday=. ?2

if. -. monday do.

<monday

end.

 

#sample NB. number of experiments

1000000

0.4996

sample+&#;sample NB. total number of awakenings

 

It's probably worth noting here that the number of heads divided by the number of awakenings would be about 0.3 -- but Sleeping Beauty was not asked to guess whether the coin was heads on Wednesday.

 

=={{header|Julia}}==

Run the Sleeping Beauty Problem experiment `repetitions` times, checking to see

how often we had heads on waking Sleeping Beauty.

CREDENCE = sleeping_beauty_experiment(1_000_000)

println("Results of experiment: Sleeping Beauty should estimate a credence of: ", CREDENCE)

Wakenings over 1000000 experiments: 1499534

Results of experiment: Sleeping Beauty should estimate a credence of: 0.33374768428058316

 

=={{header|Mathematica}}/{{header|Wolfram Language}}==

SleepingBeautyExperiment[reps_Integer] := Module[{gotheadsonwaking, wakenings, coinresult},

gotheadsonwaking = 0;

]

out = N@SleepingBeautyExperiment[10^6];

{{out}}

<pre>Wakenings over 1000000 experiments: 1499714

 

=={{header|Nim}}==

 

const N = 1_000_000

 

echo "Wakenings over ", N, " experiments: ", wakenings

 

{{out}}

=={{header|Pascal}}==

{{trans|Phix}}

program sleepBeau;

uses

end;

writeln(Format(fmt,[iterations,wakings,heads/wakings*100]));

{{out}}

<pre>Wakings over 1000000 repetitions = 1499741

 

=={{header|Perl}}==

use warnings;

 

 

my $trials = 1_000_000;

{{out}}

<pre>Wakenings over 1000000 experiments: 1499816

 

=={{header|Phix}}==

<span style="color: #008080;">constant</span> <span style="color: #000000;">iterations</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">1_000_000</span><span style="color: #0000FF;">,</span>

<span style="color: #000000;">fmt</span> <span style="color: #0000FF;">=</span> <span style="color: #008000;">"""

<span style="color: #008080;">end</span> <span style="color: #008080;">for</span>

<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">fmt</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">iterations</span><span style="color: #0000FF;">,</span><span style="color: #000000;">wakings</span><span style="color: #0000FF;">,</span><span style="color: #000000;">heads</span><span style="color: #0000FF;">/</span><span style="color: #000000;">wakings</span><span style="color: #0000FF;">*</span><span style="color: #000000;">100</span><span style="color: #0000FF;">})</span>

{{out}}

<small>(You'll get the exact result less than 1% of the time!!)</small>

=={{header|Python}}==

===Procedural===

 

def sleeping_beauty_experiment(repetitions):

CREDENCE = sleeping_beauty_experiment(1_000_000)

print("Results of experiment: Sleeping Beauty should estimate a credence of:", CREDENCE)

Wakenings over 1000000 experiments: 1499765

Results of experiment: Sleeping Beauty should estimate a credence of: 0.333542254953276

===Functional===

 

 

from random import choice

if __name__ == '__main__':

main()

{{Out}}

<pre>1500188 wakenings over 1000000 experiments.

=={{header|Quackery}}==

 

 

[ say "Number of trials: "

10 round vulgar$ echo$ cr ] is trials ( n --> n/d )

 

 

{{out}}

=={{header|R}}==

There's nothing complicated here. Pretty much every language that resembles C is going to use virtually the same solution.

{

wakeCount <- headCount <- 0

headCount/wakeCount

}

{{out}}

<pre>[1] 0.3335838</pre>

=={{header|Raku}}==

 

my $gotheadsonwaking = 0;

my $wakenings = 0;

}

 

{{out}}

<pre>Wakenings over 1000000 experiments: 1500040

 

=={{header|Red}}==

 

experiments: 1'000'000

]

print ["Awakenings over" experiments "experiments:" awakenings]

{{out}}

<pre>

=={{header|REXX}}==

When using Regina REXX, &nbsp; the seed specified &nbsp; (for '''random''') &nbsp; was &nbsp; '''46'''.

parse arg n seed . /*obtain optional arguments from the CL*/

if n=='' | n=="," then n= 1000000 /*Not specified? Then use the default.*/

exit 0 /*stick a fork in it, we're all done. */

/*──────────────────────────────────────────────────────────────────────────────────────*/

{{out|output|text=&nbsp; when using the input of: &nbsp; &nbsp; <tt> , &nbsp; 46 </tt>}}

<pre>

=={{header|Ruby}}==

{{trans|Python}}

coin = [:heads, :tails]

gotheadsonwaking = 0

puts "Results of experiment: Sleeping Beauty should estimate

a credence of: #{sleeping_beauty_experiment(1_000_000)}"

{{out}}

<pre>Wakenings over 1000000 experiments: 1499604

</pre>

=={{header|Swift}}==

var heads = 0

var wakenings = 0

}

print("Wakenings over \(experiments) experiments: \(wakenings)")

 

{{out}}

=={{header|Vlang}}==

{{trans|Go}}

import rand.seed

 

pc := sleeping_beauty(1000000)

println("Percentage probability of heads on waking = $pc%")

 

{{out}}

=={{header|Wren}}==

{{libheader|Wren-fmt}}

import "/fmt" for Fmt

 

 

var pc = sleepingBeauty.call(1e6)

 

{{out}}