﻿import matplotlib.pyplot as plt
import numpy as np
import numpy.random as rd
from scipy.integrate import odeint
from math import *
a=1
b=2
y_init=0
t=np.linspace(0,10,1000) # le temps ici va de t=0 à t=10, il y a 1000 points de mesure entre ces deux instants
def F1(y,t):
	return (b-a*y)     
y_sol =odeint(F1,y_init,t)
plt.plot(t,y_sol)
plt.show()
