前提・実現したいこと
pythonで数値解析の刻み幅をfor文で回して一斉に値を出力しようと思ったのですが,エラーが出ました.
発生している問題・エラーメッセージ
y_ans[n] = y IndexError: list assignment index out of range
該当のソースコード
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.3 y2 = 0.3 h_list = [2, 1, 0.5, 0.2, 0.1] for h in h_list: print('h = %d' %h) 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 a = 1/4 b = 3.141592653589793/6 c = 0.6 def func_f(x , y): return a * y + c * np.exp(a * x) * np.cos(x + b) def func_integ(x): return c * np.exp(a * x) * np.cos(x + b) 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
試したこと
特にないです.
補足情報(FW/ツールのバージョンなど)
ここにより詳細な情報を記載してください。
回答3件
あなたの回答
tips
プレビュー