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}}