\(y=(x+iz)^2+1\) の 虚数解


判りにくいが
・ピンク :\(y=(x+iz)^2+1\) の実部
・オレンジ:\(y=(x+iz)^2+1\) の虚数部
・灰色  :\(y=0\) 平面
・青い球 :\(y=(x+iz)^2+1=0\) となる点(=虚数解)

作り方

[1] 実数部曲面を作成する

・c pink ・・・ 色をピンクに設定
・p surface np.linspace(-1.5,1.5,30) np.linspace(-1.5,1.5,30) ((x+1j*z)**2+1).real
 ※ numpy では虚数単位は 1j , 実数部は ().real と記す。
・surface

[2] 虚数部曲面を作成する

・c orange ・・・ 色をオレンジ色に設定
・p surface np.linspace(-1.5,1.5,30) np.linspace(-1.5,1.5,30) ((x+1j*z)**2+1).imag
 ※ numpy では虚数部は ().imag と記す。
・surface

[3] \(y=0\) 平面を作成する

・c white ・・・ 色を白(灰色)に設定
・p surface np.linspace(-1.5,1.5,30) np.linspace(-1.5,1.5,30) np.zeros((30,30))

[4] 虚数解に球を置く

・c blue ・・・ 色を青に設定
・p clear
・p 0 0 1
・p 0 0 -1
・distribute sphere 0.15