【python】pythonによる科学計算を行った結果をcsvやtextで出力したい

Question

実現したいこと

下記のようなpythonによる科学計算を行った結果をcsvやtextで出力したいです。
良い方法はありますでしょうか。

実行したコード

import scipy as sp
import matplotlib as mpl
import matplotlib.pyplot as plt
from matplotlib.pyplot import plot,show,xlabel,ylabel,title,legend,grid,axis,rcParams,tight_layout
from scipy import real,imag,pi
from RI import WLx, epAg, epAu, RIAu, RIAg　　　#別ファイルで定義した屈折率データ

n1 = 1       # 周辺媒質の屈折率
n2 = RIAg    # 粒子の屈折率
r=25         # 粒子の半径
k = 2 * pi / WLx    # 波数の配列
alpha = 4 * pi * (r**3) * (n1**2) * (n2**2 - n1**2) / (n2**2 + 2 * n1**2)  # 分極率の計算
Csca = k**4 / (6 * pi) * abs(alpha)**2   # 散乱断面積
Cabs = k * imag(alpha)                   # 吸収断面積
Qsca = Csca / ((r**2) * pi)              # 散乱効率
Qabs = Cabs / ((r**2) * pi)              # 吸収効率

plt.figure(figsize=(8,6))        # figure size
plot(WLx,Csca/10**4, label=r"$C_{{\rm sca}}$",linewidth = 3.0, color='black')        # real part of the refractive index of gold
plot(WLx,Cabs/10**4, label=r"$C_{{\rm abs}}$",linewidth = 3.0, color='gray')   # real part of the refractive index of gold
xlabel("Wavelength (nm)",fontsize=22)       # x-axis label
ylabel(r"Csca, Cabs(×10000 nm^2)",fontsize=22)          # y-axis label
title(r"Csca, Cabs",fontsize=22)   # Title of the graph 
grid(True)                                 # Show Grid 
axis([300,700,0,5])                       # Plot Range
plt.tick_params(labelsize=20)  # scale fontsize=18pt
legend(fontsize=20,loc='lower right')
tight_layout()
show()

上記コードの結果

![]

Answer 1

具体的に何を保存したいのかわかりませんが、たとえば下記のようなやり方があります。
csvモジュール以外にも、numpyやpandasなど色々ありますので、調べてみて下さい。

import csv
with open("result.csv", 'a', newline="") as f:
        writer = csv.writer(f)
        writer.writerow(WLx)
        writer.writerow(Csca)

Comments

[tallman0202]

質問が不明瞭で申し訳ありません。
教えていただいたコードでcsv化できました。
pandasのデータフレーム等も検討してみようと思います。
ご回答有難うございました。

Answer 2

とりあえず、科学計算を行った結果をprint()の中に入れて実行してみればどうでしょう

Comments

[tallman0202]

教えていただいた通りprint(Csca)をしてみたところ、数値が出ました。
printするだけでも数値が得られるのですね。
ご回答有難うございました。