前提・実現したいこと
エラーが出ました.
ソースが長くなってしまって申し訳ないです.
発生している問題・エラーメッセージ
sp = scipy.interpolate.InterpolatedUnivariateSpline(x_ans , y_ans) ^ SyntaxError: invalid character in identifier
該当のソースコード
import math import numpy as np import matplotlib.pyplot as plt import scipy.interpolate EPS = 0.00000001 x = 0.0 xmax = 10.0 y1 = 0.0 y2 = 0.0 h = 2.0 num = int((xmax - x) / h) + 1 x_ans = [0.0] * num y1_ans = [0.0] * num y2_ans = [0.0] * num y_ans = [0.0] * num n = 0 def func_f(x , y): return 2.0 * x def func_integ(x): return x * x print(" X Y1 Y2 Y E1 E2") while (True): y = func_integ(x) y_ans[n] = y if (y < EPS and y > - EPS): print('%8.4lf\t%8.10lf\t%8.10lf\t%8.10lf\t%8.20lf\t%8.20lf' %(x, y1, y2, y, 0.0, 0.0)) else: print('%8.4lf\t%8.10lf\t%8.10lf\t%8.10lf\t%8.20lf\t%8.20lf' %(x, y1, y2, y, math.fabs((y1 - y) / y) * 100.0, math.fabs((y2 - y) / y) * 100.0)) #Euler method y1_ans[n] = y1 y1 += h * func_f(x , y1) #Runge-Kutta method y2_ans[n] = y2 k1 = func_f(x , y2) k2 = func_f(x + h / 2.0 , y2 + k1 * h / 2.0) k3 = func_f(x + h / 2.0 , y2 + k2 * h / 2.0) k4 = func_f(x + h , y2 + k3 * h) y2 += (h / 6.0) * (k1 + 2.0 * k2 + 2.0 * k3 + k4) x_ans[n] = x x += h n = n + 1 if (x > xmax): break sp = scipy.interpolate.InterpolatedUnivariateSpline(x_ans , y_ans) sx = np.linspace(0.0 , xmax , 1000) sy = sp(sx) plt.plot(sx , sy , color = 'blue' , linestyle = '-') sp = scipy.interpolate.InterpolatedUnivariateSpline(x_ans , y1_ans) sx = np.linspace(0.0 , xmax , 1000) sy = sp(sx) plt.plot(sx , sy , color = 'green' , linestyle = '--') sp = scipy.interpolate.InterpolatedUnivariateSpline(x_ans , y2_ans) sx = np.linspace(0.0 , xmax , 1000) sy = sp(sx) plt.plot(sx , sy , color = 'red' , linestyle = '-.') plt.grid() plt.show()
試したこと
なにも試していません
補足情報(FW/ツールのバージョンなど)
ここにより詳細な情報を記載してください。
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2020/07/06 16:28