科学,数学,PCで可視化 © ULproject

2025/12/7(日)   等速円運動 ( 3 回目 )   ( uniform circular motion )   ■ 等速円運動 ▼ ベクトル 3重積 A×(B×C) = A×(ByCz - BzCy, BzCx-BxCz, BxCy-ByCx) = (Ay(BxCy-ByCx)   -   Az(BzCx-BxCz),    (Az(ByCz-BzCy)   -   Ax(BxCy-ByCx),    (Ax(BzCx-BxCz)   -   Ay(ByCz-BzCy)) =   (AyBxCy-AyByCx-AzBzCx+AzBxCz   + ( AxBxCx - AxBxCx ) ,   … (=0を足す)     AzByCz-AzBzCy-AxBxCy+AxByCx   + ( AyByCy - AyByCy ) ,   … (=0を足す)       AxBzCx-AxBxCz-AyByCz+AyBzCy   + ( AzBzCz - AzBzCz )) … (=0を足す) =   (Bx(AxCx+AyCy+AzCz)   - Cx(AxBx+AyBy+AzBz),      By(AzCz+AxCx+AyCy)   - Cy(AzBz+AxBx+AyBy),      Bz(AxCx+AyCy+AzCz) - Cz(AxBx+AyBy+AzCz)) = (A･C)B-(A･B)C A×(B×C) = (A･C)B-(A･B)C   ▼ 角速度 ω   ω   = (0, 0, ω)   R = (x, y, z) = (rcosωt, rsinωt, 0) V = ω ×R = (ω y z- ω z y, ω z x- ω x z, ω x y- ω y x) = (-ω z y, ω z x , 0 ) = (-ω y, ω x , 0 )   A = ...

202 5 / 12/2(火 ) N88-BASICでnのn乗の和 (2回目)   ( summation )   ■   nのn乗の和の1の位 ▼   問題 Σ n=1 10 n 100   = 1 100   + 2 100   + 3 100   + … + 10 100   の 1の位はいくつか   ▼   参照   nのn乗の和の1の位 (2回目)   より Σ n=1 10 n 100   = 1 100   + 2 100   + 3 100   + … + 10 100   ≡ 3 (mod 10) (1の位は3)     ■   多桁演算 (正の整数の和積) ▼   設計 Σ n=1 10 n 100   = 1 100   + 2 100   + 3 100   + … + 10 100   < 10 100 ･ 10 = 10 101   (102桁) より 文字変数 (最大255文字)で1文字1桁に対応させて足りるので 正の整数 255桁(文字)同士の和積を実装し計算する   n = 1～10までの各結果の下2桁を表示し 最後に Σ n=1 10 n 100   の結果と計算時間を表示しています   ▼   結果 Σ n=1 10 n 100   = 1 100   + 2 100   + 3 100   + … + 10 100   = 10000265616025915085580987810979861492643890241086 46749449153743033415896167749949249632000590464813 3 … (101桁)で 1の位は3という結果になりました     V L, NL ,XL -BASICとdlg～.zip( sum002 .bas)は このブログ (以下のリンク)からダウンロードできます https:/...

202 5 / 11/27 (木 ) nのn乗の和の1の位 (2回目)   ( summation )   ■   nのn乗の和の1の位 ▼   問題 Σ n=1 9 n 100   = 1 100   + 2 100   + 3 100   + … + 10 100   の 1の位はいくつか   ■   考察 ▼   m n (1≦n≦100)を10で割った余りを調べる           5,6,7,8,…,100 (1～10のn乗を10で割った余りを以下に示す)      n   =1,2,3,4,   1 n   ≡ 1       (mod 10)   2 n   ≡ 2,4,8,6 (mod 10) … (2 4 8 16 (3)2 (6)4 …)   3 n   ≡ 3,9,7,1 (mod 10) … (3 9 27 81 (24)3 …)   4 n   ≡ 4,6     (mod 10) … (4 6 24 …)   5 n   ≡ 5       (mod 10) … (5 25 …)   6 n   ≡ 6       (mod 10) … (6 36 …)   7 n   ≡ 7,9,3,1 (mod 10) … (7 9 3 1 7 …)   8 n   ≡ 8,4,2,6 (mod 10) … (8 4 2 6 8 …)   9 n   ≡ 9,1     (mod 10) … (9 1 9 …) 10 n   ≡ 0    ...

2025/11/21(金)   等速円運動 ( 2 回目 )   ( uniform circular motion )   ■ 等速円運動 ▼ 等速円運動 の軌跡 φ φ(x, y, z) = x 2   + y 2   - r 2   と置くと φ(x, y, z) = 0 は原点を中心とする半径 rのxy平面上の円を表す   ▼ 勾配 grad gradφ = ∇ φ  = (∂/∂x,∂/∂y,∂/∂z) φ = (∂/∂x,∂/∂y,∂/∂z) (x 2   + y 2   - r 2 ) = ( ( ∂/∂x )(x 2   + y 2   - r 2 ) , ( ∂/∂y )(x 2   + y 2   - r 2 ) , ( ∂/∂z )(x 2   + y 2   - r 2 ) ) = ( 2x ,   2y ,   0)   これは点 (x, y)での法線の方向を示すベクトルになっている様です   ■ 結果 ▼ 軌跡 φ φ(x, y, z) = x 2   + y 2   - r 2   = 0   ▼ 勾配 grad gradφ = ( 2x ,   2y ,   0)  … 点 (x, y)での法線方向  

2025/ 11 / 17 ( 月 ) NHKの集金でテレビを持っていないのでと帰ってもらった 翌週来たので理由を尋ねるとテレビを買ったかもしれない とのこと 買う予定はないと告げるとまた来ますと言って帰った それから引っ越すまでの3年間毎週2回ずつピンポンされた この嫌がらせが合法なの納得できない 週に2回も30分以上部屋を見せろなど言われ 見せてもなかなか帰らない 迷惑でしかない 出ないと毎回数分間ピンポンを鳴らし続けられ ドアをたたかれうるさすぎる これが合法なので警察も手が出せない 出ないと数分間ピンポンを鳴らし続けられ ドアをたたかれ これが毎週2回引っ越すまで3年間続きました 世間「国力が低下しているのに博士の職がありません」 政府「博士号所持者を増やします」 世間「教員不足なのに教員免許持ちの職がありません」 政府「単位減らして教員免許所持者を増やします」 政府はなぜ現在余っている人を活用しないのか不思議 ローマ字入力で「わたしわ～」とか「こんにちわ」 とか打っている人をみるとイライラします かな入力だと頭の中で発音しながら 「わたしは～」とか「こんにちは」と打てるので かな入力派です 日本の法律では 00:00(12:00am)は午後12:00 12:00(12:00pm)は午前12:00 だが 持っているデジタル時計の表示は 00:00(12:00am)は午前12:00 12:00(12:00pm)は午後12:00 で法律と違うがこちらが正しい と思うので法律をかえるべき ジャックダニエル シングルバレル セレクト 感想「すごく澄んだ味」 「経営者は社会人経験がないので常識がない」を 地で行く人がSNSを騒がせていますが 他の経営者は怒った方が良い ふつう1/1-12/31はすべての日を 24時間すべて含むと考えますので 1年を表しています -[00～01～11～][12～13～23～]24- -[am～～～～～ ][pm～～～～～][am -[12～01～11～][12～01～11～]12- 午後 ][午前～～～～   ][午後～～～～] 日本の午前午後は分かりにくいので法律を (am/pmに合わせて)変えた方が良いと思います 人間は自然の一部です 自然がなくなれば残るものは何もありません 当然人間も存在できなくなります 人間は...

2025/11/15(土)   等速円運動 (1回目)   ( uniform circular motion )   ■ 等速円運動 ▼ 位置R、速度V、加速度A θ :角度 ω:角速度 t :時間   R = (x, y, z) = (rcosθ, rsinθ, 0) = (rcosωt, rsinωt, 0)   V = (v x , v y , v z ) = ((d/dt)rcosωt, (d/dt)rsinωt, (d/dt)0) = (-ωrsinωt, ωrcosωt, 0) = (-ωy, ωx, 0)   A = (a x , a y , a z ) = ((d/dt)(-ωrsinωt), (d/dt)ωrcosωt, (d/dt)0) = (-ω 2 rcosωt, -ω 2 rsinωt, 0) = (-ω 2 x, -ω 2 y, 0)   ▼ ナブラ ∇ ∇ = (∂/∂x,∂/∂y,∂/∂z)   ▼ 回転rot rotR = ∇ ×R = (∂/∂x,∂/∂y,∂/∂z) ×(x, y, z) = (( ∂/∂ y)z - ( ∂/∂ z)y, ( ∂/∂ z)x - ( ∂/∂ x)z, ( ∂/∂ x)y - ( ∂/∂ y)x) = (0, 0, ( ∂/∂ x)rsinωt - ( ∂/∂ y)rcosωt) = (0, 0, 0)   rotV = ∇ ×V = (∂/∂x,∂/∂y,∂/∂z) ×(v x , v y , v z ) = (( ∂/∂ y)v z   - ( ∂/∂ z)v y , ( ∂/∂ z)v x   - ( ∂/∂ x)v z , ( ∂/∂ x)v y   - ( ∂/∂ y)v x ) = (0, 0, ( ∂/∂ x)ωx - ( ∂/∂ y)(-ωy)) = (0, 0, ω - (-ω)) = (0, 0, 2ω)   rotA = ∇ ×A = (∂/∂x,∂/∂y,∂/∂z) ×(a x , a y , a z ) = (( ∂/∂ y)a z   - ( ∂/∂ z)a y , ( ∂/∂ z)a x   - ( ∂/∂ ...