perldl> use PDL::Complex
perldl> use PDL::LinearAlgebra

perldl> $matrixM = cplx pdl [ [ 1+1*i, 2+1*i], [ 1-2*i, 2-1*i] ];

perldl> p $matrixM

[
 [1 +1i  2 +1i]
 [1 -2i  2 -1i]
]


perldl> $matrixB_assigned_value = cplx pdl [ 5+8*i, 10-5*i ];

perldl> p $matrixB

[
 [ 5 +8i]
 [10 -5i]
]

perldl> msolve($matrixM,$matrixB)

[
 [3+1i]
 [2+1i]
]
 
[
 [   1   -2i     2   -1i]
 [-0.2 +0.6i   1.8 -0.4i]
]
 [2 2] 0
perldl> scalar msolve($matrixM,$matrixB)

[
 [3+1i]
 [2+1i]
]

perldl> $matrixM->minv x $matrixB

[
 [3     +1i]
 [2     +1i]
]