高機能電卓の情報

2016年1月10日日曜日

賀正 2016年

何とも素っ気ないタイトルで申し訳ない。
明けましてお目出度いのはテメェのオツムのみだと困ってしまいますが、それでも鄙びた辺りから、か細い声でも謹んで新年を賀する所です。
松の内も終わりまして、本blogも少しずつ作業を始めようと思います。

正月早々から余り大ネタもないので、最近考えている事を「年賀タオル」の様にポツポツと。

1. HP50GのMultiple Equation Solvreがスゴかった

先日、HP50GのEquation Library を使ってみたのですが、結構面白いんですヨ。
一部の関数電卓にもある、ビルトインされた「公式集」の様なものではありますが、グラフ電卓の強みを活かして、図が表示される項目があります (一部だけではありますが)。
メニューから項目を選んでやると、登録された式を元にしたSolver環境が動き、そのまま数値解を計算出来るという便利ツールなのですが、どうも、50GのMultiple Equation Solver (MES)を利用したものらしい。Users guide を見ると、ヘロンの公式、余弦定理と複数の式を登録してMESを使う例が出ております。

Using the Multiple Equation Solver (MES) - HP50G Users Guide
(page 7-9)

大まかなMESの使い方は、

+ 方程式リストをEQ、SolverのタイトルをTITLE、方程式の独立変数リストをLVARI という変数に記憶させておき、
+ MINIT, MITM, MSOLVR の順にコマンドを実行する

という感じ。コマンド実行時には、スタックに必要なパラメタを積んでおくなどの準備がありますが、詳細はUsers guideを御山椒。

MESの何がスゴイかというと、複数の方程式を登録して、それらを元に数値解を自動で計算する所もありますが、コマンドを複数使うだけでSolver環境が使えるので、プログラムとして構成する事で「専用のSolverプログラム」に出来るのです。
最近の関数電卓には、方程式1本のSolverが用意されていますが、50GのMESは、更に先を行っていたのでした。

先日、JulyさんAmazonの所のHP Primeは在庫が無くなってしまいましたが、50Gは未だ11台残っております。気になる方はお買い求めを。

2. 電卓内蔵のPeriodic Table考

HP50Gには周期表が用意されております。最近の電卓でも、一部には内蔵されているものがあります。
CASIOのfx-J900/700では、118番目の元素まで記録されているとか。

先日、113番目の元素について、日本の研究者が3度に渡り生成する事に成功した事から、命名権を取得したという話題がありました。
HP50Gの周期表は、作られた時期が古く、106番目までの元素までしか収録されておりません。名前も、104,105,106 については「Un-nil Quadium」、「Un-nil Pentium」、「Un-nil Hexium」と、暫定名であります。
これは、ベースとなっているのがHP48SX向けに用意されたアプリケーションで、1995年頃に作成されているので、古いのも仕方の無い所です。
fx-J900/700の周期表は十分に新しく、113番目の元素も収録されているのですが、命名権を取得したという段階なので、未だ名前が決まった訳ではなく、コレばかりはどうしようもありません。その内、113番元素の名前も決まるのでしょうから、その時になったら周期表も改変されるのでしょうか。まあ、市中在庫が結構ある上、ファームウェアの改変も必要であり、簡単には進まないかも。
そこで、新しい高機能電卓を出して、周期表のアップデートを図るという手を期待したい。
最近の電子辞書では、コンテンツの追加やアップデートが出来る製品があると言います。また、ファームウェアのアップデートが可能な様に、マスクROMからFlashメモリにした電卓製品も多く普及して居ります。新しい高機能電卓では、ファームウェアのアップデート機能によって色々な改良が為され、永く使いつづける事の出来る製品である事を期待したいものです。

3. CASIOの「スマートウォッチ」発表

藤堂様からもコメントをお寄せ戴いておりますが、先日、アメリカで開催されたCES 2016で、かねてよりスマートウォッチ市場に参入を表明していたCASIOがAndroidスマートウォッチの販売を表明しました。
いわゆる「フラッグシップ」という扱いで、高めの価格。CASIOとしてはそれだけ自信を持っている製品なのでしょうが、実入りの少ない当方としては「チプカシ」(チープ・カシオ)の様な腕時計に目が行ってしまう。それは蛇足ではありますが、スマートフォンの普及する今日、腕時計の需要は減る傾向にある様子。しかし、AppleがApple Watchを投入した事で、高級腕時計の需要が喚起された模様。そこでCASIOも高級腕時計としてのスマートウォッチを投入する決断に至ったのは、素晴らしい事ではあります。

電卓の方面では昨年、電卓事業50周年記念モデルとして「S100」を発売しました。キーボードなどに工夫を込め、耐久性も高い高級電卓です。これなどは「事務用電卓のフラッグシップ」なのでしょう。高機能電卓は、fx-CG20やfx-FD10 pro辺りがフラッグシップに相当するのか ?
すると、次は高機能電卓の「チプカシ」を期待したい。チプカシは、見た目は安いながらも十分な機能を持ち合わせている「CASIOの普及価格腕時計」であり、高機能電卓でも、こんな感じの製品を期待したいのであります。

過日、CASIOの電卓投票では高機能電卓・fx-5800Pが堂々の3位となりました。高機能電卓もそれだけ需要があるのですから、後継機では、もっと面白い仕掛けを用意して欲しく思うのです。

2015年11月28日土曜日

祝 ! [e-Gadget] 開設2周年

CASIO fx-5800Pの情報を中心に、CASIO BASICの情報を積極的に発信されておられますやす様のblogが開設2周年との事で、蔭ながらお祝い申し上げます。

e-Gadget - プログラム関数電卓
http://egadget.blog.fc2.com/

やす様のblogでは、最近はCASIO BASIC　の構造の詳細を述べられており、その情報を元にSentaro 様が高速動作版 CASIO BASICインタプリタ「C.Basic」の開発を進められております。既にLifeGame がウレシイくらいの速度で動いているらしく、とても興味深い成果なので、fx-9860GIIのユーザーの方は、是非とも御覧頂き度。

以前のポストで「10月にもCASIOがスマートウォッチを発表する」と書きましたが、どうも来年の話になっていた様です。
また、CASIOは電卓50周年のキャンペーンもやっており、既に年表も出していますが、そこには50年記念電卓「S100」も出ているので、年内の新しい高機能電卓は期待出来そうもありません。もうすぐ年末ですから「夢」は来年に持ち越し、にしておきましょう。尤も、CLASSWIZ は昨年12月に発表があったので、ウレシイ「裏切り」があって欲しい気もします。

他方、経営難が続くシャープですが、こんなTweetを出していました。

https://twitter.com/SHARP_JP/status/659206364967780352

ぜ…全国の高校生よ高専生よそしてその親御さんのみなさま……このオラにほんのちょっとずつだけ弊社の電子辞書とか関数電卓を採用する圧力を…学校にかけてくれ…!!!

高校生向けの関数電卓だそうです。

http://www.sharp.co.jp/calc/lineup/education.html

惜しむらくはプロ電ではない所でありますが、今のSharpの企画力では、こんな所から始めるのが精一杯なのかも知れません。しかし、発売は来年2月となっております。それまで持つのか !?
いや、それこそ「来年の事を言うと鬼がナントカ」とか申しますから、CASIOの新製品と共に、静観します。

さて、CASIO機ばかりではマズいので、HP機の話題も。JulyさんのAmazon売り場のHP50G, HP Primeですが、最近になって猛烈にバイが進んでいる様子。特にHP Primeは残り5台。もうすぐクリスマスという事で、高機能電卓のお買い求めを。

「すり抜け」確率の研究

火砕流に遭遇する確率に関するTweet

https://twitter.com/HayakawaYukio/status/658584031123107841

1万年で火砕流90％は、人生100年とすれば、一生の間に遭遇する確率（すなわち火砕流で死ぬ確率）は0.9％だ。そんなに小さくない

「1万年で90%」という事から、100年では90%/100 = 0.9% という目安を出しておりますが、確率計算の上では、どうなのでしょうか。
早川センセは、発生(火砕流に遭遇)する確率をp=一定としている様ですから、ここでも、その前提で計算を進めようと思います。

「1万年で90%の確率」とは言いますが、そのうちの「人生100年」は、1万年の1/100期間であります。
そこで、この「人生100年」を1回の施行とすれば、「1万年」は100回分の施行と考えてよいでしょう。
ここまで単純化して、いよいよ本題です。
100回(=1万年)の施行で90%となる発生確率の場合、1回(=人生100年)の試行による確率がどの程度になるのか、という問題に還元される事となりましょうや。

これを考える為に、複数回試行の確率を考えます。まず，最初は2回の試行から始めましょう。

p=2/3 という確率で生じる事象があるとして、1回で生じる確率は2/3であり、生じない確率は1/3である事に議論の余地はないでしょう。それでは、2回ではどうでしょうか ?

最初の1度目で生じる確率は2/3, 次に生じる確率も2/3 でありますが、ここには、次の組み合わせがあります。

1度目	2度目	確率
生じる	生じる	2/3 * 2/3 = 4/9
生じる	生じない	2/3 * 1/3 = 2/9
生じない	生じる	1/3 * 2/3 = 2/9
生じない	生じない	1/3 * 1/3 = 1/9

2回の事象に於ける「組み合わせ」では、実に8/9の場合が「生じる」を含んでいる勘定になります。
すなわち、2回の施行において「全く生じない」場合は1/9であり、いつの時点で起こるかを特定しない「起こる」確率は、全体の「全く生じない」場合を除いた場合として考える事が出来、その確率は 1-p = 1-1/9 = 8/9 と計算できます。

さて、この部分をもう少し考えてみます。
2回の施行で「全く起こらない」確率は、1/3*1/3 = 1/9 でしたが、つぎの様に考えてもいいでしょう。

n回の試行で「全く起こらない」確率は、p = (1/3)^n
n回の試行で「1度以上起こりうる」確率は、q = 1-p = 1-(1/3)^n

これをより一般的な表現とするならば、1回の試行で発生する確率を P として、

n回の試行で「全く起こらない」確率は、p = (1-P)^n
n回の試行で「1回以上起こりうる」確率(pの余事象)は、q = 1-p = 1-(1-P)^n

という、少し込み入った表式を得ます。
ここでは複数回の施行で発生しなかった事を「n回の試行をすり抜けて来た」と考え、「すり抜けの確率」と呼称します。

この辺りの考えをまとめて、件の「火砕流に遭遇しない確率」を、今一度考えなおそうというコンタンです。
100回の施行で90%の遭遇確率とは、100回の間、遭遇せずに「すり抜けて来た」確率の余事象ですから、

0.9 = 1-(1-P)^100

となります。但し、Pは「1回で火砕流に遭遇する確率」です。関数電卓をチョイチョイと叩けば、Pの値はたちどころに得られます。

1-P = (1-0.9)^(1/100) = 0.977237221...

P = 1-0.977237221 = 0.022762779, およそ2.28% という値が得られました。
0.9%よりも大きい値で、センセの仰る通り「そんなに小さい値ではない」とは言えます。
一方、危険率を5%に採ってしまうと、「2.28%は滅多に無いんだから、無視しても良かんべ」という事になってしまう。危険率の設定というものを、少し考えてしまう話題です。

昔懐かしい、プロ電のプログラムの技巧

かなり昔の事ですが、プロ電にプログラムを仕込んで「ゲームマシン」にするという本があり、その中に面白い技巧が出ておりました。

1. 「角度スイッチ」をジョイスティック代わりにつかう
2. GRADモードの活用で、乱数ルーチンのステップ数を稼ぐ

まずは、1の方から

1. 「角度スイッチ」をジョイスティック代わりにつかう

今でも、大抵の関数電卓では「DEG」(度)、「RAD」(ラジアン)、「GRAD」(グラード、グラジアン)の3つの角度モードが設定可能ですが、昔のプロ電の中には、三角関数の角度モードをスライドスイッチで切り替えるものがあったそうです。残念ながら、今となっては何と言う機種かは判りません。
このプロ電用に作られたゲームの内容も今となっては不明ではありますが、スクロール(シューティング)ゲームだったのだと思われます。画面(恐らく1行表示)を見ながら、角度モードのスライドスイッチをカチカチと切り替え、自機の操作を行っていたらしい。
スライドスイッチをジョイスティック代わりにするという発想が驚きですが、これにより、キー入力による実行停止がなく、リアルタイム性が向上したのだと推察されます。

角度モードを切り替えるスライドスイッチの状態をどうやって検出するのでしょうか ?
答は簡単で、適当な三角関数値を計算させ、角度モードによって違う値が得られるので、分岐して処理する、という具合です。

2. GRADモードの活用で、乱数ルーチンのステップ数を稼ぐ

GRADモードでは、直角を100とした角度指標となっていて、普段、余り使われない様に思われるGRADモードですが、こんな利用法があった、と言います。
逆三角関数、arc-sin, arc-cos をGRADモードで計算すると、0から100までの数値が返ります。0〜99までの乱数値が必要な場合、内蔵の乱数(0〜1未満)を実行、その逆三角関数値を計算し、最後にround(数値丸め)を実行する事で、0〜99までの乱数値が得られるという手法でした。
通常、こうした場合には乱数値を100倍し整数部を採るという操作を行いますが、100倍するという操作で「×」、「1」、「0」、「0」と、4ステップ余計に喰ってしまいます。メモリの少なかったプロ電では有効だったのかも知れません。

この技法は、38ステップのプログラミングが出来る機種(CASIO fx-3800P)に向けたものでした。今から思うに、38ステップのプログラムでナニが出来たのか、大変考えさせるものがあります。
プログラムで使える命令についても、制御機構は簡単な条件判断(レジスタ・メモリの0判定)と「文頭へ戻る」だけです。それらは、電卓に用意されていた数値積分機能に寄与するもので、積分を計算したい数式を「プログラム」によって実装し、数値積分を行っていたのだと思われます。数値積分は、微分方程式の数値解を得るのにも使われておりました。
そこで、更に「これ使ってゲームも作れるのではないか」という発想で、GRADモードを利用する技巧が編み出されたに相違ありません。

興味深い技巧ではあります。しかし、逆三角関数を通してしまう事で、乱数値の分布に偏りが出たりしないのか、という疑問があります。グラフ電卓で分布を描いて、調べてみましょう。

大抵のプロ電には、一様乱数の機能があります。これをarc-sinに通して、標本を取得します。電卓の実験例なので、100個くらいのデータが丁度良さそうです。
得られる標本の数値範囲は角度モードによりますが、ここでは分布の傾向をみるだけなので、どの角度モードでも構いません。当方、普段はRadianにしているので、Radianならば0～PIの範囲となります。そして、この0～PIの範囲を10の区分にして、棒グラフを描いてみる事にしました。一例では、確かに偏った分布となっております。

更に、統計的にこの分布が偏っているのか、調べる事も出来ます。
推計学の基礎が教える所によりますと、この標本集団から得られるχ自乗値は、近似的にχ自乗分布に従う、との事です。
高機能電卓でも、機能が充実しているものならば、χ自乗分布のCDF(累積分布関数)値を計算できます。これを使って検定も出来るのですが、計算方法は「χ自乗検定」を調べれば割合簡単に判ります。この例の場合、一様乱数と近しいかどうかを検定する事になり、一様分布の「期待値」としては、どの区分も同じ値になっている、として期待値を設定すればよい筈です。

2015年10月9日金曜日

10月になりました

10月になりました。

Julyさんの所の電卓、全て完売になっており、再入荷の目処は立っておりません。ウーン。

CASIOが「新しいスマートウォッチを投入する」としたのが、今月です。さて、何が出るのやら ? ついでにオモロイ電卓を期待したい。

CASIOの電卓も50週年とかで、キャンペーンを実施中だそうです。

カシオ電卓50周年
http://casio.jp/dentaku/sp/50th/

新しい高機能電卓を期待したい所です。

さーて、何か考えないとなァ ... 。

2015年9月26日土曜日

Julyさんの所のHP Primeが「おいでおいで」をしている ...

先ずは、最近は入手しにくくなりつつあるHP電卓マニュアルなどのPDFの一部について、リンク(直リン含む)を載せておきます。

HP35S Training modules
http://h30499.www3.hp.com/t5/Calculators/HP35S-Training-Modules-for-specialties/td-p/5854845?notmigrated#.VfKfk7O-feQ

HP50G Training modules
http://www.ele.uri.edu/faculty/vetter/Other-stuff/HP-calculators/HP-50g/index.html

HP42S manual
http://h10032.www1.hp.com/ctg/Manual/c01094960.pdf

HP15C manual
http://h10032.www1.hp.com/ctg/Manual/c01095005.pdf

HP12C manual
http://h10032.www1.hp.com/ctg/Manual/c00787560.pdf

HP12C クイックスタートマニュアル
http://www1.hp.com/ctg/Manual/c01798125

少し前の事ですが、某県知事が「高校女子に三角関数の勉強など不要」と公言し、インターネットでは「何を言うか、女性蔑視だ」などという意見が多く出回りました。そうした意見を受け、訂正を述べたとの事。県知事自身が「俺も使った事ないしネ」と言ったとかナントカ。

「女子に三角関数必要ない」 - Reuter
https://jp.reuters.com/article/2015/08/28/idJP2015082801001432

https://twitter.com/Stakesh/status/637048757142458368

https://www.facebook.com/photo.php?fbid=1522492048040753

この「訂正」は、結局の所「男女問わず、高校生程度のガキにサインコサインの勉強なんて、いらねんじゃね」という事になっております。
高齢に達し県知事になった所で人生を振り返ってみれば、確かに三角関数の計算など無縁だったのでしょう。大学出ても、計算を必要としない人は少なくなく、確かに申す通りなのかも知れません。

日本で少なく売られている関数電卓について、土地家屋診断士試験対応電卓などを謳っているものもあり、三角関数が生活にまったく不要なのか、というとそうでもなさそうです。
遠山啓・著「数学入門 (下)」(岩波新書)には、洋服の縫製について、袖の展開を行うと接続面にsin curveが現れる、と書いてありました。工業高校では間違いなく必要であり、高専なんかでは入学してからすぐに微積分を教えるとか。
商業高校ではサインコサインの計算には余り縁がないのかも知れませんが、複利計算を考えると、指数・対数の計算が顔を出します。
「高校生程度のガキに三角関数は不要」とは限らないのだと思う所です。

過日、Twittreで、こんな話題がありました。(リンクを持っていないので、記憶に基づいて文意を書いております。御容赦)

「せんせー、なんで三角関数なんて勉強すんの」というガキには「それは、こうしてバカなガキに三角関数を教えるお仕事が存在するからだ」と言うとガキは納得したので、おススメです

ネタとしては、なかなか面白い答えですが、現実では、バカなガキが言うのではなく、県知事がこんな事を言うております。
こういう県知事が出てくる事こそ「この国の教育が荒廃した」という証左なのですが、こうしたイカレテいる知事の話は、これ以上相手にしていると気が触れそうなので止めます。

CASIOの社員ですらも「高機能関数電卓の高機能を使う必要を感じていない」とは以下の記事でも書かれております。

異様に欲しくなった電卓の話カシオCLASSWIZ - 週刊ASCII
http://weekly.ascii.jp/elem/000/000/339/339010/

関数電卓の基本的機能さえあれば十分、という意識なのでしょう。
ただ、こういう意識で関数電卓、高機能電卓を使っているとなると、確かに需要は減っていく。

所が、海外ではAggressiveに、関数電卓、高機能電卓が売られているのです。TIはトップシェアで、CASIO,HPと続いているらしいのですが、Sharpでさえも海外では高機能電卓を出している !

Sharpも海外ではEL-9950なんて製品を出しているそうで、当方としては「EL-9950という電卓がある」という事に驚いたのでした。これは、過去に発売されていた「EL-9900」の改良新機種らしく、機能的にはTI-83+に相当する模様。

[シャープ] SHARP EL-9950関数電卓ライセンスシールを提供グラフィックコンバータ（海外直送品）- Amazon.co.jp
https://www.amazon.co.jp/%E3%82%B7%E3%83%A3%E3%83%BC%E3%83%97-EL-9950-SHARP-EL-9950%E9%96%A2%E6%95%B0%E9%9B%BB%E5%8D%93%E3%83%A9%E3%82%A4%E3%82%BB%E3%83%B3%E3%82%B9%E3%82%B7%E3%83%BC%E3%83%AB%E3%82%92%E6%8F%90%E4%BE%9B%E3%82%B0%E3%83%A9%E3%83%95%E3%82%A3%E3%83%83%E3%82%AF%E3%82%B3%E3%83%B3%E3%83%90%E3%83%BC%E3%82%BF%EF%BC%88%E6%B5%B7%E5%A4%96%E7%9B%B4%E9%80%81%E5%93%81%EF%BC%89/dp/B00T7U4A4G/ref=pd_sim_sbs_229_8?ie=UTF8&refRID=0SSSXVYTXG92RNAAE4SE#productDetails

Sharpは、海外ではこんなのを出して居るのですね、何で日本では出して呉れないの ?

一方、以下の記事には興味深い事があります。世界的に見るとSharpの電卓シェアは、実の所かなり少ない。

4 電卓は今後どう進化していくのか? - 【レポート】「いつでも、どこでも、だれにでも」使える計算機を! シャープ、電卓発売から50年の節目 - mynavi
https://news.mynavi.jp/articles/2014/03/18/sharp_calculator/003.html

電卓は現在でも、全世界で1億3,000万台の需要がある。シャープは現在、年間700万台を出荷しており、そのうち国内が100万台、海外が600万台となる。

(国内における)電卓戦争の一方の覇者・Sharpでさえも、世界的には大した売上がなかった。そのSharpでさえも、海外では高機能電卓を販売している。結局の所、国内においては高機能電卓の需要がない、という事に尽きるのでしょう、残念ながら。

先日、Julyさんの在庫が復帰したHP Primeですが、入荷数が多く無かったのか、既に残り2台となっておりました。
他のHP電卓も大分さばいていて「もうHP電卓の販売は終了なのか ?」と思ったのですが、価格が少々高めのリアル店舗・東映無線さんでは在庫が潤沢な様子なので、Julyさんでも当面販売は続くのかなぁ ?

https://www.toeimusen.co.jp/~dp/

東映無線さんでは、HP Prime以外にも50Gなどの在庫が十分あるみたいですから、「どうしても50Gが欲しい」方も、今暫くは迷えそうです。(Julyさんから引いているので日本語クィックガイド付きの筈ですが、念の為、お問い合わせ戴き度)

当方はHP Primeの「日本語冊子付き」を待っているのですが、頻繁なファームウェアのアップデートと、それに伴うPDFマニュアルの更新もあって、紙マニュアルは印刷と製本が追いつかないのかも知れず、日本語冊子付きは当分出ない可能性があります。

2015年8月27日木曜日

今回は移植ネタ「酔っ払いの虫」

先日、親分からつぎの話題を紹介されたのですが、関連したものを色々と調べている内に、記事のアップが遅れてしまいました。

楽屋裏 - 酔っ払いの虫 (Drunk Bug) - e-Gadget - プログラム関数電卓
http://egadget.blog.fc2.com/blog-entry-445.html

以下では、この記事を「原作」と記しております。まずは、TI-83+への移植を試みました。

TI-83+版のコード

PROGRAM:DRUNK
:ClrDraw
:94/2→V:62/2→W
:While 1
:Pxl-Change(W,V)
:2*randInt(0,1)-1→S
:2*randInt(0,1)-1→T
:If V+S≧0 and V+S≦94 and W+T≧0 and W+T≦62
:Then
:V+S→V:W+T→W
:End
:End

randInt(Begin, End)は、Begin,Endの間の整数乱数を返す関数で、randInt(0,1)は、0もしくは1を返します。
2*randInt(0,1)-1 の部分で、-1, +1の整数乱数にしております。
原作ではMODを使っておりましたが、まったく同じでは芸がないので、少しひねってみました。

今回、Pxl-Change()を使って判ったのですが、不規則ながらもグラフィクスコマンドを使う所で、一部変数の値が書き換わってしまうケースがあります。CASIO BASICと同じ様な感じです。
マニュアルには「ただし、変数X,Y,T,θ,Tはグラフを描く際、更新されることがあるので、グラフ関係以外の値を代入することは避けたほうが賢明です」(ページ A-50)と明記されておりました。

原作の興味深い所としては、酔歩の「移動の等方性」を担保するため、斜め方向の移動のみを行っている所です。
斜め方向の移動は、縦、横の移動と較べると、√2倍の距離の移動になります。縦、横の移動と斜め方向の移動とでは、単位時間での移動速度が異なる事になりますが、縦、横の移動、もしくは斜め方向の移動だけならば、移動の等方性が確保されます。原作の場合、斜め方向の移動のみを行う事で、移動の等方性を確保し、更に視覚的にも面白い表現を実現しています。

直交座標系では、縦、横と斜め移動の等方性について考える事が必要になりますが、四角ではなく六角形を敷き詰めたマス目空間(Hexagonal field)では、黙っていても移動の等方性が確保されます。少し工夫をする事で、直交座標系のピクセル空間にHexagonal fieldを擬似的に表現できます。

PROGRAM:HXDRNK
:ClrDraw
:94/2→V:62/2→W
:While 1
:Pxl-Change(W,V)
:Pxl-Change(W,V+1)
:Pxl-Change(W+1,V)
:Pxl-Change(W+1,V+1)
:randInt(0,5)→D
:(iPart(D/2)-1)*2→T
:2(fPart(D/2)*2)-1→S
:If T=0
:Then
:S*2→S
:End
:If V+S≧0 and V+S≦94 and W+T≧0 and W+T≦61
:Then
:V+S→V:W+T→W
:End
:End

図が必要な事もあり、詳しい内容は、いずれ折を見て述べましょうか。

更に、HP Prime向けコードも。

// title : Diagonal random walk for HP Prime
// begin : 2015-08-02
// note  : 

EXPORT drunk()
BEGIN
// Clean the screen (G0)
RECT();

LOCAL i, p;
LOCAL x, y, v, w;

x := 320/2;
y := 240/2;

WHILE GETKEY() == -1 DO
  p := B→R(GETPIX_P(x, y));
  IF p<>0 THEN
    PIXON_P(x, y);
  ELSE
    PIXOFF_P(x, y);
  END;
  v := RANDINT*2-1;
  w := RANDINT*2-1;
  IF x+v>=0 AND x+v<320 AND y+w>=0 AND y+w<240 THEN
    x := x+v;
    y := y+w;
  END;

//  WAIT(0.1);

END;

//  key wait loop
REPEAT UNTIL GETKEY() == -1;
FREEZE;

END;

実行中に適当なキーを押すと、停止します。更にもう一度キーを押す事で、プログラムの終了となります。
エミュレータで動かすと「爆速」ですが、そんな時にはコメントになっている「WAIT(0.1)」を使ってみて下さい。

お次はHexagonal field版

// title : Hexagonal random walk for HP Prime
// begin : 2015-08-02
// note  : 

EXPORT hexdrunk()
BEGIN
// Clean the screen (G0)
RECT();

LOCAL i, p, d;
LOCAL x, y, v, w;

x := 320/2;
y := 240/2;

WHILE GETKEY() == -1 DO
  p := B→R(GETPIX_P(x, y));
  IF p<>0 THEN
    PIXON_P(x,   y);
    PIXON_P(x+1, y);
    PIXON_P(x,   y+1);
    PIXON_P(x+1, y+1);
  ELSE
    PIXOFF_P(x,   y);
    PIXOFF_P(x+1, y);
    PIXOFF_P(x,   y+1);
    PIXOFF_P(x+1, y+1);
  END;
  d := RANDINT(5);
  w := (IP(d/2)-1)*2;
  v := 2*(FP(d/2)*2)-1;
  If w==0 THEN
    v := 2*v;
  END;
  IF x+v>=0 AND x+v<320 AND y+w>=0 AND y+w<240 THEN
    x := x+v;
    y := y+w;
  END;

//  WAIT(0.1);

END;

//  key wait loop
REPEAT UNTIL GETKEY() == -1;
FREEZE;

END;

HP PrimeのPPL (Prime Programming Language)での注意ですが、「変数への代入には := を使うべし」です。
実は、ピクセル値の値を取得する所で「p = B→R(GETPIX_P(x, y));」とやって、かなりハマってしまいました。
これでは代入ではなく「比較」になってしまうのですが、構文上は問題がないのでそのまま実行出来てしまい、エラー発見に手こずる事になります。

小中学生の皆様はもう新学期だそうで、今更「夏休みの自由研究」でもないのですが、来年の自由研究のネタに使えるかも知れません (鬼が笑うとりますがな)。

登録: 投稿 (Atom)