#=============================================================================================
# Project name: emtws01_test
# Project version: 1.0.0.4
#---------------------------------------------------------------------------------------------
# Generated by Aqualis (algorithm and equation analyzer for lightwave simulation)
# Aqualis version: 186.0.3.0
# Generated date: 2026/04/11 14:07:30
#=============================================================================================
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
i = 0
j = 0
i0001 = 0
z0001 = 0
z1001 = numpy.array([])
z1001_size = numpy.array([-1])
z2001 = numpy.array([[]], dtype=int)
z2001_size = numpy.array([-1,-1])
z2001_size[0] = 4
z2001_size[1] = 4
z2001 = numpy.zeros(z2001_size[0]*z2001_size[1], dtype=numpy.complex128).reshape(int(z2001_size[0]),int(z2001_size[1]))
z1001_size[0] = 4
z1001 = numpy.zeros(z1001_size[0], dtype=numpy.complex128)
for i in range(0, z2001_size[0], 1):
  for j in range(0, z2001_size[1], 1):
    #ハンケル関数計算テスト
    z0001 = jv(0, 1.0E-1*(i+2*j+1))-yn(0, 1.0E-1*(i+2*j+1))*1.0j
    z2001[i,j] = z0001
z1001[:]=0
for i in range(0, z1001_size[0], 1):
  for i0001 in range(0, z2001_size[1], 1):
    z1001[i] = z1001[i]+z2001[i,i0001]*i0001
z1001 = solve(z2001, z1001)
for i in range(0, z1001_size[0], 1):
  print("%27.12e%27.12e" %((z1001[i]).real,(z1001[i]).imag))
z1001_size[0] = -1
del z1001
z2001_size[0] = -1
z2001_size[1] = -1
del z2001