==========  Corrigé de l'exercice donné le 7 janvier  ====


rm(list=ls(all=TRUE))
epsilon=rnorm(500)
plot(epsilon,type="l")               # y_1
 
x=sin(1:500)
y_6=x+epsilon
plot(y_6,type="l")                   # y_6
 
 
y_5=x*epsilon
plot(y_5,type="l")                   # y_5


s_fichier=read.table("effets_saison.csv")
head(s_fichier)
attach(s_fichier)

s=V1                      

s_part=s[1:60]
epsilon_part=epsilon[1:60]
y_2=s_part+epsilon_part
plot(y_2,type="l")                   # y_2
 
 
 
t=1:60
y_3=1.05*t+s_part+epsilon_part
plot(y_3,type="l")                   # y_3
 
 
tt=1:500
y_4=1.05*tt*s*epsilon
plot(y_4,type="l")                   # y_4