使用np.linspace定義網格並評估該網格上的函數

Question

我定義了以下時間網格：

T = 10000
tmin = 0
tmax = 10**20
t = np.linspace(tmin, tmax, T)

所以我的步長約為10 ** 16。 我想根據該步長評估從tmin到tmax的函數。 以下是對tmin到10,000的評估還是對整個網格的評估？

timedep_PD = np.zeros([t,Y])
for t in range(tmin,T):
    timedep_PD[t]= probdensity_func(x,t,0)

我如何以（tmax-tmin / T）的步長評估從tminx到tmax此函數，同時仍以timedep_PD[t]的形式存儲

Answer 1

為了簡單起見：

T = 10
tmin = 0
tmax = 10**2
t = np.linspace(tmin, tmax, T)

timedep_PD = np.zeros((len(t),T)) ##len(t) is T : np.zeros((T,T))
dt = ((tmax-tmin)/T)
k=0
for t in np.arange(tmin,tmax,dt):
    timedep_PD[k]= np.sin(t)
    print timedep_PD[k]
    k+=1

我們也可以使用np.meshgrid