python - Python: 4 次ルンゲクッタ法

Question

from math import sin
from numpy import arange
from pylab import plot,xlabel,ylabel,show
def answer():
    print('Part a:')
    print(low(x,t))
    print('First Graph')
    print('')


def low(x,t):
    return 1/RC * (V_in - V_out)

a = 0.0
b = 10.0
N = 1000
h = (b-a)/N
RC = 0.01
V_out = 0.0

tpoints = arange(a,b,h)
xpoints = []
x = 0.0

for t in tpoints:
    xpoints.append(x)
    k1 = h*f(x,t)
    k2 = h*f(x+0.5*k1,t+0.5*h)
    k3 = h*f(x+0.5*k2,t+0.5*h)
    k4 = h*f(x+k3,t+h)
    x += (k1+2*k2+2*k3+k4)/6

plot(tpoints,xpoints)
xlabel("t")
ylabel("x(t)")
show()

だから私はコード化された4次のルンゲクッタ法を持っていますが、私が当てはめようとしている部分は、[2t]が偶数の場合はV_in(t) = 1、[2t]が奇数の場合は-1と言う問題です。

また、次の式を返すかどうかもわかりません: return 1/RC * (V_in - V_out)

問題は次のとおりです。

問題 8.1

あなたが私を助けていただければ幸いです！