F(x)の各xでの値をtxtに書き込みたい

Question

関数F(x)=x∫[x~∞] special.kv(5/3,y)dy　の各xの値をtxtに書き込みたい。

from scipy.special import kv
import matplotlib.pyplot as plt
from scipy.integrate import quad
import numpy as np

xs = np.linspace(0, 10, 1000)
f = lambda z: kv(5/3,z)
F = [quad(f,x,np.inf)[0]*x  for x in xs]

x=0.0
h=0.001

with open("test3.txt","w") as f:
    while x <= 10:
        x += h
        print("{:,5f} {:7f}".format(x,F),file=f)
        print(x,F)

エラー

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-22-e3787445c636> in <module>
     14     while x <= 10:
     15         x += h
---> 16         print("{:,5f} {:7f}".format(x,F),file=f)
     17         print(x,F)

ValueError: Invalid format specifier

追加変更

from scipy.special import kv
import matplotlib.pyplot as plt
from scipy.integrate import quad
import numpy as np

xs = np.arange(0.0,10,0.1)
f = lambda z: kv(5/3,z)
F = [quad(f,x,np.inf)[0]*x for x in xs]


with open("test4.txt","w") as f:
        print(xs,F,file=f)

でtxtに起こせはするけれど,
[x=0 0.1 0.2~][F=0.91 0.90 0.89~]
になって見づらいので
[X=0.1 F=0.90 x=0.2 F=0.89 ~]
にしたいです。

Answer 1

まず提示エラーの原因は誤字です。{:,5f}→{:.5f}としてください。
次に、

[X=0.1 F=0.90 x=0.2 F=0.89 ~]

にしたいです。

であれば以下のように書けばよいです。

    for x, fval in zip( xs, F):
        print("{:.5f} {:7f}".format(x,fval),file=f)

Comments

[hataki7]

なるほど。とても勉強になります。ありがとうございます。