cat("TP1: \n")

cat("Ex 1: \n")

x=c(1,3,5,7,9)
y=c(2,3,5,7,11,13)
z=c(9,3,2,5,9,2,3,9,1)

cat("x+2:",x+2,"\n")
cat("y*3:",y*3,"\n")
cat("length(x):",length(x),"\n")
cat("sum(x>5):",sum(x>5),"\n")
cat("sum(x[x>5]):",sum(x[x>5]),"\n")
cat("sum(x>5|x<3):",sum(x>5|x<3),"\n")
cat("y[3]:",y[3],"\n")
cat("y[-3]:",y[-3],"\n")
cat("sum(y>3&y<11):",sum(y>3&y<11),"\n")
cat("y[x]:",y[x],"\n")
cat("(y>7):",(y>7),"\n")
cat("y[y>7]:",y[y>7],"\n")
cat("sort(z):",sort(z),"\n")
cat("sort(z,dec=TRUE):",sort(z,dec=TRUE),"\n")
cat("rev(z):",rev(z),"\n")
cat("order(z):",order(z),"\n")
cat("unique(z):",unique(z),"\n")
cat("duplicated(z):",duplicated(z),"\n")
cat("table(z):",table(z),"\n")
cat("rep(z,3):",rep(z,3),"\n")

cat("\n Ex 2: \n")

v0=c(51.4,64,71,82,98.7,pi)
v1=rep(c(5,6,3),5)
v2=rep(1:10,1:10)
v3=c("Killian","Nicolas","Pierre")
v4=c(1i,8+5i,9-7.5i)

cat("v0:",v0,"\n")
cat("v1:",v1,"\n")
cat("v2:",v2,"\n")
cat("v3:",v3,"\n")
cat("v4:",v4,"\n")

cat("\n Ex 4: \n")

x=c(2:6)
cat("Mode de x:",mode(x),"\n")
a=c(x[1],x[2],x[3])
cat("a:",a,"\n")
b=c(x[1],x[3],x[5])
cat("b:",b,"\n")
cat("a+b:",a+b,"\n")
cat("x+a:",x+a,"\n") #R recycle
cat("\nx*2:",x*2,"\n")
cat("a*b:",a*b,"\n") #Multiplication donnée à donnée
cat("ind<-which(x%%2==0):",ind<-which(x%%2==0),"\n") #Renvoie le vecteur non composé des éléments divisibles par 2
cat("x[which(x%%2==0&x%%3==0)]:",x[which(x%%2==0&x%%3==0)],"\n") #Renvoie le vecteur composé des éléments divisibles par 2 et par 3
cat("x[which(x%%2==0|x%%3==0)]:",x[which(x%%2==0|x%%3==0)],"\n") #Renvoie le vecteur composé des éléments divisibles par 2 ou par 3
cat("sum(x):",sum(x),"\n")
x0=c(NA,x)
cat("x0:",x0,"\n")
cat("sum(x0):",sum(x0),"\n")
cat("ls()",ls(),"\n")

cat("\n Ex 5:\n")

y=rep(0:10,2)
cat("matrix(y,nrow=2):\n")
matrix(y,nrow=2)
cat("matrix(y,nrow=2,ncol=11):\n")
matrix(y,nrow=2,ncol=11)
cat("matrix(y,nrow=2,byrow=TRUE):\n")
matrix(y,nrow=2,byrow=TRUE)

A=matrix(seq(0,20,2),nrow=4,ncol=3,byrow=TRUE)
cat("A:\n")
print(A)
cat("\n")

B=matrix(1:12,nrow=4,ncol=3)
cat("B:\n")
print(B)
cat("\n")

cat("1ère colonne de A:",A[,1],"\n")
cat("2e ligne de A:",A[2,],"\n")

C=rbind(A[2,],A[4,])
cat("C:\n")
print(C)
cat("\n")

cat("colSums(A):",colSums(A),"\n")
cat("rowSums(A):",rowSums(A),"\n")

cat("A*B:\n")
print(A*B)
cat("\n")

#cat("A%*%B:\n")
#A%*%B
#cat("\n")

cat(" Ex 8:\n")

nl=c("Marrons","Noisettes","Bleus")
nc=c("Bruns","Chatains","Blonds")
v=c(68,15,25,119,54,113,33,24,141)
X=matrix(v,nrow=3,ncol=3)
dimnames(X)=list(nl,nc)
cat("X:\n")
print(X)
cat("\n")

cat("\n Ex 9:\n")

A=matrix(0,nrow=5,ncol=5)
B=abs(col(A)-row(A))+1

cat("A:\n")
print(A)
cat("\n")

cat("B:\n")
print(B)
cat("\n")

cat("det(B):",det(B),"\n")

cat("Solution du système:",solve(B,c(1,2,2,3,2)),"\n\n")

cat("TP2: \n")

cat("\nEx 1: \n")

x=seq(-3*pi,3*pi,0.1)
y=sin(x)
#plot(x,y)
#plot(x,y,type="l",col="red")
#plot(x,y,xlab="Titre - x",ylab="Titre - y",main="Titre")
#plot(x,y,type="p",col="black")
#plot(x,y,type="b",col="tomato")
#plot(x,y,type="o",col="aquamarine",pch=1)

cat("\nEx 2: \n")

x=seq(-4*pi,4*pi,0.1)
y=cos(x)
#plot(x,y,main="Graphe de la fonction cosinus",col="green")

cat("\nEx 3: \n")

cat("3 parmi 10:",choose(10,3),"\n")

cat("\nEx 5: \n")

x=seq(-5,5,0.1)
#plot(x,dnorm(x),type="l",main="Densité de la loi normale centrée réduite")
#curve(expr=dnorm,xlim=c(-3,3),main="Densité de la loi normale centrée réduite")
#segments(x0=c(-4,1),y0=dnorm(1),x1=1,c(dnorm(1),-1),lty=2,col="blue")
#curve(expr=pnorm,xlim=c(-3,3),add=TRUE)

cat("\nEx 7: \n")

#curve(expr=dpois(x,3),xlim=c(-3,3),main="Densité de la loi de Poisson de paramètre 3")
#curve(expr=ppois(x,3),xlim=c(-3,3),main="Fonction de répartition de la loi de Poisson de paramètre 3")

#curve(expr=dexp,xlim=c(-3,3),main="Densité de la loi exponentielle de paramètre 1")
#curve(expr=pexp,xlim=c(-3,3),main="Fonction de répartition de la loi exponentielle de paramètre 1")

#curve(expr=dchisq(1,6),xlim=c(-3,3),main="Densité de la loi du chi-deux en paramètre 1 à 6") #FAUX
#curve(expr=pchisq(1,6),xlim=c(-3,3),main="Fonction de répartition  de la loi du chi-deux en paramètre 1 à 6") #FAUX