忍者ブログ

new3room関連ブログ

"new3の部屋"の補完ブログ。 興味有るものを右の”分類”から選んでね:お勧めは 工学系(電子、頭の体操)。 お笑系(落語、テレビ、頭の体操、もろもろ、AV)。 ローカルネタ(赤穂、食、気候、千種・観光)。 乗物(自転車、運動、交通)。修正履歴はホームページの更新履歴(管理人覚書)。 何か連絡あるなら、 >分類>雑記帳(4)>とりあえず何か言いたい場合ここにコメントしてok をクリックして、この記事の開いた中のコメント欄へどうぞ。 非公開機能停止、チャット機能運用するとパンクしそうなので停止。

10進BASIC

十進BASICのホームページ

上より10進basic言語を導入してみた。
インストールしないで、展開利用版。
実行すると、windowsが実行ブロックするという表示がでるが、詳細ボタンをおして実行する。
ヘルプが表示されないのは、展開ディレクトリトップにある
BASIC.chm
を、ダブルクリップ。開いたダイアログの実行を監視するオプションを外す。

知恵袋のパイに関する質疑に関連して、またモンテカルロ法によるΠの計算をしてみる。
3.14358457179229
 3.1435847858962 
400万回リピート、数十秒後の計算結果は3.144。まだ足りん?乱数の精度(というかランダム度というか)?

RANDOMIZE
FOR i=1 TO 4000000
   LET x=RND*2-1
   LET y=RND*2-1
   IF x^2+y^2 <1 THEN LET en=en+1 
   PRINT en/i*4
NEXT i
END
PR

コメント

お名前
タイトル
文字色
メールアドレス
URL
コメント
パスワード Vodafone絵文字 i-mode絵文字 Ezweb絵文字

モンテカルロ法によるΠの収束wo

10回の移動平均を実施しながら、Πの収束をみてみました。
100万回。
RANDOMIZE
DIM ref(10)
SET WINDOW 0,1000,0.1,-0.1
SET POINT STYLE 2
SET POINT COLOR 2
REM DRAW GRID0
DRAW AXES0
FOR i=1 TO 1000000
LET x=RND*2-1
LET y=RND*2-1
IF x^2+y^2 <1 THEN LET en=en+1
LET ref(MOD(i,10)+1)=en/i*4-PI
LET sumref=0
FOR j=1 TO 10
LET sumref=sumref+ref(j)
NEXT j
PRINT en/i*4,ref(MOD(i,10)+1),sumref
PLOT POINTS: i/1000,sumref
NEXT i
END

トラックバック

ブログ内検索

カレンダー

08 2024/09 10
S M T W T F S
1 2 3 5 6 7
8 9 10 11 12 13 14
15 16 17 18 19 20 21
22 23 24 25 26 27 28
29 30

ワイはこんな奴やで

愛称:
new3 :”しんさん”と呼んでネ
こっちも読んでね:
職業:
体重が気になる田舎者の遊び人
趣味:
電子・自転車・落語・パズル
自己紹介:
かに座B型で、"B型自分の取扱説明書"によると,”変人といわれるとなぜかうれしい”らしい。
自称”ぽち”。ご主人様に尻尾を振って、責任ないまま、生きて行きたい。

ご意見ご感想は

第三者非公開の場合、下部↓リンク欄の”new3の部屋(こっちも見てね)”を開いて、トップページの下部の、”関連リンク”ページの先頭部分にヒント。

  公開してよい連絡は↓コメント記入で

バーコード

P R