IndexError: индекс 10 за пределами оси 0 с размером 10

Question

Я численно настраиваю сетку сетки для x-сетки и x-вектора, а также сетку времени, но снова я установил массив для x (позиция), который должен быть только от 0 до 20, а t (время) будет от 0 до 1000, таким образом, чтобы решить уравнение теплоты. Но каждый раз, когда я хочу, в течение, например, я делаю несколько шагов 10, я получаю сообщение об ошибке:

"Traceback (most recent call last): 
File "/home/universe/Desktop/Python/Heat_1.py", line 33, in <module> 
x[i] = a + i*h 
IndexError: index 10 is out of bounds for axis 0 with size 10" 

Вот мой код:

from math import sin,pi 
import numpy 
import numpy as np 

#Constant variables 
N = int(input("Number of intervals in x (<=20):")) 
M = int(input("Number of time steps (<=1000):")) 

#Some initialised varibles 
a = 0.0 
b = 1.0 
t_min = 0.0 
t_max = 0.5 

# Array Variables 
x = np.linspace(a,b, M) 
t = np.linspace(t_min, t_max, M) 


#Some scalar variables 
n = []       # the number of x-steps 
i, s = [], []     # The position and time 

# Get the number of x-steps to use 
for n in range(0,N): 
    if n > 0 or n <= N: 
     continue 

# Get the number of time steps to use 
for m in range(0,M): 
    if m > 0 or n <= M: 
     continue 

# Set up x-grid and x-vector 
h =(b-a)/n 
for i in range(0,N+1): 
    x[i] = a + i*h 

# Set up time-grid 
k = (t_max - t_min)/m 
for s in range(0, M+1): 
    t[s] = t_min + k*s 

print(x,t) 

Answer 1

Вы индексировать вне диапазона:

for s in range(0, M+1): 
    t[s] = t_min + k*s 

Изменить на:

for s in range(M): 
    t[s] = t_min + k*s 

И это работает.

t Вы создаете с длиной M:

t = np.linspace(t_min, t_max, M) 

Таким образом, вы можете получить доступ только M элементы t.

Python всегда начинает индексировать ноль. Поэтому:

for s in range(M): 

будет делать M петли, в то время как:

for s in range(0, M+1): 

будет делать M+1 петли.