Блог

Я пытаюсь изменить этот Matlab для цикла:

 n =100;
u_lsg = zeros (n ,n) ; 
for i =1:1: n
    for j =1:1: n
        u_lsg(i , j) =( i*h)*(j*h )*(1 -( i*h ))*(1 -( j*h ));
    end
end
 

в петлю питона.

Что у меня есть до сих пор:

 n =100
u_lsg = np.zeros((n ,n))
for i in range (1,n 1,1):
    for j in range (1,n 1,1):
        u_lsg[i,j]= (i*h)*(j*h)*(1-(i*h))*(1-(j*h))
 

Но это дает IndexError: index 100 is out of bounds for axis 1 with size 100 .

Что я сделал не так? И из того, что я понимаю end , в python нет команды, но тогда как закончить цикл?

Matlab 1-indexed (индексы массива начинаются с 1 ), Python (и NumPy) 0-indexed (индексы массива начинаются с 0 ), поэтому просто измените диапазоны соответственно:

 n = 100
u_lsg = np.zeros((n, n))
for i in range(n):
    for j in range(n):
        u_lsg[i, j] = (i*h)*(j*h)*(1-(i*h))*(1-(j*h))