computes the solution to a complex system of linear equations A*X=B, using LU decomposition with partial pivoting and row interchanges
function C_solve2 import Complex; input Complex A[:, size(A, 1)]; input Complex B[size(A, 1), :]; output Complex X[size(A, 1), size(B, 2)]; end C_solve2;