import matplotlib.pyplot as plt
import numpy as np

N=1000
tf=1
h=tf/N

omega0=10
omega=100
E=1

def f(x,t):
    return ......................

def Euler(N,tf,h,f):
    X=[1]
    T=[0]
    for i in range(N-1):
        T.append(.............)
        X.append(...............)
    return X,T
    
X,T=Euler(N,tf,h,f)

Entree=[]
for i in range(N):
    Entree.append(E*np.cos(omega*T[i]))

plt.plot(T,X,'r+-',label='s(t)')
plt.plot(T,Entree,'b-',label='e(t)')
plt.xlabel("t(s)")
plt.ylabel("x")
plt.legend()
plt.show()