고등수학

원주율 계산 (가우스-루장드르의 산술-기하 평균법)

筑波大学の高橋さんらが達成した2兆5769億8037万桁までの円周率計算に使われていたアルゴリズムであるAGM法（ガウス・ルジャンドルアルゴリズム、算術幾何平均法）を確認します。 a0=1,b0=1/√2,t0=1/4,p0=1を初期値として、 an+1 = (an + bn) / 2 bn+1 = √(anbn) tn+1 = tn - pn * (an+1 - an)^2 pn+1 = 2*pn として、 π≒(an + bn)^2 / (4*tn) で計算しています。

원주율의 계산 (라마누잔과 추도노프스키)

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

GCJ 2008 - R1A - C

高精度計算を利用して、 Google Code Jam 2008 - Round 1A - Problem C. Numbers http://code.google.com/codejam/contest/dashboard?c=32016#s=p2 を力技で解いてみます。 n = 69 までは整数部の精度が保証されているので、 Small dataset（2 ≦ n ≦ 30）はこれで通過できるのかもしれません。

디리클레의 이타 함수

ディリクレのイータ関数です。

페르 방정식의 최소 해

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

포아송 난수

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

쌍둥이 소수

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

Lychrel 수 찾기

56をひっくり返して65にして、足すと121のように回文数になります。このような操作を繰り返して、例えば57→57+75=132→132+231=363で2回で回文数に、89は24回で、10911は55回で回文数になります。何回繰り返しても回文にならないのがLychrel numbersです。196が最小の候補。

5차 방정식 풀이 (복소수)

5次方程式をDKA法により計算します。係数は複素数も使えます。

6차 방정식 풀이(복소수)

6次方程式をDKA法で計算します。係数は複素数も使えます。

n차 방정식의 근(복소수, n은 20까지)

n次方程式、a(n)*x^n +a(n-1)*x^(n-1)+...+a1*x+a0 = 0をDKA法で計算します。 nは最大20までです。20より下の次数を計算するときは、不要な高次の 係数は0にしておいてください。係数は複素数も使えます。

실수의 근사 분수

Stern-Brocot木を探索し、与えられた正実数の近似分数を求めます。

사원수의 사칙연산

四元数

페렐리만 수열 (제곱해도 원래대로 돌아오는 수, 하늘을 향해 이어지는 수)

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

나눗셈의 수 (영화: 『기적이 준 수식』에서)

分割数p(n)は、自然数nを順序問わず自然数の和（自分含む）として何種類で表せるか？というものです。 ラマヌジャンとハーディの映画、奇蹟がくれた数式で大きく取り上げられていたので、その公開記念として 厳密な分割数p(n)とラマヌジャン・ハーディの漸近公式を計算します。nは250までです。

Jesus Guillera의 원주율 π 계산 공식

Jesús Guilleraさんがとんでもなく早く収束する円周率Πの公式を見つけています。試してみましょう。 ラマヌジャンの公式とチュドノフスキー兄弟の公式とも比べると面白いかも。

4변수 상미분방정식을 룬게-쿠타법(4단계 4차)으로 계산

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변수 뉴턴-러프슨법

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

오일러의 φ 함수 (토센트 함수)

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

램버트의 W 함수 (제k 해)

ランベルトのW関数は z=w*exp(w)の逆関数 w=W(z)です。主要解の他に第k解まで計算できます。 主要解の場合はk=0としてください。

푸리에 급수의 계산

フーリエ級数を計算します。関数f(x)（範囲は-L<=x<=L,　周期2L)を入力して係数を積分で求めます。

아커만 함수 Ack(m,n)

アッカーマン関数Ack(m,n)を計算します。mが大きくなるとものすごい巨大数になることが知られています。 残念ながらm=4,n=2までしか計算できません、、、（このサイトで扱える限界）

합수의 판별 (타넬의 정리)

合同数は、辺の長さがすべて有理数である直角三角形の面積のことですが、 その判定をします。判定したい数は平方因子を持たない（n=q*k^2と書けない）としてください。 タネルの定理（Tunnell's theorem)を使って判定しています。

실수 열린 구간의 근사 분수

与えられた正実数開区間の近似分数を探索します。

라그랑주 보간법(최대 11개 점)

x座標が全て異なるn+1点を通るn次の多項式をラグランジュ補間で求めます。 入力する点は、nが11より小さくて不要な場合でも0を入れておいてください。

삼각함수 표와 큰 그래프

三角関数 sin(x)とcos(x)とtan(x)の表とグラフ

리우빌의 약수의 약수의 개수가 합의 제곱이 3제곱의 합이라는 정리

Fermat's Libraryで出ていた以下を確認します。 Liouville proved the following property of the divisors of the divisors of a number: 1) Start with a number (e.g. 10) 2) Write down its divisors (1,2,5,10) 3) Write down the number of divisors of each divisor (1,2,2,4) 4) (1+2+2+4)²=81=1³+2³+2³+4³

구면 삼각형의 면적

球面三角形の面積を求める計算ができます。

1/(1+x²) (x=[0,1])의 그래프와 원

y=1／(1+x^2) と、円の方程式　のx=[0,1]区間の積分値はどちらも、π/4（1/4円)となり、等しくなります

원주율 π의 임의의 자릿수 값 (Simon Plouffe 제공)

Simon Plouffeさんが”A formula for the nth decimal digit or binary of π and powers of π”という論文で円周率πの任意の桁の値を計算する公式を出しています。これを確認します。ただこのサイトでは50桁までが限界です。

임의의 피타고라스 수 생성

入力１と入力２を指定することにより任意のピタゴラス数を生成します。ただし、m>nはでなくてはいけません。

임의의 피타고라스 수 생성 (모든 패턴)

ピタゴラス数を表に表示

n!의 끝자리에 0이 몇 개나 있는지

n!の末尾に0がいくつあるか

소수 판별!

素数判定

현상 이론에서 탄생한 새로운 원주율 공식

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

실수의 유리수 근사

実数の有理数近似

2차 정사각 행렬·실수 대칭·뉴턴 근사

2次正方行列が実対称行列であるときの数値的対角化です。 2次正方行列であり、対称行列なので、「絶対値が最大の非対角要素」を探索する手順が省けるので作業はシンプルになります。 また、『cos2θとsin2θ』の一次式に直すと周期も圧縮され二分法による区間発見と短縮も良くなります。 小さくした区間内でニュートン近似します。 精度は区間分割数を毎回100、区間の幅を0.1^15、ニュートン近似では0.1^20に落してあります。

2차원 벡터에 회전 행렬을 곱합니다.

２次ベクトルに回転行列を掛ける計算をします。 A,Bは２次元ベクトルであるとします。 Rは２次正方行列であり、回転行列です。 B=R*Aを計算します。

3개의 2차 정사각 행렬의 곱

2次正方行列の三つの積を計算します。 三つの2次正方行列A,B,Cの積を AB=A*B ABC=(A*B)*C と計算し、結果を出力します。(左に右を掛けます。)

2차 정사각 행렬을 이용한 2차 식 계산

2次正方行列xでの、2次式の計算をします。 2次正方行列による多項式を考えます。 2次正方行列xについての2次の整式です。 式の形式をy=aL*x*x*aR + bL*x*bR + cとします。 (*は行列同士の積の記号、+は行列同士の和。 積は左の行列に右の行列を順に掛けます。) aL,aR,bL,bR,c,x,yは全て2次正方行列です。

2차 정사각 행렬에 회전 행렬을 곱합니다.

2次正方行列に回転行列を掛ける計算をします。 2次正方行列Aを回転行列Pの逆行列Qに右から掛け行列QAを計算し、 次に回転行列PをQAに右から掛け、QAPを計算します。 Pの要素は(Qの右の数字は行番号・列番号)、 P11=cosθ P12=-sinθ P21=sinθ P22=cosθ とします(θはラジアン)。

2차 정사각 행렬의 3차 식 계산

2次正方行列の3次式を計算します。 y,x,aL,aR,bL,bR,cL,cR,dは2次正方行列とします。 式の形式を y=aL*x*x*x*aR + bL*x*x*bR + cL*x*cR + d とします。 「普通の数」のxについての3次式y=ax^3+bx^2+cx+dを真似しています。 L,Rはxのべきに対して左、右の位置を意味します。 *は行列同士の積で、左の行列に右の行列を掛けます。 複数の*結合は順次、左に対して右を掛けて行く事に成ります。 最後に各項を足します。

2차 정사각 행렬의 고유방정식의 계수와 해를 구합니다.

2次正方行列の固有方程式を計算します。 固有方程式の係数を表示します。 固有方程式は2次方程式ですが、判別式Dの値を計算します。 また、判別式の平方根の値を表示します。 固有方程式の解を出力します。 重解のときも解を２つ表示します。判別式の値が0であれば重解です。

3차 정사각 행렬 3개의 곱

3次正方行列の三つの積を計算します。 三つの3次正方行列A,B,Cの積を AB=A*B ABC=(A*B)*C と計算し、結果を出力します。(左に右を掛けます。)

3차 정사각 행렬과 관련된 계산을 해보겠습니다. 그 5

3次正方行列に関連した計算をします。 3次正方行列の固有方程式の係数書き出しと、解の計算を行います。 固有方程式の形式をa(λ^3)+b(λ^2)+cλ+d=0とします。 (1)3次正方行列の固有方程式の係数を求めます。 (2)λの最高次数(3)の係数(-1)で各係数を割った処理をし、表示します。 (3)固有方程式の近似解を表示します。 (4)検算として式に近似解を代入した結果等を表示します。 (5)係数b,c,dが整数のとき因数定理を使い、解が見つかった場合には結果を表示します。

2차 정사각 행렬의 제곱근을 구합니다.

２次正方行列の平方根を計算します。 ある２次正方行列Mが与えられたとき、M=X^2を満たす２次正方行列Xを求めます。ここではXをMの平方根(2乗根)と呼ぶことにします。普通の数のルート(根号記号√)に当たる演算を考えることになります。

3차 정사각 행렬과 관련된 계산을 해보겠습니다. 그 4.

３次の正方行列に関連した計算をします。 1)固有方程式の係数を出力します。固有方程式の形式はk0λ^3+k1λ^2+k3λ+k4=0とします。 2)小行列式の値を計算します。 3)余因子の値を計算します。 4)余因数行列を計算します。 5)サラスの図により行列式の値を計算します。 6)サラスの図により計算した行列式の値が0でなければ逆数を計算します。

4차 정사각 행렬과 관련된 계산을 해보겠습니다.

４次正方行列に関連した計算をします。 ここでは小行列式の値、余因子、余因数行列を計算します。 また、参考に入力した４次正方行列の行列式の値とその逆数を出力します。

벡터의 정규화(2차원)

二次元ベクトルを単位ベクトルに規格化(正規化)します。 規格化前のベクトルvの要素をv1,v2とします。 ベクトルの大きさをrとします。 規格後のベクトルuの要素をu1,u2とします。

3개의 4차 정사각 행렬의 곱을 계산합니다.

4次正方行列3個の積を計算します。 三つの4次正方行列A,B,Cを考え、まずAB=A*Bを計算し、次にABC=ABC=AB*Cを計算します。 AB,ABCも4次正方行列です。結果であるABCを出力します。(*は行列の積)

2행 2열 행렬의 행렬식을 구합니다. (요소가 복소수 분수인 경우)

２行２列行列式の計算をします。 ここでは複素数の分数が要素である場合を扱います。 aは1行目1列目の要素とします。 bは1行目2列目の要素とします。 cは2行目1列目の要素とします。 dは2行目2列目の要素とします。 入力ではa,b,c,dは分数であるとして入力します。 a,b,c,dの分母には0でない数値を入力します。 a,b,c,dの分子・分母とも複素数の範囲で入力します。

벡터의 정규화(3차원)

三次元ベクトルを規格化(正規化)します。 規格化するベクトルvの成分をv1,v2,v3とします。 vの大きさをrとします。 規格化したベクトルuの成分をu1,u2,u3とします。

벡터의 정규화(4차원)

四次元ベクトルを規格化(正規化)します。 規格化するベクトルvの成分をv1,v2,v3,v4とします。 vの大きさをrとします。 規格化したベクトルuの成分をu1,u2,u3,u4とします。

3차 방정식의 계수를 구합니다. (분수, 복소수)

３次固有方程式(特性方程式)の係数を計算します。 ３次正方行列での固有値問題を考えるので ３行３列の行列式からの固有方程式導出となります。 ３次固有方程式の形式を -(λ^3)+m*(λ^2)+n*(λ^1)+定数項とします。 最高次数(λ^3)項の係数は-1のままとします。 (λ^2)項の係数、(λ^1)項の係数、定数項を計算し出力します。 固有方程式の形式は展開した後の形で表示します。

소수 생성(밀스의 상수)

任意の自然数ｎを入力してください。

소수 개수 (윌슨의 정리)

任意の自然数nを入力してください。 16以上の数は対応しておりません。

소수 생성식(윌슨의 정리)을 만들려고 한다면

a

2차 함수의 정점

2次関数の頂点を求めます。

스칼라 3중 곱을 계산합니다.

スカラー三重積の計算をします。 xyz空間座標系を右手系で考えます。 ベクトルbとcの外積をbcとします。 aとbcの内積を計算し、スカラー三重積の値を求めます。

코사인 정리 - 수학Ⅰ

a,b,cosの値を用いてcの値を求めます

3차 실수 대칭 행렬의 대각화

実対称行列について数値的対角化法を行います。 ここでは３次正方行列の実対称行列を対角化します。

벡터의 내적에 대한 코사인 cosθ의 값.

ベクトルa,bの内積の、余弦cosθの値を計算します。

벡터의 내적에서 cosθ의 값(3차원 벡터)

二つのベクトルa,bの内積の余弦cosθの値を計算します。

4차 다항식의 계수

4次固有多項式の係数導出をします。 式の形式をk4*λ^4+k3*λ^3+k2*λ^2+k1*λ+k0 とします。

벡터의 방향 코사인 계산

３次元ベクトルの方向余弦を計算します。 ベクトルをaとします。始点は原点です。

벡터의 원점을 지나는 직선 위의 수직 투영

ベクトルの正射影の計算をします。

n차원 구의 부피와 표면적

n次元球体の体積の一般化です

고유방정식의 계수 도출(5차)

固有方程式の係数導出を行います。 ここでは五次正方行列の場合を扱います。 (五次の特性方程式)

사원수의 n제곱을 계산합니다.

四元数のn乗を計算します。 式の形式をq=a+bi+cj+dk, y=q^nとします。 i,j,kは各虚数単位です。 a,b,c,dは実部・各虚部です。

사원수 계산을 해보겠습니다.

四元数の計算をします。 ２つの四元数Q1とQ2の絶対値、逆元、積Q1Q2,Q2Q1を計算します。

레비-치비타 기호(3자)

レヴィ・チヴィタの記号の値を表示します。 ここではεとします。 数字のセットを{1,2,3}とします。 添字i,j,kによるセット{i,j,k}へ置換されるとき、 偶置換か、奇置換か、その他か、で値が変わります。 εの値とi,j,kを並べて列挙します。

Levi-Civita의 기호(4자)

レヴィ・チヴィタの記号の値を表示します。 ここではεとします。 数字のセットを{1,2,3,4}とします。 添字t1,t2,t3,t4によるセット{t1,t2,t3,t4}へ置換されるとき、 偶置換か、奇置換か、その他か、で値が変わります。 εの値とt1,t2,t3,t4を並べて列挙します。 但し、ε=+1,-1の場合だけ表に出力し、ε=0の組み合わせは省略します。

Levi-Civita 기호(5자)

レヴィ・チヴィタの記号の値を表示します。 ここではεとします。 数字のセットを{1,2,3,4,5}とします。 添字t1,t2,t3,t4によるセット{t1,t2,t3,t4,t5}へ置換されるとき、 偶置換か、奇置換か、その他か、で値が変わります。 εの値とt1,t2,t3,t4,t5を並べて列挙します。 但し、ε=+1,-1の場合だけ表に出力し、ε=0の組み合わせは省略します。

행렬식의 부호와 첨자(3차)

行列式展開時の符号と添字を出力します。 三次の場合です。 整頓したものは「サラスの公式」として知られています。 符号と要素の添字n・σ(n)を列挙します。 n=1,2,3です。

행렬식의 부호와 첨자(4차)

行列式展開時の符号と添字を出力します。 四次の場合です。 符号と要素の添字n・σ(n)を列挙します。 n=1,2,3,4です。

행렬식의 부호와 첨자(5차)

行列式展開時の符号と添字を出力します。 五次の場合です。 符号と要素の添字n・σ(n)を列挙します。 n=1,2,3,4,5です。

사원수를 이용한 1차 식의 계산

四元数による一次式の値を計算します。 式の形式を y=aL*x*aR+b とします(*は積の記号)。 y,x,aL,aR,bは四元数です。 積は左に対して右を掛けて行きます。

사원수를 이용한 이차식 계산

四元数による二次式の値を計算します。 式の形式を y=(aL*x*x*aR)+(bL*x*bR)+c とします(*は積の記号)。 y,x,aL,aR,bは四元数です。 積は左に対して右を掛けて行きます。

사원수를 이용한 3차 식의 계산

四元数による三次式の値を計算します。 Lnはx^nの左側の係数です。 Rnはx^nの右側の係数です。 Mは定数項です。 積は左に対し右を掛けて行きます。 各n次項では、Lnにxを右からn回掛け、さらに右からRnを掛けます。 最後に定数項も含めて和算を行います。 これを多項式の値として出力します。

8원수의 곱을 계산합니다(두 개).

八元数の積を計算します。 二つの八元数a,bの積の計算です。

sinc 함수

1. 定義と役割sinc（シンク）関数は、\(\frac{\sin (x)}{x}\) で定義される関数です。通信理論や信号処理（サンプリング定理）において、波形を滑らかに補完するために不可欠な数式です。 2. プログラムの要点 0除算の回避: \(x=0\) のとき、数式上は \(\frac{0}{0}\) となり計算不能（エラー）になります。極限値の採用: 数学的に \(\lim _{x\rightarrow 0}\frac{\sin (x)}{x}=1\) であるため、コード内では if 文を使って \(x=0\) のときに強制的に \(1\) を返すように制御しています。 3. 挙動 \(x=0\) で最大値 \(1\) をとります。\(x\) が離れるにつれ、波打ちながら \(0\) に収束していきます。

Dilogarithm(이중 로그 함수) ※실수(x≦1)에만 해당

Dilogarithm(二重対数関数)の値を計算します。 x∈(-∞,1]で計算可能です。 x>1ではLi(2,x)は実数にはならず、実装していないので、エラーとして0が返ります。

감마 함수 수치표 출력

ガンマ関数の数値表を出力します。 複素平面をイメージしてください。 ガンマ関数への入力数値sは複素平面の線分上から行います。 始点となる点Pの複素数の実部は、0でない正の数値であるとします。 終点となる点Qの複素数の実部も、0でない正の数値であるとします。 線分をn等分し、始点を第０番目、終点をn番目の入力として、 ガンマ関数の数値を計算します。 出力形式は表形式です。