Pythonにてforループの計算で出たエラーを解決したい。

Question

Pythonにてforループの計算で# Q[i] の式を時間指定でtemp_new[i]の式に入れて計算したいのですが、以下のエラーが見つかりました。
何か解決策等あれば教えて頂けるとありがたいです。

(base) C:\Users\knishimura\Desktop\レーザ熱分布\00>python 00.py
File "00.py", line 48
for n in range(1, nt+1):
^
SyntaxError: invalid syntax

import matplotlib.animation as animation
import numpy as np
import math

# input parameter
den = 2330.0
cp = 678.0
cond = 83.7
temp_init = 20.0
lx = 100e-6
nx = 101
tend = 400e-9  
dt =   0.000000001  #1.0e-9
tout = 0.000000005  #5.0e-9
a=300
alpha = cond / (den * cp)
dx = lx / (nx - 1)
nt = int(tend / dt)
nout = int(tout / dt)

#initial condition 
temp = np.full(nx, temp_init) 
time = -200.0
temp_new = np.zeros(nx)

# Boundary condition
temp[0] = temp[1] 
temp[nx-1] = temp[nx-2] # Neumann @ x=Lx

# graph data array
ims = []
fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
gx = np.zeros(nx)
for i in range(nx):
    gx[i] = i * dx


# FTCS
Q = np.zeros(nx)
for i in range(1, nx-1):
    # Q[i]
    Q[i] = (math.exp(-(dt/50e-9)**2)*math.exp(-a*dx*i))/(den*cp*(50e-9)*math.sqrt(math.pi))  #単位 [℃]

for n in range(-199, nt+1):
    for i in range(1, nx-1):
        temp_new[i] = Q[i]+temp[i] + dt * alpha * (temp[i+1] - 2.0 * temp[i] + temp[i-1]) / (dx * dx)
       
    # update
    for i in range(1, nx-1):
        temp[i] = temp_new[i]

   
    # Boundary condition
    temp[0] = temp[1] 
    temp[nx-1] = temp[nx-2] # Neumann @ x=Lx
   

    time += dt

    if n % nout == 0:   
        print('n: {0:7d}, time: {1:8.1e}, temp: {2:10.6e}'.format(n, time, temp[nx-1]))
        im_line = ax.plot(gx, temp, 'b')
        im_time = ax.text(-200, 200, 'Time = {0:8.1e} [s]'.format(time))
        ims.append(im_line + [im_time])  

# graph plot
ax.set_xlabel('x [m]')
ax.set_ylabel('Temperature [C]')
ax.set_xlim(1e-6, 100e-6)
ax.set_ylim(0, 2500)
ax.grid()
anm = animation.ArtistAnimation(fig, ims, interval=200)
anm.save('animation.gif', writer='pillow')
plt.show()

Answer 1

前の行の括弧が閉じてないので、式の続きと見なされて文法エラーです。

Comments

[ab-3]

回答ありがとうございます。
Q[i]のところで訂正してみましたが、同様のエラーがでました。
訂正箇所が間違っていますでしょうか？

[otn]

・訂正の仕方が間違っている
・訂正の仕方は正しいが別のエラーが出てそれを「同様のエラー」と勘違いしている
などいろいろ、考えられます。

[ab-3]

回答ありがとうございます。＃FTCS内のtemp_new[i]の式を初期時間の-200nsから200nsまで時間幅1nsで計400s間計算したいのですが、48行目のfor n in range(-199, nt+1):のところでSyntaxError: invalid syntaxが出ています。現在のコードでなにか間違い等見受けられたら教えて頂きたいです。

[otn]

現時点で質問文に書いてあるコードだと、そのエラーにはならないので、実行していると思っているコードが違うんでしょうね。

Answer 2

そのforの前の行、カッコの数があってません

Comments

[ab-3]

回答ありがとうございます。
Q[i]のところで訂正してみましたが、同様のエラーがでました。
訂正箇所が間違っていますでしょうか？

[y_waiwai]

そのエラーを提示しましょう。
エスパーならあなたのやらかしてることをわかるんでしょうけど、そーではありませんｗ

[ab-3]

ご指摘ありがとうございます。＃FTCS内のtemp_new[i]の式を初期時間の-200nsから200nsまで時間幅1nsで計400s間計算したいのですが、48行目のfor n in range(-199, nt+1):のところでSyntaxError: invalid syntaxが出ています。現在のコードでなにか間違い等見受けられたら教えて頂きたいです。