8 replies to this topic

[Tritonio]

Here's a little program that tries to estimate π using a Monte Carlo method.

Be patient. :-) After an hour it should have reached 3.14.

0→I
0→O
Lbl 1
If Ran# ²+Ran# ²<1:Then
1+I→I
Else
1+O→O
IfEnd
Locate 1,1,(I÷(I+O))×4
Goto 1

Edited by Tritonio, 12 November 2014 - 02:00 AM.

[somebody1234]

Why did you make this?

(Just curious, it looks really inefficient)

[Tritonio]

For fun. It is one of the least efficient methods but it's soooo simple.

[Krtyski]

If you use Isz I instead of I+1->I, also use => instead of If / Else / IfEnd, then time to get result would be a bit shorter. please note do not start from I=0, O=0, but from I=1, O=1.

 

1→I:1→O
Lbl 1
IRan# ²+Ran# ²<1=>Isz I:Isz O
Locate 1,1,(I÷(I+O))×4
Goto 1

Edited by Krtyski, 15 November 2014 - 05:13 PM.

[Tritonio]

If you use Isz I instead of I+1->I, also use => instead of If / Else / IfEnd, then time to get result would be a bit shorter. please note do not start from I=0, O=0, but from I=1, O=1.

 

1→I:1→O
Lbl 1
IRan# ²+Ran# ²<1=>Isz I:Isz O
Locate 1,1,(I÷(I+O))×4
Goto 1

 

Indeed you can shorten it more. Although you made a mistake. You left lsz O running on every iteration this way so the locate line should be like this:

Locate 1,1,(I÷O)×4

 

Also why should I start from 1?

[Krtyski]

Oh, you are right.

 

In scripting 

A=>B:C

C is always processed. This is indeed wrong.

 

With 'Isz O' when O comes to zero (0), a command next to Isz O is skipped. Isz O is not equivalent to O+1->O. That's why O should not start with zero, but it should start with 1, 0.1, 0.001 so on...

Edited by Krtyski, 15 November 2014 - 07:23 PM.

[Krtyski]

I played around with following a but modified program;
 
1->I:1->O:1->D
Lbl 1
If Ran#²+Ran#²<1
Then Isz I
Else Isz O
IfEnd
Isz D
Locate 1,1,D
Locate 1,2,(I÷(O+O))x4

Getkey=0=>Goto 1

 

By hitting any key the calculation stops. I tried several times. When D is about 1500, the result was about between 3.13 to 3.15,

As the theory says accuracy is like this. It's interesting.

 

[Tritonio]

Oh, you are right.

 

In scripting 

A=>B:C

C is always processed. This is indeed wrong.

 

With 'Isz O' when O comes to zero (0), a command next to Isz O is skipped. Isz O is not equivalent to O+1->O. That's why O should not start with zero, but it should start with 1, 0.1, 0.001 so on...

 

I think it first increases the variable and then checks if it's zero. So it should work... I mean it works for me at least...

[Krtyski]

At last I enjoyed with following souce;

 

0->I:0->O:1->D

Lbl 1

Ran#²+Ran#²<1=>IszI

Isz O:Isz D

Locate 1,1,D

Locate 1,2,(I÷O)x4

 

Yes, the first time of Isz I gives I=1, so the jump never happen, as you said.

Please note the last Locate display (I÷O)X4, so this works correctly.