なぜPerlで書くのか？

自費出版掲載のスクリプトは主にPerlで書く予定です
で、僕は何故Perlで書くのか？

・Perl5.xはほぼ言語仕様が確定している
Rubyを勉強してRubyで書いてもいいんですけど
Rubyはver upするとスクリプトのメンテナンスが必要な気がします
WEB上だとメンテナンス出来ますが、書籍上だとできませんよね
そういう意味で自費出版用のプログラミング言語としてPerlの良さがあります

・スクリプト言語である
なんでスクリプト言語にこだわっているんやろ？
HSPもスクリプト言語ですし
（ちなみに僕は主にHSPやってます）
でもスクリプト言語って手軽に書いて手軽に実行ってイメージありますよね
手軽に書きたいとは常々思ってますけど
Cも言語仕様がほぼ確定してますが、スクリプト言語じゃないから採用しないのかな？
ちなみに今は何回目かのJavaの勉強中だったりします
JavaとPerlを比べたらやっぱりPerlの方が手軽ですよね

・割とメジャーであり、研究用っぽいイメージがある
僕は主にHSPやってるんですけど、HSPってマイナーですよね
ってか研究用のスクリプトをHSPで書く人はいるのでしょうか？
研究用としてはCやLISPがメジャーなのかな？イメージですが
Cより研究用として向いていそうであり
LISPよりメジャーではあります
LISPは読みやすさ。より書きやすさに重点を置いているかんじしますよね

・色々と便利な機能がそろっている
Perlと言えば正規表現
僕は正規表現を書いたことは一度もありませんが
でも正規表現を除いても便利な機能があります
主にハッシュとリスト
CはハッシュとリストがないからPerlを採用した面も大きいです
ハッシュは使うかどうかわかりませんが
リストは使いたいです
リスト（実際は配列をこう使う予定です）に確率をぶちこんで
平均確率とか平均の点数とかを求めたいです
HSPにはハッシュとリストがないのがメインにしたくない理由です

・文字列操作に強い
CとPerlを比べて、ここでの一番のPerlの魅力は文字列操作に強いことです
Cの文字列は連結関数だとか比較関数だとか直勘的でないです
Perlは文字列操作に強いですよね
数値演算が主ではありますが、直観的に文字列操作ができるのはPerlの魅力です

・直観的に書ける
僕のプログラミング・スタイルにおいて直勘的に書けることは重要です
HSPでも直勘的な書き方がしたくて、よくラッピングします
色々な書き方が出来るのはPerlの特徴であり
その中で出来るだけ僕の直勘に近い書き方をして行こうと思っています

・第二言語としてHSP
自費出版用の文章に載せるプログラムは主にPerlで書く。と決めましたが
それでも第二言語としてHSPを使用する予定です
用途はグラフ描画です
HSPでグラフを描画するスクリプトを書いて自分でグラフを描けば
自費出版の費用も安くなるかなって
更にX軸Y軸Z軸の三次元的なグラフを3Dで描いて見やすく回転しよう。だなんて考えています
簡単に3Dの画像を書けるのはHSPの特徴の一つですよね

桁の小さい平方根の計算

とりあえず桁の小さい平方根の計算の処理は完成しました
次はBIGFloatを利用してもう少し大きい桁の平方根を計算したいです
完成まで長かったｗて、やってないから長いんですけどｗ

後、HSPでグラフを描くモジュールも作成したいです
Perlでグラフを描くのは僕の技量では無理ｗ

use strict; use warnings; my ($r,$ans); ($r,$ans)=&small_sqrt(20000.0,0); print $r." ".$ans."\n"; ($r,$ans)=&small_sqrt(0.0002,0); print $r." ".$ans."\n"; ($r,$ans)=&small_sqrt(20000.0,1); print $r." ".$ans."\n"; ($r,$ans)=&small_sqrt(0.0002,1); print $r." ".$ans."\n"; # 平方根の計算。計算する桁が小さい sub small_sqrt{ # 引数A：sqrtしたい値 my $value=$_[0]; # 引数B：これが1だとprint debug my $on_debug=$_[1]; # 返り値用 # $r:error情報。-1：数値が大きすぎる：-2数値が小さすぎる # $ans：解 my ($r,$ans); if ($value>=1.0){ # sqrtしたい値が1.0より大きい場合 ($r,$ans)=&small_sqrt_L_in($value,$on_debug); }else{ # sqrtしたい値が1.0より小さい場合 ($r,$ans)=&small_sqrt_S_in($value,$on_debug); } return ($r,$ans); } # 平方根の計算。計算したい値は1.0より大きい sub small_sqrt_L_in{ # 計算したい値 my $p=$_[0]; # これが1だとdebug my $on_debug=0; $on_debug=$_[1]; if ($on_debug==1){ printf "sqrt(".$p.")"."\n"; } # &some_divはsqrtしたい値が1以上の時呼び出す # sqrtしたい値は10の何乗の２乗か？求める my ($r,$r_A,$r_B)=&some_div($p); # &in_calc_0p1to0p0000001は実際に計算する my $ans=&in_calc_0p1to0p0000001($r_A); if ($on_debug==1){ printf $r."\n"; printf $r_A."\n"; printf $r_B."\n"; printf $ans*$r_B."\n"; } return ($r,$ans*$r_B); } # 平方根の計算。計算したい値は1.0より小さい sub small_sqrt_S_in{ # 計算したい値 my $p2=$_[0]; # これが1だとdebug my $on_debug=0; $on_debug=$_[1]; if ($on_debug==1){ printf "sqrt(".$p2.")"."\n"; } # &some_mulはsqrtしたい値が1より小さい時呼び出す # sqrtしたい値は0.1の何乗の２乗か？求める my ($r2,$r_A2,$r_B2)= &some_mul($p2); # &in_calc_0p1to0p0000001は実際に計算する my $ans2=&in_calc_0p1to0p0000001($r_A2); if ($on_debug==1){ printf $r2."\n"; printf $r_A2."\n"; printf $r_B2."\n"; printf $ans2*$r_B2."\n"; } return ($r2,$ans2*$r_B2); } # sqrtしたい値が1以上の時呼び出す # sqrtしたい値は10の何乗の２乗か？求める sub some_div{ # $r:error通知用。普通は1。10の301以上だとerror # $r_A:繰り返しで/100.0づつしていく # $r_b:10のX乗の２乗か？のXの値を入れる my ($r,$r_A,$r_B)=(1,$_[0],1.0); # $r_A=&mul_10_X_count(301); if ($r_A>&mul_10_300_count()){ $r=-1; } while ($r_A >= 100.0) { $r_A/=100.0; $r_B*=10.0; } return ($r,$r_A,$r_B); } # sqrtしたい値が1より小さい時呼び出す # sqrtしたい値は0.1の何乗の２乗か？求める sub some_mul{ # $r:error通知用。普通は1。10の301以上だとerror # $r_A:繰り返しで*100.0づつしていく # $r_b:0.1のX乗の２乗か？のXの値を入れる my ($r,$r_A,$r_B)=(1,$_[0],1.0); # $r_A=&mul_0p1_X_count(285); if ($r_A<&mul_0p1_285_count()){ $r=-2; } while ($r_A <= 1.0) { $r_A*=100.0; $r_B/=10.0; } return ($r,$r_A,$r_B); } # 計算の処理。in_calcを８回呼び出す sub in_calc_0p1to0p0000001{ # $p:この数値以下で最も大きい数値を導きたい。等しくてもいい # $ret0～$ret8:計算結果が入る。少しづつ小さな単位を計算して代入 # &in_calcの３つめの引数:計算する単位 my $p=$_[0]; my $ret01=&in_calc($p,0,1); my $ret02=&in_calc($p,$ret01,0.1); my $ret03=&in_calc($p,$ret02,0.01); my $ret04=&in_calc($p,$ret03,0.001); my $ret05=&in_calc($p,$ret04,0.0001); my $ret06=&in_calc($p,$ret05,0.00001); my $ret07=&in_calc($p,$ret06,0.000001); my $ret08=&in_calc($p,$ret07,0.0000001); return $ret08; } # 計算の処理。一桁だけ計算 sub in_calc{ # この数値以下で最も大きい数値を導きたい。等しくてもいい my $expression=$_[0]; # 計算結果を入れる my $exp01=$_[1]; # 計算の単位 my $unit01=$_[2]; # 返り値用 my $ans01=0; # 乗算の基 my @for_mul; for (0..9){ $for_mul[$_]=$exp01+($_*$unit01); } # 乗算の基を2乗 my @of_mul; for (0..9){ $of_mul[$_]=$for_mul[$_]*$for_mul[$_]; # print "i:$for_mul[$_]:i*i:$of_mul[$_]\n"; } # X以下で大きい数値を求める for (0..9){ if ($of_mul[$_]<=$expression){ $ans01=$for_mul[$_]; } } # print "解の候補:$ans01\n"; $ans01; } #printf &mul_10_X_count(0)."\n"; #printf &mul_10_X_count(3)."\n"; #printf &mul_10_300_count()."\n"; #printf &mul_10_300_count()*&mul_0p1_300_count()."\n"; #printf &mul_0p1_X_count(0)."\n"; #printf &mul_0p1_X_count(3)."\n"; #printf &mul_0p1_300_count()."\n"; #printf &mul_0p1_300_count()*&mul_10_300_count()."\n"; # 10の引数Ａ乗の数を返す sub mul_10_X_count{ my $x=1.0; my $i=$_[0]; $i--; for(0..$i){ $x*=10; } return $x; } # 10の300乗の数を返す sub mul_10_300_count{ return &mul_10_X_count(300); } # 0.1の引数Ａ乗の数を返す sub mul_0p1_X_count{ my $x=1.0; my $i=$_[0]; $i--; for(0..$i){ $x/=10; } return $x; } # 0.1の285乗の数を返す sub mul_0p1_285_count{ return &mul_0p1_X_count(285); }

裏テーマ

僕が今書いている自費出版用の書籍の裏テーマは
問題をプログラミングで解く。と
数式の解の近似値を得る。です

人には様々な問題がある
社会にも様々な問題がある
それをプログラミングを利用して解こう
そういうテーマです
今書いているやつには
人の問題、社会の問題を解くようなプログラミング例は
思いついていないので書けませんが
プログラミングで解けることは
数学的な問題や
物理的な問題だけじゃないと思うんです
それを書籍を書くことを通して
道筋をつけたい
そういうテーマです
なので今やってるルートの計算なんて
（それが終われば２乗である数になる近似値だけじゃなく
　Ｘ乗（まずはＸは整数）すればある数になる近似値を解くプログラムを書きたいです）
ググれば出てきそうですが
あえて自力で解こうと思っています

数式の近似値を得る。というテーマは
数式の解を解くって難しいんですよ
３次式の解を求める解の公式とか
やたらややこしいんですよ
でもそれは近似値ならば少数として
より把握しやすい値になるかもしれない
正確な解ってややこしいし
解き方も難しいし
なのでコンピューターの力で
強引に近似値を得てしまおう
書籍に載っているやり方を応用して
式の解を強引に得て欲しい
そういう道の道筋をつけたい
そういうテーマです

とりあえず

とりあえず、ルートを求める処理は出来てます
少々の改良と
（改良ってかスクリプトのサブルーチン化で呼び出せるように。です）
コメント書くのとサブルーチンの名前のリファクタリングが
TODOとして残ってるんですが
中途半端でも載せれば、眺めながらやる気出るかもしれないので
載せときますｗ

追記：一応、スクリプト完成したので削除ｗ

自費出版を目指すブログ

TACKは数学っぽい経済学っぽいプログラミングっぽい書籍を自費出版するのが夢です
数学と経済学とプログラムを足して３で割ったような書籍です
あんまり言ってないんですが（笑）僕の夢は作家になることです
作家は何冊も書かなければ作家とは言えなく
文章を書くことは好きですが、作家になれるほど大量に書いている訳でもなく（笑）

インタプリタ/コンパイラの作成を目指すブログ　で
このネタ書籍で使えるかも？と思ったネタを書いたので
それを機にGOOGLEのブログに投稿しながら執筆していこうかな？
と、考えるに至りました

はてなの方が見てくれそうですが
GOOGLEのブログはアフィリエイト？が使えるみたいなんで
収入になるかもしれないし
こっちで書くことにしました

幾つか寝かせてあるネタもあるので
ブログに投稿しながら
それをマイルストーンとして
ある程度の分量になるまで執筆して
そしたら自費出版しようと思いました

大学でゼミの頃、空想していた概念をまとめて自費出版する！というのは作家を目指す上での最初の僕の夢です！
最初は卒論にしようと思いましたがまだ構想があやふやで、実行するのに踏み切れませんでした
あれから１０年経つんだなぁ
大分形になって来ましたが
手を伸ばせば、遠ざかり
手を離そうとすると、近づいて来るようで
難しいですね
そろそろ概念をきちんと文章にすることを決意しようかな？と考えているところです
なんか理想の人みたいな概念みたい
手を伸ばせば、遠ざかり
手を離そうとすると、近づいて来るようで

自費出版したら絶版になるまで
書いた記事は削除する予定です

文章をよりよくするためのコメントもお待ちしています

BIGFloat

先日ルートを計算するPerlスクリプトを書いていて
Perlのdoubleの有効桁数は何桁だったかな？と思ってググったら
PerlにはBIGFloatがあるみたいですｗ
ルートの計算の後で桁の大きい少数を扱うスクリプトを書こうと思ってましたｗ

あろうがなかろうが自力で処理を書いてもいいんですけど
めんどくさいしｗバグが怖いので
素直にBIGFloat使う予定です

ルートの計算は大筋完成していて
後は整理してコメント書いて完成です

＃最近、コメント書くのめんどくさいｗ

上限と下限の条件分岐

rootの計算はdoubleを使ってるので
ものすごい大きい数から
ものすごい小さい数まで計算できるので
いらないかもしれませんが
一応、考え方として
どれくらい大きい数から
どれくらい小さい数まで計算できるか
決めておいた方がいいですよね

とりあえず、最大1e+300
最小1e-300まで計算できるように決めました
決めるったら決めるｗ

でそんなのを求めるスクリプトです

追記：doubleの下限は300桁ちょいくらいなので
doubleの桁数も考え最小1e-285に変更しました

use strict; use warnings; printf &mul_10_X_count(0)."\n"; printf &mul_10_X_count(3)."\n"; printf &mul_10_300_count()."\n"; printf &mul_10_300_count()*&mul_0p1_X_count(300)."\n"; printf &mul_0p1_X_count(0)."\n"; printf &mul_0p1_X_count(3)."\n"; printf &mul_0p1_285_count()."\n"; printf &mul_0p1_285_count()*&mul_10_X_count(285)."\n"; # 10の引数Ａ乗の数を返す sub mul_10_X_count{ my $x=1.0; my $i=$_[0]; $i--; for(0..$i){ $x*=10; } return $x; } # 10の300乗の数を返す sub mul_10_300_count{ return &mul_10_X_count(300); } # 0.1の引数Ａ乗の数を返す sub mul_0p1_X_count{ my $x=1.0; my $i=$_[0]; $i--; for(0..$i){ $x/=10; } return $x; } # 0.1の285乗の数を返す sub mul_0p1_285_count{ return &mul_0p1_X_count(285); }