N88-BASICで電離度 (3回目)

2026/2/28(土)
 
N88-BASICで電離度 (3回目)
 
(degree of ionization)
■ pHと電離度(3価の酸)
▼ 式
https://ulprojectmail.blogspot.com/2026/02/ph-3.html
水素イオン指数 (3回目)
 
より
 
▼ 水素イオン濃度
[H+]5 + Ka1[H+]4 + {Ka1(Ka2 - c) - Kw}[H+]3 
+ Ka1{Ka2(Ka3 - 2c) - Kw}[H+]2 - Ka1Ka2(3cKa3 + Kw)[H+]
- Ka1Ka2Ka3Kw = 0
 
[H2A-] = [H+]2[A3-]/(Ka2Ka3)
[HA2-] = [H+][A3-]/Ka3
[A3-] = cKa/{Ka+{Ka1Ka2+(Ka1+[H+])*[H+]}*[H+]}
α1 = ([H2A-]+[HA2-]+[A3-])/c
α2 = ([HA2-]+[A3-])/([H2A-]+[HA2-]+[A3-])
α3 = [A3-]/([HA2-]+[A3-])
 
▼ 酸の濃度
A = [H+]{[H+]{[H+]{[H+]{[H+]+Ka1}+(Ka1Ka2-Kw)}
+ Ka1(Ka2Ka3-Kw)}-Ka1Ka2Kw}-Ka1Ka2Ka3Kw
B = = Ka1[H+]{[H+]([H+]+2Ka2)+3Ka2Ka3}
c = A/B
 
▼ グラフ
横軸
濃度の指数pc = -log10 c (7が真ん中)
縦軸
水素イオン指数pH = -log10 [H+] (7が真ん中の黄色のグラフ)
電離度α1 (1目盛りが1の赤色のグラフ)
電離度α2 (1目盛りが1の紫色のグラフ)
電離度α3 (1目盛りが1の緑色のグラフ)
で描画しています
 
 
VL,NL,XL-BASICとdlg～.zip(alp003.bas)は
このブログ(以下のリンク)からダウンロードできます
https://ulprojectmail.blogspot.com
Readme.txtを読んで遊んで下さい















このブログの人気の投稿

NEWS

続きを読む をクリックでリンク表示 (かな入力,ビッグエンディアン推奨派)    n回目の次回はタイトルで検索するか目次で見つけて下さい    2026/3/28(土)  N88-BASICで電離度 (3回目)  公開しました    2026/2/17(火)  メモ (18回目)    2026/3/1(日)  VL,NL,XL-BASIC ver～28x5e  ダウンロード    2026/2/27 (金)  水素イオン指数 (3回目)    2026/2/21(土)  解析力学 (邪道編)    2025/12/17(水)  午前午後問題 (1回目)    2026/3/12(木)  掛け算の順序問題(新版) (3回目)    2025/8/25(月 )  加速度座標系について
続きを読む

N88-BASICでゲーム (1回目)

イメージ
2021/8/16(月) N88-BASICでゲーム  (1回目) by ULproject   球同士の衝突と跳返り   大文字ベクトル、小文字スカラー 文字の後の、 'は衝突後、数字は球番号   球同士の反発係数、衝突面に水平 ,垂直u,e 球の中心位置 P,半径r,質量m,速度V 衝突面の球 1方向の法線N = (P 1 -P 2 )/|P 1 -P 2 |   球同士の衝突は、 |P 1 -P 2 | ≦ r 1 +r 2   なので √{(P 1 -P 2 )･(P 1 -P 2 )} ≦ r 1 +r 2   (P 1 -P 2 )･(P 1 -P 2 ) ≦ (r 1 +r 2 )(r 1 +r 2 ) で判断できます。   ここで、ベクトルの計算ですが、 P(i) (i=0,1,2) をPのx,y,z成分と すると、 P = P 1  - P 2 は、 P(i) = P1(i) - P2(i) (i=0,1,2) P･Pは ΣP(i)P(i) (i=0,1,2) となります   BASICでは、 P(0) = P1(0) - P2(0) P(1) = P1(1) - P2(1) P(2) = P1(2) - P2(2) D = P(0)*P(0)+P(1)*P(1)+P(2)*P(2) R = R1 + R2 IF D <= R*R THEN 重なり又は接触 となります。   反発させるときは、球が接触している 間は再び反発しないようにしないと 球が離れなくなります。   次に反発についてですが、衝突面に 水平な成分は、球の回転に関係します が、ここでは球の回転は考えずに 垂直成分と同じように反発係数で 減速させる事にします。   衝突面に垂直な成分 (v 1 は球 1の速度、v 1 'は衝突後) 運動量 p = mv 運動量保存則 p 1 +p 2 =p 1 '+ p 2 'より p 1 +p 2  = m 1 v 1 '+ m 2 v 2 ' 反発係数 e = -(v 1 '- v 2 ')/(v 1  - v 2 )より ev 1  - ev 2  = -v 1 '+ v ...
続きを読む

d88ファイルの変換 (1回目)

2022/8/27(土) d88ファイルの変換 (1回目)   d88をbas(ascii)ファイル(PC-88)に変換する方法   ファイル変換は ～.d88 → ～.bas　　　　　　　  … ① ～.bas(binary) → ～.bas(ascii)　… ② の手順で解説します (Windows前提)     ① N88-BASICディスクイメージファイル～.d88ファイル から個別の～.basファイルを取り出す方法   Sasaji様のページ http://s-sasaji.ddo.jp/bml3mk5/#l3diskex から L3ディスクエクスプローラ(L3 Disk Explorer) を探して ダウンロード(Downloads) - (V0.6.0) [2022-05-05更新] をクリックしたページから   Windows 32ビット 実行ファイル (32bit Binary) l3diskex_win32_060.zip (2.2MB, 2022-05-05) をクリックしてダウンロードし 展開します   l3diskex.exe を走らせて、～.d88ファイルを読み込ませると 中身を見ることが出来、エクスポートで 個別のファイル(～.basなど)を取り出す事が 出来ます     ② バイナリーファイル～.bas(binary)を アスキーファイル～.bas(ascii)に変換する方法   TSAPPS開発室 様 のページ https://letsplay-retropc.tsapps.net/pc8801/n88basic-source-editor から ダウンロードとインストールをクリックした ページから   ダウンロードはここをクリック！ をクリックして n88srced_v103.zip をダウンロードして展開し n88basic_srceditor.exe をクリックし アスキーファイルでは無いを選んで 読み込み SAVEするとアスキーセーブされます   グラフィックキャラクターを含まない ファイルはメモ帳で編集可能ですが   グラフィックキャラクターを含む 場合は n88basic_s...
続きを読む