#=============================================================================================
# Project name: test 
# Project version: 1.0.0.68 */
#---------------------------------------------------------------------------------------------
# Generated by Aqualis (algorithm and equation analyzer for lightwave simulation) #
# Aqualis version: 185.0.3.0 #
# Generated date: 2025/06/19 12:56:32 #
#=============================================================================================
import numpy
import math
import cmath
import copy
import struct
import re
from scipy.linalg import solve, svd, eig, lu
from scipy.special import jv, yn
import sys
ic0001 = 0
ic0002 = 0
z0001 = 0
z10001 = numpy.array([])
z10001_size = numpy.array([-1])
z20001 = numpy.array([[]], dtype=numpy.complex128)
z20001_size = numpy.array([-1,-1])
z20001_size[0] = 4
z20001_size[1] = 4
z20001 = numpy.zeros(z20001_size[0]*z20001_size[1], dtype=numpy.complex128).reshape(int(z20001_size[0]),int(z20001_size[1]))
z10001_size[0] = 4
z10001 = numpy.zeros(z10001_size[0], dtype=numpy.complex128)
for ic0001 in range(0, z20001_size[0], 1):
  for ic0002 in range(0, z20001_size[1], 1):
    # ハンケル関数計算テスト 
    z0001 = jv(0, 1.0E-1*(ic0001+2*ic0002+1))-yn(0, 1.0E-1*(ic0001+2*ic0002+1))*1.0j
    z20001[ic0001,ic0002] = z0001
z10001[:]=0
for ic0001 in range(0, z10001_size[0], 1):
  for ic0002 in range(0, z20001_size[1], 1):
    z10001[ic0001] = z10001[ic0001]+z20001[ic0001,ic0002]*ic0002
# ===連立方程式の求解==================================================================== 
z10001 = solve(z20001, z10001)
# ===end 連立方程式の求解================================================================ 
for ic0001 in range(0, z10001_size[0], 1):
  print("%27.17e%27.17e" %((z10001[ic0001]).real,(z10001[ic0001]).imag))
del z10001
del z20001