MATLAB,Python,Scilab,Julia比較 第5章 その18【複雑な定積分⑧】

バックナンバーはこちら。
https://www.simulationroom999.com/blog/compare-matlabpythonscilabjulia5-backnumber/

はじめに
登場人物
【再掲】フーリエ係数に至る道
【再掲】力業で解きたい数式
Juliaコード
処理結果
まとめ

はじめに

フーリエ係数に至る道。
今回は偶関数と奇関数を利用した数学パズルを力業で解く方法をJuliaで実現する。

登場人物

博識フクロウのフクさん

イラストACにて公開の「kino_k」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=iKciwKA9&area=1

エンジニア歴8年の太郎くん

イラストACにて公開の「しのみ」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=uCKphAW2&area=1

【再掲】フーリエ係数に至る道

太郎くん

まずは、フーリエ係数に至る道を再掲。

偶関数
奇関数
関数の内積
三角関数の加法定理
三角関数の積和公式
重要な極限値
三角関数の直交性
フーリエ係数

フクさん

今回は偶関数と奇関数を利用した数学パズルを力業で解く方法をJuliaで実現する。

【再掲】力業で解きたい数式

太郎くん

力業で解きたい数式は以下。

\(
\begin{eqnarray}
&&\displaystyle\int_{-2}^2\Big(x^3\cos(2x)+\frac{1}{2}\Big)\sqrt{4-x^2}dx\\
&=&\displaystyle\int_{-2}^2\Big(x^3\cos(2x)\sqrt{4-x^2}+\frac{1}{2}\sqrt{4-x^2}\Big)dx
\end{eqnarray}
\)

フクさん

無限次元ベクトルであることを前提として、
今回の定積分を表現すると以下になる。

\(
\begin{eqnarray}
\Bigg(
\begin{bmatrix}
-2^3&\dots&2^3
\end{bmatrix}\odot
\begin{bmatrix}
\cos(2(-2))&\dots&\cos(2(2))
\end{bmatrix}
\begin{bmatrix}
\sqrt{4-(-2)}\\
\vdots\\
\sqrt{4-(2)}\\
\end{bmatrix}\\+
\displaystyle\frac{1}{2}
\begin{bmatrix}
\sqrt{4-(-2)}&\dots&\sqrt{4-(2)}
\end{bmatrix}
\begin{bmatrix}
1\\
\vdots\\
1
\end{bmatrix}
\Bigg)\Delta x
\end{eqnarray}
\)

Juliaコード

フクさん

Juliaコードは以下。

using Printf

N=10000000; # 次元数
L=2-(-2);   # 定積分範囲
dt=L/N      # Δx

x=range(-2,2,length=N)'; # x軸作成(行ベクトルにするために転置)
y=(x.^3 .*cos.(x)*sqrt.((x.^2).+4)' + sqrt.(-x.^2 .+4)*ones(N,1)/2)*dt;
@printf("%.15f\n",y[1])