関数型プログラミング言語 F# (fsharp)の数値計算ライブラリ Math.NET Numericsを使う(1) 複素行列を定義して一次方程式や逆行列、行列式などを計算する。Polyglot Notebooksで扱うのが便利。

関数型プログラミング言語をやってみようと思ってF#を容易にJupyter Notebook形式で使えるPolyglot Notebooksを導入した。

さっそく使ってみる。実はC#より簡単に数値計算ライブラリMath.NET Numericsが使えるのだった。

https://numerics.mathdotnet.com/

C#のときにやった例題をF#でもやっていく。

まずはこれ。

Visual C# (C_sharp)の数学ライブラリ Math.NET Numericsを使う(1) 複素行列を定義して一次方程式や逆行列、行列式などを計算する。 

Polyglot NotebooksにNuGetで入れるにはこうする。

#r "nuget: MathNet.Numerics.FSharp, 5.0.0"

そしてopenしてから、

open MathNet.Numerics.LinearAlgebra

行列とベクトルを定義する。

let A = matrix [[ 2.0; 1.0; 1.0]

                      [ -1.0; 1.0; -1.0]

                      [ 1.0; 2.0; 3.0]] : Matrix<Complex>

let b = vector  [2.0; 3.0; -10.0] : Vector<Complex>

最後のMatrix<Complex>などが無ければDouble型になる。

連立方程式の計算、逆行列、行列式は以下で計算できる。

let x = A.Solve(b)

A.Inverse()

A.Determinant()

まとめてこちらに結果を。

黄色い波線が数字にあるのはDoubleで書いているのにComplexと定義しているからで、複素数を入れれば消えます。

C#より完結に書けるのはいいな。