au光はマンション入口までが光であり、そこから俺んちの部屋には既存の電話線、つまりメタリックケーブルに信号を載せておるおるおーーるず。そこで下り100Mbpsとか上り30Mbpsとか言うとる。数値はテキトー。
だが、G.fastという技術規格だと下りが600Mbpsぐらい可能となるのであります。というので、待っている。
いや、またないでauひかりは解約してフレッツ光にすればいいんだけどさ。これはフレッツ業者に確認したら俺んちの集合住宅(マンションとも言うらしいが恥ずかしいので書かない)では可能だ。配管があるので光ファイバーをそこに這わせるとか言う。管理人に聞いたら既に使っている事がいるって去年言うとったし。
つーわけで、auひかりを解約して、Nifty + KDDIの組み合わせを解約してかね、So-net + フレッツ光で新規契約してキャッシュバック狙いを計画しておるおるおーーるず（笑）
だが、KDDIのG.fastが今後の半年ぐらいで実現したら、それでもいいなあ（笑）
どゔたろ。
と思ってKDDIに聞いた。auひかりに聞いた。メールで聞いた。電話で聞いた。そしたら、わかりませんって（笑）
新聞発表では2018年の10月に、つまり去年の10月にVDSLはG.fast規格に切り替えていく、それで下り600Mbps超が狩野舞子だ、叶姉妹だ、可能となるのです、キリッ)
というような新聞発表は単なる宣伝に過ぎず実態は・・・・。くそ、どっかで調べるか。5chで、ちみぃらのauのVDSLはG.fastに切り替わりましたか、教えてケロって書くか。
というわけで、価格コムを見ると100Mbpsは数少なく、多くが1Gbpsになっている。NTT東西、ソネット等が活躍している。au, KDDIは遅れている印象だ。
まあしかし、いいや、どうでも。
2019年09月16日
VDSLの次のG.fastというの・・・期待しておるのだが
なんだ、タマネギ男って?
https://www.nikkansports.com/general/news/201909160000141.html
韓国のソウル中央地検は16日、文在寅（ムンジェイン）大統領の側近、曺国（チョグク）法相の親族に絡む疑惑を巡り、曺氏の親戚である30代の男について横領や証拠隠滅教唆などの容疑で裁判所に逮捕状を請求した。
タマネギ男って・・・・法相のことかいなっと。
法相の親族が不祥事かい。どうなることやらだが。大統領に対する批判が増えるであろう。まーどうなることやら、とっこっちゃないでござるが、タマネギと言えば、タマネギをスライスして生で食うと痛風に効くという民間療法があってな・・・・
試したけど効果はなかった（笑） 痛くなったらロキソプロフェンですね。テッカーー!!
一次元調和振動子のシュレディンガー方程式の数値解に対して知見が広がった てへ
櫻井捷海 パーソナルコンピュータを用いた量子力学入門 裳華房 1989年 を参考にしてプログラムをいじっていたのである。著者の時代はパソコンはNECの国民機とそれの互換機がありDOSが主流だった。
本のBASICのコードをCで書いたり、gFortranで書いたり、Pythonで書いたりしたが結果はいずれも右側で発散が生じた。
一次元調和振動子のシュレディンガー方程式は2階微分方程式であり単純化された次のようなものだ。
y’’ = (t^2 – E) y
これを4次R-K法で数値計算で解く。この方程式は固有値・固有関数を求める問題でもあり数値解の中から固有値・固有関数を満たすものを探す。固有値固有関数を満たす場合には境界条件が満たされる、すまわち遠方での波動関数は0になる。そこには粒子はないからだ。
そこで左端の初期値を設定して右側に計算していく。固有値・固有関数を満たす場合には右側で0に収束するのである・・・・はずである…・収束してぇ!!
そう、収束するのだが、さらに計算を続けると発散するのである。そのグラフは過去に何度か載せた。
私のプログラムは本のBASICを移植しただけである。gFortran, Cではそうだ。4次R-K法は単純な計算式だし。PythonではScipy というモジュールのodeint()を使った。これは4次R-K法の改良版だと思う。
そして、wxMaximaでも上の微分方程式の数値解を出して見た。wxMaximaには4次R-K法の関数が用意されている。そこでも結果は右側で発散した。
この数値計算で発散するのはどうしてか? うーーむ・・・・4次R-K法の限界か? うーーむ。
本では右側で収束すると直ぐに計算を打ち切っているのでグラフに出てこない。また、サイト記事で「物理のかぎしっぽ」にC言語で似たような計算しているのがあったが、これも右側で発散している。
微分方程式の数値解で右側に発散が見られる場合に、その発散は解ではなくアルゴリズムによるものなのか、それとも解であるのか? そこの区別をどうするか。ノータリンのワタシにはわからないである（笑）
というときに別の解法を見つけた次第である。
https://qiita.com/sci_Haru/items/edb475a6d2eb7e901905#%E5%95%8F%E9%A1%8C%E8%A8%AD%E5%AE%9A
一次元シュレディンガー方程式の数値解を求めている。手法は、ヌメロフ法と書いてある。初耳である。不勉強が身に染みるである(´・ω・｀)
結果のグラフの例である。Pythonである。numpy + matplotlibを使っている。Scipyは使っていないである。
このグラフを見る限り、横を20に増やしても発散しない。。。。気がする。
というわけで、ヌメロフ法とは? 上のサイトにリンクがあったのでクリックしたである。
https://qiita.com/sci_Haru/items/338a5bb68ce17189917b
この方程式の解法として，ヌメロフ法と呼ばれる非常に簡単かつ効率的な陽解法のアルゴリズがあ る。この方法は4次のルンゲ・クッタ法[1]よりも1次だけ精度が高い[2]。
なんて書いてあるである。すばらしい!! 1次だけ精度が高いというのであれば、4次R-K法を改良したバージョンもあるはずだなあ。それだと発散しないのかなあ。どうだろなあ。
ただ、2階微分方程式で1階尾分を含まない方程式で使えるというので・・・・なによ、それ。
では、別問題。Scipyにはodeint()の他にも積分器があるって話だが、ようするに各種のアルゴリズムを揃えているのか? もし各種のアルゴリズムがあるのなら試してみたいである。
自分で書くの・・・・やだーーー。。。。。コピペだけどさ（笑）
以上、数値計算の入門あたりでウロウロしたりオロオロしたりしているワタクシの紹介でした。
おお、Niftyが !!
Nifty + KDDI 光(VDSL)を使っている。5年めだ。5年前にキャッシュバック狙いでYahoo ADSLから乗り換えた。ADSLは遅いぞ!! 安かったけど。
https://setsuzoku.nifty.com/niftyhikari/cam/nojima.htm?utm_source=nifty&utm_medium=TOP&utm_campaign=top_priority
で、当時はNIftyは富士通系列だったと思うが今はノジマが親会社になっている。そこで家電のノジマとのコラボだ。相乗効果だ。8%割引だ。9月末までだ、どうだ。これこそノジマがNiftyを傘下に置いた理由だ。テキトーだ。
で、昔のインターネット草創期の有名なプロバイダーだったNiftyはノジマに、NEC系列のBiglobeは今はKDDI傘下になっているではないか。ドッ日ー!! ドッヒーー。
Niftyはパソコン通信の頃から知っている。昔の昔、300bpsのモデムの海賊版キットを作って8bit PCのSharp X1のBASICでポート叩いて電話かけた。ピーヒョロヒョロの音が耳に心地よく流れる文字も読めた（笑） 哀愁の、郷愁、懐かしい時代を思い出し涙がボロボロ・・・でもないな。てへ
https://www.nifty.co.jp/company/outline/
しかし、儲からないから富士通系が手放したと考えるとノジマに変わっても儲からない構造は変わらないのではないか。そこだな、問題は。まー今のところは正常に運営が続いているようだし。そのうち様子がおかしくなったら、とっとと別のプロバイダーにキャッシュバック狙いで乗り換えだなあ、うむうむ。
ところでNiftyの社長は野島という若い人だが、ノジマの社長は野島というジーサンみたいだけど親子なのか。ノジマって野島という人が作ったからノジマか。するとニトリは二鳥という人がつくったからニトリか。二鳥か似鳥か煮鶏か知らんが。
おや、PCデポの社長も野島という人だな。一族なのか。まーどうでもいいけど。
で、Biglobeだが。NECを飛び出してから今はKDDIの完全子会社である。うーーむ。
そこだな。KDDIはau netというプロバイダー事業をやっている。6年ぐらい前はそこ使っていたし。そして、KDDI傘下にジュピターテレコムもおる。Jcomだ。ケーブルテレビの。そこはプロバイダー事業もやっておるおるおーるず。
では、KDDIがau net, Big lobe, Jcomとプロバイダーを3つも持っていてメリットは何さ? 相乗効果というものは何かあるのか? どうよ、KDDI ぬぬ
なんちて。
2019年09月15日
一次元調和振動子のシュレディンガー方程式の数値解で四苦八苦して苦しー!!
y'' = (v(t) - E)*y これ。無次元化した後の微分方程式。2階微分方程式である。初期値を与え、4次R-K法で解く。数値計算する。v( t) = t^2
櫻井捷海 パーソルコンピュータを用いた量子力学入門 裳華房 1989年 を参照のこと
その本ではBASICでプログラムを書いておった。それを自分でc言語で、gFortranで、wxMaximaで、Pythonで・。・・・やってみた。Pythonでの結果の例を次に示す。 
右側で発散しているのは無視する・・・としていたのだが気になった。wxMaxima, c, gFortranでも同じように右側では発散する。計算は左から始める。
で、これは4次R-K法でも誤差の蓄積によって発散に至るのだろう・・・・と思っていたが実はそうではないのではないか、と思い始めている。いや、誤差の蓄積なんだろうけどさ。
この微分方程式は解析的に解がでているので、それだと発散しないのだろ・・・確認したわけではないが。
y'' = (v(t) - E)*y この微分方程式を y’’ = y という形に無理やり単純化すれば、解は簡単で、y = e^x だ。それ指数関数なのでxが増えると発散するよ。
なのでxがある領域を超えると、この微分方程式の解は発散する・・・・と思った。解析的な解が発散しないのであれば数値解では何が原因でそうなるのか解明したい。
固有値・固有関数としてはE = 5は解析的にも解であって数値解でもグラフは左と右で収束していて。発散の前の0担っているところで計算を止めるといいので（笑）
この微分方程式は初期値を x= 0で0に近い値に設定して数値解を出せば発散するんでね? と思った。ようするに4次R-K方だろうが修正したアルゴリズムだろうが発散するのだ。と思う。
x = 0で初期条件を0近くにして計算したらこうなった。
あるてころで、解がy = e^xの形に近づいているのでそうなるのではないか? 正負はともかくさ。
ここらの話は教科書てあれば右側で収束した領域で計算を止めているし、サイト記事もそういうのがヒットした・・・きがするぞ。
だが、左側から計算を始めて右側も収束した後も計算を続けると発散するのだ。wxMaxima, C, gFortranでも同じ結果。
発散まで計算した例は「物理のかぎしっぽ」に載っていた。Cで計算したとか。
というわけで、疑問は発散するのが本来の解ではないのか? いや、解析解が発散はしないのだから数値計算上の問題だろ、その発散は・・・・ とか。
で、右側で微分方程式がy'' = (v(t) - E)*y から y’’ = y という形に近似できる領域では発散して当然だと、さっきビール飲んでいて気づいたのワタクシ。
発散のように見えるのが本来の解であるという可能性はどうなのか?
ボテンシャルが対称な場合には波動関数も対称になる・・・として右側で一度収束したら計算は止めるって妥当か?
この微分方程式の数値解で発散しない計算結果って、どっかにある?
というわけで今の所、ワケワカメなのであります。賢人・先人の解説を探しても右側をずーと計算した例が「物理のかぎしっぽ」しか見つからなかった次第であります。
科学技術計算の達人に知り合いはいないしよ。
もっとも、ここらが理解できたとしても、それがどしたのって話だけどね。オレ、引きこもりでアル中 予備軍だし（笑）
分かりが遅いボクのために例題はたくさん・・という本やサイトを探す Python
クラスがまるで理解できない。分かりが遅いボクのためには例題がたくさん必要だ。たくさん見ているうちに、ある時点で、あ、そういうことか!! と気づくとその一瞬が理解の瞬間だ。
サイトにしても本にしてもクラスを説明するために短いコードを出す。初心者向けとしてはいいのだろう。だが、何かがちんぷんかんププんがブン。
要するに、なんでこんなことするの、別の方法でもできるじゃん、関数定義してやれば簡単じゃん、・・・・・とかね。
原理的なところの説明のために短いコードを用意する著者の努力はわかるが、そんな約立たない例題ではなくてさ。多少、長くてもいいから実用になりそうな例を持って来いやーー てへ。
クラスの説明をカラスがカーと鳴くとか、カラスは黒いとか、そんな説明のものを見ると、オレとは次元が違うわな、こいつらは・・・・と思う次第である。
というわけだが、本を探す。図書館で。中古はAmazonで送料込みで500円未満を探す。
考えればわかることだがオレがプログラムを書くとしてもクラスを使うほどのものは書かない。せいぜい、100行前後。長くても500行ぐらいだ。クラスは使う側であり作る側ではないのだ、オレは・・・・・Python + Scipy + matplotlib で科学技術計算の入門みたいなところでアーーソブってのが趣味だし（笑）
クラスを理解して操れるようになろうって妄想であり、高望みであった ぎゃっふーーん!!
2019年09月14日
2階常微分方程式の数値解の練習 Python + Scipy + matplotlib
一次元調和振動子のシュレディンガー方程式の数値解だが、元の微分方程式は無次元化して単なる2階常微分方程式である。これの数値解を計算するが、ここで固有値・固有関数を求めたいわけだ。
そこでScipyのodeint()を使うと簡単に計算が可能となるのであります、ペカーー!! 
解析解は1,3,5,7,9….とわかっているので確認。8.0近辺に解はないって確認した。ポテンシャルが対称なので、ここではx^2 という単純なポテンシャル。その場合には固有値・固有関数を満たす波動関数は原点を通る。あるいは波動関数の導関数が原点で0である。なんちて。
このグラフは原点を通らないし、右側で発散しているので境界条件も満たさない。よって解ではない。てへ。計算は左から始めている。
櫻井捷海(さくらい かつみ) パーソナルコンピュータを用いた量子力学入門 裳華房 1989年
30年前の本だぜ・・・ちみぃら生まれてないがなーーー
BASIC, Turbo Pascalでのプログラム例が載っている。それをPythonでやってみたという話ね。
4次R-K法のアルゴリズムは書かず。odeint()という用意された関数を使うでし!! Pythonが科学技術計算において意外と高速なのはScipy, Numpy等のモジュールがCやFortranで作られたものだからである・・・・とかどっかに書いてあった。Python自体の四則演算を使うのではない。
Python + Scipy + matplotlib が使えるとグラフ作成機能付き高級電卓の代わりになるよん（笑）
肝心の量子力学についてはシュレディンガー方程式が2階線形常微分方程式である、という事を知ったのでいいや、もう。上の本は入門なのでね。その次の本が出ているのかどうかは知らない。
ただ、この手の本は欲しいよね・・・・。C/C++でもJAVAでもPythonでもFortranでもいいけどさ、入門の次を紹介するような本が。
物理の法則とか現象とかの視覚化ってのが初学者には有効であるぞよ。Mathematicaを使った視覚化ってのはどっかの大学の先生のサイトにあったけど。
2019年09月13日
二刀流は終わりか・・・・しーーん
https://headlines.yahoo.co.jp/hl?a=20190913-00000135-spnannex-base
大谷は怪我が多いなあ。これでは二刀流は無理・・・・二刀流どころかプロ野球選手としての寿命そのものも怪しいと思う。打者か投手の一つに絞れば良いという話ではなくなって来るぞ。
というわけで今季は終わり。来季はどうなるか。回復次第だ。投手でも打者でもいいけど試合に出て活躍しておくれよ。と思う次第である。なむ
考えるとわからなくなるっての、面白いよね しーーん
中性子星ってものがある・・・・らしい。ブラックホールの話では必ず出てくる中性子星。これ、中性子がほとんどである・・・・と言う場合に何を想像するか。
まず原子核の内部の中性子は核子どおしが強い相互作用で結びついている。だから中性子星も強い相互作用で固まっているだろう。そして電気力が働かないからいくらでも接近して詰め込むことができる。重力によって詰め込まれ、収縮して中性子星になるのであろう。
強い相互作用 + 重力による収縮で固まる。
では、中性子星は壊れるのではないか。重力によって押しつぶされて中性子星は壊れる。つまり、中性子は壊れる。壊れてuクォークとdクォークがバラバラになる。グルーオンも開放される。
そうならないのはなんでか? 知るかよ。そこだな。
太陽のような水素ガスの星は重力で壊れないように中心部で核融合が起きて収縮に耐えている。この核融合は見方を変えると物質の製造工場の役割を果たす。中心部に鉄ができると核融合は止まる。すると重力による収縮に耐えられなくなる。どっひゃーと鉄が壊れて膨大なエネルギーが放出される。超新星爆発である。
その後、条件に依って中性子星ができる。その中性子星は壊れることはないのか。
原子は重力による収縮のエネルギーで壊れるが陽子や中性子は耐えるのか。
陽子や中性子も壊れてクォークがバラバラと放出される・・・・とかならんのか。
陽子や中性子を壊すにはどれだけのエネルギーがあればいいのだ? 知るかよ。
というような妄想をすると楽し～であるよ。
ただ、ビッグバンの話があると初めはクォークがバラバラでクォークのスープ状態であると書いてある。つまりクォークが単独で存在するのだ。ちょ～エネルギーがチョーチョーデカイ場合には。その後、膨張とともにエネルギーが下がりクォークが複合粒子を作り始める。
では、宇宙にクォークがバラバラになるぐらいの、ちょ～チョーエネルギーが高い状態はどこにもないのかいなっと。知るかよ。あったらどうなるってんだよ。知らんがな。
とにかく、宇宙論の入門書は妄想すると楽Ｃーですわ。どうせ理屈はわからんし（笑）
GmailでNiftyを設定する場合に、どんづまり・どんづまり・くそったれ
俺んちの場合、このような設定ではだめだった。ユーザー名とパスワードが違いますと出て進まず。
上の場合、ユーザー名は@nifty.comを削除。そうしたら俺んちでは設定が出来た。
Android 7。
半日かかったぞ、バカタレの鎌足!!
