tonagai さんの自作式

10年間で脳卒中を発症する確率を計算します（国立がん研究センターによる）。

ヴァージニア・アプガーさんが考案した新生児の健康状態を判定するアプガースコア（APGAR score アプガーテスト、アプガー指数）を計算します。

感染症の流行を表す最もシンプルなモデル、SIRモデルを計算します。

SIRモデルに引き続き、感染症の流行を表すモデル、SEIRモデルを計算します。

NHKチコちゃんに叱られるで、関西大学・社会学部の保田先生が導出した方程式です。

国立研究開発法人医薬基盤・健康・栄養研究所の山田陽介室長らがアメリカやイギリス、オランダなどの研究者と共同で、科学雑誌「サイエンス」に発表した、人は一日に体内からどれくらいの水分を失うのか、予測式を計算します。

あなたの飼っている犬や猫の年齢を、人間に換算すると何歳になるか計算します。

余りが計算できるカシオの電卓MP-12Rが人気ということで、同じことができる自作式を作ってみました。

完璧なゆで卵を作るための方程式です。

数学セミナー2020年6月号の伊藤大雄さんの記事「雨に走れば公式　雨の中を歩くと走るで濡れる量の違い」に出てきた雨に走れば公式(Runnin' in the Rain Formula, RiR)を計算します。

から揚げなどの冷凍食品や赤ちゃん用のミルクのパッケージには個数とレンジ時間、体重とミルク量などが書いてありますが欲しい個数や体重が書いてないことがあります。

生年月日を一桁ずつ足したものが運命数。

浅草寺のおみくじは凶が多いと言われています。

乱数でズンとドコを5個発生します。

2021年共通テスト数学IA第二問はピッチとストライドについてでした。

Official髭男dismの115万キロのフィルムは、35㎜フィルムが1秒24コマで、フィルムを1秒当たり1.5フィート（1フィートは0.3048m)使うことから来ています。

2023年M-1グランプリ決勝戦での「さや香」のネタ、見せ算です。

出現確率が1%のガチャを100回やったら1回くらい当たるとか思ってない？ 全然違いますよ。

数当てゲームです。

生まれた月を4倍し、それに9を足し、それ全部を25倍し、それに生まれた日を足した数を計算してください。

あなたの年齢を当てます。

Steamで話題の『NKODICE（んこダイス）』もどきです。

5文字の英単語を当てるWordleにはまってますが、それを数字で行うゲームを作りました。

TVアニメ「しかのこのこのこのここしたんたん」のオープニング曲「シカ色デイズ」が流行っていますが、これをキヨシゲーム（ズン・ズン・ズン・ズンドコ）やNKODICE（んこダイス）のように乱数を使ったゲームにしてみました。

2つの自然数で、自分自身を除いた約数の和が同じになるものを友愛数といいます（博士の愛した数式でも登場しました）。

(x^2+9*x*y-y^2)^3 + (12*x^2-4*x*y+2*y^2)^3=(9*x^2-7*x*y-y^2)^3 +(10*x^x+2*y^2)^3 というラマヌジャンが発見した恒等式でx=1,y=0ならラマヌジャンが入院中にタクシーの ナンバー1729を聞いて1729=1^3+12^3 = 9^3+10^3と即座に答えたものが再現できます。

正則連分数を計算します（n項まで）

10進数をn進数(n＜＝34の任意の整数）に変換します。

海辺にn人の美女が並んでいる。

お好きな数を思い浮かべます（3桁か4桁がお勧め）。

フィボナッチ数列・ペル数列・パドヴァン数列・ペラン数列の計算と、その隣接項の比を計算します。

ジェイムズ・ジョイスのユリシーズにちなむ数列、a_n=(n^(n^n)の桁数）を計算します。

お好きな6つの数字を書きます。

好きな正の整数を思い浮かべて、その各桁を2乗したものを足していくと1になるか、 4,16,37,58,89,145,42,20のサイクルに落ちます。

積を素早く計算する方法として、和に直すために16世紀以前はProsthaphaeresisという、cos(a)cos(b)=(cos(a+b)+cos(a-b))/2が使われていたそうです。

11桁の数字からチェックディジットを計算します。

自然数 n の約数 d の x 乗の総和を値として持つ関数である約数関数σx(n)を計算します。

整数aとpとが互いに素であり、合同式 x²≡a (mod p)が解をもつとき、 aは p を法として平方剰余といいます。

Valeriy G. Narushin, et al. "Egg and math: introducing a universal formula for egg shape". Ann. N. Y. Acad. Sci., 2021.の紹介を彩恵りり (科学系ニュース解説アカウント)さんがTweetしていたのでカシオの高精度計算サイトでもやってみました。

任意の正の整数 n をとり、 n が偶数の場合、n を 2 で割る n が奇数の場合、n に 3 をかけて 1 を足す という操作を繰り返すと1になる、というコラッツ予想を図示します。

昔、イギリスのマクドナルドではチキンマックナゲットのセットに入っている個数は6個、9個、20個でした。

2022年共通テストの数学I・数学Aの問題4は不定方程式 （ディオファントス方程式） 5^4*x-2^4*y=1, 5^5*x-2^5*y=1, 11^5*x-2^5*y=1の解を求めるものでした。

日経サイエンス　2025年1月号の特集　和算再発見の佐藤賢一さんの記事、算聖　関孝和の実像に出てきた円弧の２乗の近似式と正確な値を比べてその精度を見ます。

筑波大学の高橋さんらが達成した2兆5769億8037万桁までの円周率計算に使われていたアルゴリズムであるAGM法（ガウス・ルジャンドルアルゴリズム、算術幾何平均法）を確認します。

円周率を計算する公式として知られている、ラマヌジャンの公式とチュドノフスキー兄弟の公式を比較します。

ペル方程式（Pell's equation） X^2 - N・Y^2 = ±1　の最小整数解(X,Y)と一つ大きい解を求めます。

ポアソン分布 P(X=n) = exp(-λ)λ^n / n! に従う乱数を生成します。

差が2であるような2つの素数を双子素数と言います。

56をひっくり返して65にして、足すと121のように回文数になります。

5次方程式をDKA法により計算します。

6次方程式をDKA法で計算します。

n次方程式、a(n)*x^n +a(n-1)*x^(n-1)+...+a1*x+a0 = 0をDKA法で計算します。

天に向かって続く数（加藤文元さん、中井保行さん著）に出てくるペレリマン数列（2乗するともとにもどる数）の計算です。

分割数p(n)は、自然数nを順序問わず自然数の和（自分含む）として何種類で表せるか？というものです。

Jesús Guilleraさんがとんでもなく早く収束する円周率Πの公式を見つけています。

dx/dt=f1(x,y,z,w,t) dy/dt=f2(x,y,z,w,t) dz/dt=f3(x,y,z,w,t) dw/dt=f4(x,y,z,w,t) という4変数の常微分方程式をルンゲクッタ法（4段4次）で計算します。

2変数のニュートン・ラフソン法を使って、f(x,y)=0, g(x,y)=0 を満たすx,yを計算します。

正の整数 n に対して、1 から n までの自然数のうち n と互いに素なものの個数をφ(n)と書き、オイラーのφ関数またはトーシェント関数と言います。

ランベルトのW関数は z=w*exp(w)の逆関数 w=W(z)です。

フーリエ級数を計算します。

アッカーマン関数Ack(m,n)を計算します。

合同数は、辺の長さがすべて有理数である直角三角形の面積のことですが、 その判定をします。

x座標が全て異なるn+1点を通るn次の多項式をラグランジュ補間で求めます。

Fermat's Libraryで出ていた以下を確認します。

Simon Plouffeさんが”A formula for the nth decimal digit or binary of π and powers of π”という論文で円周率πの任意の桁の値を計算する公式を出しています。

インドの弦理論の研究者Arnab Priya SahaさんとAninda Sinhaさんが論文Field Theory Expansions of String Theory Amplitudesの中で新しい円周率の公式を編み出しました。

平成26年センター試験　数学II・数学B　第6問は階乗（N!）の素因数分解を行う問題でした。

自然数M,NはM>Nであり、最大公約数が1であるとする。

京都・伏見の御香宮にある算額の問題です。

関数f(x)の1階微分・2階微分を計算してグラフにします。

N人のうち、k人以上の誕生日が一致する確率を計算します。

誕生日が一致する確率をグラフ化します。

2つの正規分布X1,X2の比（商）、X1/X2の分布を計算します。

ピアソンのχ^2検定（適合度の検定）を行います。

ピアソンのχ^2検定（独立性の検定）を行います。

1つの平均の検定を行います。

二群の比率の差の検定を行います。

対応のない場合の2群の差のt検定を行います。

1つの母平均の検定時に、効果量（Δ=(μ-μ0)/σ　平均の差が標準偏差の何倍か？）と有意水準を与えたとき、必要なサンプルサイズを計算します。

2つの平均の差の検定時に、効果量（Δ=(μ1-μ2)/σ　平均の差が標準偏差の何倍か？）と有意水準を与えたとき、必要なサンプルサイズn(n=n1=n2,それぞれのサンプルサイズの意味)を計算します。

検出力曲線（2つの平均の差μ1-μ2場合の1-β）を計算します。

検出力曲線（1つの平均の場合の1-β）を計算します。

相関係数の検定を行います。

単回帰分析 (y=ax+b)の回帰係数a,b（傾きと切片）のt検定と信頼区間を計算します。

選挙の出口調査の得票率から当選確実かどうか？を調べるような母比率の信頼区間を推定します。

コンドルセの陪審定理としても知られており、市川憂人さんの小説「神とさざなみの密室」でも出てきた正当性確率の計算をします。

y=(sqrt(cos(x))*cos(400*x)+sqrt (abs(x))-0.4)*(4-x*x)^0.1はハートを描く関数です(Googleのグラフ表示でおなじみ)。

人参をどの位置で切れば半分の体積になるか？のような問題、円錐台を2分割にする高さを求める計算式です。

三角形ABCの頂点座標を入力すると、三角形の五心（重心・内心・外心・垂心・傍心）の座標を計算します。

正多角形の頂点を計算します。

楕円の長軸半径、短軸半径、中心角を入力すると楕円の弧長L、弦長d、矢高hを計算します。

弧長（円弧の長さ）L、弦長d、矢高（円弧の高さ）h、半径rのどれか2つに値を入力して、残りの2つを0と入力すると（空白にはしないでください）、その残りの2つおよび中心角を計算します。

楕円の円周と一方の半径（長軸半径でも短軸半径でも可）からもう一方の半径を計算します。

超楕円（ラメ曲線） |x/a|^n + |y/b|^n = 1　(n=2の時は通常の楕円） の円周と面積を計算します。

相対論的量子力学の基礎方程式であるDirac方程式による水素原子(Z=1)のエネルギー準位を計算します。

楕円体 x^2/a^2 + y^2/b^2 + z^2/c^2 = 1 に対する反磁界係数　Nx, Ny, Nzを計算します。

3次元で星がランダムに分布しているとき、その重力の影響をある一点でみるときの分布がホルツマーク分布です。

球の抗力係数CD（Drag coefficient）をレイノルズ数Reを使って計算します。

カオスを示す有名なロジスティック写像 x_n+1 = a*x_n*(1-x_n)で 初期値が少し変わっただけでどのように値がずれていくかを見ます。

洪水などでドアのところまで水がやってくると、それほど水没している深さが深くなくても もうドアが開かなくなります。

標高が高くなると水の沸点が下がることは知られていますが、その近似計算を行います。

媒質（樹脂など）に粉などを混ぜたときの有効誘電率を計算します。

水と氷の複素誘電率の周波数特性（誘電分散）を計算します。

夏休みの宿題、提出する課題、申請書類、研究の進捗、などなど人間は締め切りギリギリまで先延ばしをします。

お湯を注ぐ音と水を注ぐ音が違って聞こえるのは、水の粘性が温度によって大きく異なるからだそうです。

新型コロナウイルスで話題になったパルスオキシメーターの原理、酸素と結びついたヘモグロビン（HbO2)と酸素を離したヘモグロビン(Hb)のモル吸光係数の波長による差をグラフとして表示します(計算ボタンを押した後、グラフボタンを押してください）。

電磁波の波長,反射率,透過率,電力半減深度,表皮厚みを計算します。

宇宙人がいる星がいくつあるか？というドレイク方程式を計算します。

PREM(Preliminary Reference Earth Model)による地球内部の密度分布を図示します。

PREM(Preliminary Reference Earth Model)による地球内部・外部の重力加速度を計算します。

電圧定在波比（VSWR）から反射損失（リターンロス）に変換する。

ストリップラインの特性インピーダンスを計算します。

反射損失（リターンロス）から電圧定在波比VSWRに変換する。

サーミスタのB定数をB=ln(Rb/Ra)/(1/(273.15+Tb)-1/(273.15+Ta))で計算します。

NTCサーミスタの温度特性を R = R0 * exp (B * (1/T - 1/T0)) によって計算します。

マイクロストリップラインの特性インピーダンスを計算します。

導波管入力標準ゲインホーンアンテナ（ピラミダルホーン）のアンテナゲイン(dBi単位）を計算します。

サーミスタの温度係数を 1/T = a + b*ln(R) + c*ln(R)^3 で計算します。

サーミスタの温度特性をSteinhart-Hart式 1/T = a + b*ln(R) + c*ln(R)^3 に基づいて、三点の温度・抵抗を入力するとa,b,cパラメータを算出します。

負荷インピーダンスと特性インピーダンスから反射係数を求めます。

部品などの寿命の解析に使われるワイブル分布 f(t)=(m/η)*(t/η)^m *exp[-(t/η)^m]のパラメータm,ηおよび平均寿命μを算出します。

円柱導体の高周波のインピーダンスを計算します。

矩形マイクロストリップアンテナの設計を行います。

導体は表面粗さで高周波の電気伝導率（抵抗率）が変わります。

低電圧でも人体に電流が流れると危険です。

T型回路のSパラメータを計算します。

π型回路のSパラメータを計算します。

HammerstadとJensenモデルでマイクロストリップラインの特性インピーダンスを計算したのち、導体損失(Morganモデルによる粗さ込み）と誘電体損失を計算し、特性インピーダンスの虚数部分も考慮してSパラメータに直しています。

ストリップラインの特性インピーダンスを計算したのち、導体損失(Morganモデルによる粗さ込み）と誘電体損失を計算し、特性インピーダンスの虚数部分も考慮してSパラメータに直しています。

レーダ方程式　Pr＝Pt*G^2*λ^2*σ/[(4π)^3*R^4]を用いて、レーダ最大検知距離を計算します。

負荷インピーダンスZL（ωの関数で入力）に対し、特性インピーダンスZo, 長さL、有効比誘電率εeff, 減衰定数αの伝送線路がついたとき、どういう入力インピーダンスになるかを計算します。

積層セラミックコンデンサ（MLCC)の市場での寿命予測に使われる経験式です。

水晶振動子の等価回路からインピーダンスの大きさ（Log単位）と位相（度単位）を計算します。

T型回路とπ型回路のインピーダンス３つを相互に変換します。

温度と湿度を高温・高湿にして加速試験を行う際、相対湿度を用いたモデルであるHallberg-Peckモデルを計算します。

インピーダンスが整合できる減衰器（アッテネータ）を設計します。

ABテストで必要なサンプルサイズを計算します。

パンサー尾形さんが出演しているNHK　笑わない数学の暗号理論に出てきたRSA暗号を検証します。

損益分岐点/損益分岐点比率/限界利益/限界利益率を計算します。

数学マーケティングに用いられる、購入回数の確率を表す負の二項分布モデル（NBDモデル）を計算します。

数学マーケティングに用いられる、購入回数の確率を表す負の二項分布モデル（NBDモデル）に出てくる確率分布を決めるKをニュートン・ラフソン法によって計算します。

数学マーケティングに用いられる、カテゴリー内の自社ブランド購入確率を表すディリシュレーNBDモデル（一般にはディリクレ）を計算します。

数学マーケティングに用いられる、カテゴリー内の自社ブランド購入確率を表すディリシュレーNBDモデル（一般にはディリクレ）を使って購入率と購入回数を計算します。

数学マーケティングに用いられる、カテゴリー内の自社ブランド購入確率を表すディリシュレーNBDモデル（一般にはディリクレ）に出てくるパラメータSを自社ブランドの購入率実績より算出します。

最近いつ買ったか？最近いつ訪れたか？を予測するガンマ・ポアソン・リーセンシー・モデルを計算します。

最近いつ買ったか？最近いつ訪れたか？を予測するガンマ・ポアソン・リーセンシー・モデルに現れる パラメータm,kを準ニュートン法で計算します。

イグノーベル賞を受賞したインド象の表面積を求める式です。

Daily Mailの記事、”Mathematicians create a formula for the perfect Christmas Day, and say the key is turning OFF the TV”からです。

1から順番に数を増やしていって3で割り切れる時はFizz,5で割り切れる時はBuzz,3でも5でも割り切れる時はFizzBuzzと言います。

女子の点数を低くするなどの不正で有名な東京医大ですが、その入試点数をシミュレーションします。

映画、アルキメデスの大戦で菅田将暉くんが演じる天才数学者、櫂直が導き出した鉄の総量と戦艦の建造費の関係式です。

調べたい年(デフォルトで今年)を入れると、節分で恵方巻きを丸かぶりする方角を表示します。