\documentclass[10pt]{article}
\usepackage[french]{babel}
\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage[version=4]{mhchem}
\usepackage{stmaryrd}
\usepackage{bbold}
\usepackage{graphicx}
\usepackage[export]{adjustbox}
\graphicspath{ {./images/} }

\begin{document}
\section*{Conception : ESSEC BS \\
 MATHÉMATIQUES 2 APPLIQUÉES \\
 FILIÈRE ÉCONOMIQUE ET COMMERCIALE VOIE GÉNÉRALE}
Vendredi 25 avril 2025, de 14 h. à 18 h.

La présentation, la lisibilité, l'orthographe, la qualité de la rédaction, la clarté et la précision des raisonnements entreront pour une part importante dans l'appréciation des copies.\\
Les candidats sont invités à encadrer dans la mesure du possible les résultats de leurs calculs.\\
Aucun document n'est autorisé. L'utilisation de toute calculatrice et de tout matériel électronique est interdite. Seule l'utilisation d'une règle graduée est autorisée.\\
Si au cours de l'épreuve, un candidat repère ce qui lui semble être une erreur d'énoncé, il la signalera sur sa copie et poursuivra sa composition en expliquant les raisons des initiatives qu'il sera amené à prendre.

Dans ce sujet on s'intéresse à un problème d'arrêt optimal c'est à dire au problème du choix d'un moment pour entreprendre une action spécifique, afin de maximiser un gain attendu ou de minimiser un coût attendu.

Des problèmes d'arrêts optimaux peuvent être trouvés dans les domaines des statistiques, de l'économie et des mathématiques financières (par exemple dans la tarification des options américaines).

On peut modéliser le problème étudié comme suit. Supposons que nous recevions une suite finie de nombres réels, un par un. Ces réels sont des réalisations de variables aléatoires indépendantes de lois connues à l'avance mais pas nécessairement identiques. Nous ne pouvons garder qu'un seul nombre de la suite. À chaque observation, nous pouvons soit sélectionner le nombre actuel, soit pousser notre chance et passer à l'observation suivante. Notre objectif est de maximiser la probabilité de sélectionner le nombre maximal de la suite.

Dans tout le sujet on considère un espace probabilisé \((\Omega, \mathbb{P}, \mathcal{A})\). Toutes les variables aléatoires réelles et événements qui interviennent dans cet énoncé sont définies sur cet espace.

On rappelle que si \(A\) désigne un événement, \(\mathbf{1}_{A}\) est la variable aléatoire qui vaut 1 sur \(A\) et 0 \(\operatorname{sur} \bar{A}\).

La partie 2 utilise des résultats de la partie 1. La partie 3 est indépendante des deux premières et seules les deux dernières questions 19 et 20 de la partie 4 utilisent des résultats établis dans les parties précédentes.

Un aide-mémoire Python se trouve à la fin de l'énoncé.\\
Pour les scripts et fonctions Python, on supposera que les instructions suivantes ont été exécutées :

\begin{verbatim}
import numpy as np,numpy.random as rd,matplotlib.pyplot as plt
\end{verbatim}

\section*{Partie I - Des résultats généraux}
\begin{itemize}
  \item[1.] Montrer que pour tout \(t\) réel, \(\mathrm{e}^{t} \geqslant 1+t\) et que, pour tout \(t>-1, \ln (1+t) \leqslant t\).
  \item[2.] On définit la fonction \(f\) sur \([0,1]\) par \(f(t)=(1+t) \mathrm{e}^{-t}-t\).
  \begin{itemize}
    \item[a)] Étudier les variations de \(f\) et montrer qu'il existe un unique \(\alpha \in] 0,1[\) tel que \(f(\alpha)=0\) et que pour tout \(t \in[0,1], f(t)>0 \Longleftrightarrow t<\alpha\).
    \item[b)] Écrire un programme Python qui renvoie une valeur approchée de \(\alpha\) à \(10^{-3}\) près.
    \item[c)] Montrer que pour tout \(t \in[0,1], \mathrm{e}^{t} \leqslant 1+2 t\).\\
En déduire que \(\alpha \geqslant \frac{1}{\sqrt{2}}\).
  \end{itemize}
  \item[3.] Soit \(U\) une variable aléatoire à densité, à valeurs dans [0,1[, qui suit la loi uniforme et \(p \in] 0, \alpha[\). On pose \(\beta=f(p)\).\\
On définit les variables aléatoires \(X\) et \(Y\) par :
\[
X=1_{[\beta<U \leqslant \beta+p]} \text { et pour tout } \omega \in \Omega, Y(\omega)=\min \left\{k \in \mathbb{N} / U(\omega) \leqslant \sum_{i=0}^{k} \frac{p^{i}}{i!} \mathrm{e}^{-p}\right\}
\]
  \begin{itemize}
    \item[a)] Écrire une fonction Python, minimum(x,p) qui renvoie le minimum de l'ensemble \(\left\{k \in \mathbb{N} / x \leqslant \sum_{i=0}^{k} \frac{p^{i}}{i!} \mathrm{e}^{-p}\right\}\) lorsque \(x \in[0,1[\) et \(p \in] 0,1[\).\\
En déduire une fonction simuly (p) qui réalise une simulation de \(Y\).
    \item[b)] Montrer que \(Y\) suit la loi de Poisson de paramètre \(p\).
    \item[c)] Soit \(k\) un entier, \(k \geqslant 2\). Montrer que si \([Y=k]\) est réalisé alors \([X=0]\) l'est. En déduire que \(\mathbb{P}([X=0] \cap[Y=k])=\frac{p^{k}}{k!} \mathrm{e}^{-p}\).
    \item[d)] Montrer que \([X=0] \cap[Y=1]=\emptyset\). En déduire que
\[
\mathbb{P}([X=1] \cap[Y=1])=p \mathrm{e}^{-p}, \quad \mathbb{P}([X=0] \cap[Y=0])=\beta, \quad \mathbb{P}([X=1] \cap[Y=0])=p\left(1-\mathrm{e}^{-p}\right)
\]
    \item[e)] Montrer que \(\mathbb{P}(X \neq Y)=1+p-(1+2 p) \mathrm{e}^{-p}\), puis que \(\mathbb{P}(X \neq Y) \leqslant 2 p^{2}\).
  \end{itemize}
  \item[4.] Inégalité de Boole - Soit \(k \in \mathbb{N}^{*}, B_{1}, \ldots, B_{k}\) des événements. On pose \(T=\sum_{i=1}^{k} \mathbf{1}_{B_{i}}\).
\end{itemize}

\begin{itemize}
  \item[]
  \begin{itemize}
    \item[a)] Montrer que \(\mathbb{P}(T \geqslant 1)=\mathbb{P}\left(\bigcup_{i=1}^{k} B_{i}\right)\).
    \item[b)] En utilisant une inégalité du cours, que l'on énoncera précisément, en déduire que :
\[
\mathbb{P}\left(\bigcup_{i=1}^{k} B_{i}\right) \leqslant \sum_{i=1}^{k} \mathbb{P}\left(B_{i}\right)
\]
  \end{itemize}
  \item[5.] Soit \(X\) et \(Y\) deux variables aléatoires à valeurs dans \(\mathbb{N}\). On pose
\[
\delta(X, Y)=\sum_{k=0}^{+\infty}|\mathbb{P}(X=k)-\mathbb{P}(Y=k)| \text { et } d(X, Y)=\mathbb{P}(X \neq Y)
\]
  \begin{itemize}
    \item[a)] Justifier que la série définissant \(\delta(X, Y)\) est bien convergente.
    \item[b)] Soit \(k \in \mathbb{N}\) tel que \(\mathbb{P}(X=k) \geqslant \mathbb{P}(Y=k)\).\\
Montrer que \(|\mathbb{P}(X=k)-\mathbb{P}(Y=k)| \leqslant P([X=k] \cap[Y \neq k])\).
    \item[c)] En déduire que pour tout \(k \in \mathbb{N}\),
\[
|\mathbb{P}(X=k)-\mathbb{P}(Y=k)| \leqslant P([X=k] \cap[Y \neq k])+P([X \neq k] \cap[Y=k])
\]
    \item[d)] En conclure que \(\delta(X, Y) \leqslant 2 d(X, Y) \leqslant 2\).
  \end{itemize}
\end{itemize}

\section*{Partie 2 - Une inégalité d'après Hodge et Le Cam}
On conserve les notations de la partie 1.\\
On considère une suite \(\left(U_{k}\right)_{k \in \mathbb{N}^{*}}\) de variables aléatoires à densité indépendantes qui suivent la loi uniforme sur \([0,1[\).

Soit \(n \in \mathbb{N}^{*}\) et \(p_{1}, \ldots, p_{n}\) des réels appartenant à ]0, 1[. Pour tout \(k \in \llbracket 1, n \rrbracket\), on définit, comme \(X\) et \(Y\) dans la question 2 de la partie 1, \(X_{k}\) et \(Y_{k}\) avec \(U_{k}\) et \(p_{k}\) à la place de \(U\) et \(p\).

On pose \(\lambda=\sum_{k=1}^{n} p_{k}, S_{n}=\sum_{k=1}^{n} X_{k}\) et \(T_{n}=\sum_{k=1}^{n} Y_{k}\).\\
On souhaite établir l'inégalité : \(\delta\left(S_{n}, T_{n}\right) \leqslant 4 \sum_{k=1}^{n} p_{k}^{2} \quad(L C)\).

\begin{itemize}
  \item[6.] Montrer que si l'un au moins des \(p_{k}\) est supérieur ou égal à \(\alpha\) alors \((L C)\) est vérifiée.
  \item[-] On suppose dans la suite de cette partie que pour tout \(k \in \llbracket 1, n \rrbracket, p_{k}<\alpha\).
  \item[7.] Justifier brièvement que \(X_{1}, \ldots, X_{n}\) (respectivement \(Y_{1}, \ldots, Y_{n}\) ) sont indépendantes.
  \item[8.] Quelle est la loi de \(T_{n}\) ? Si les \(p_{k}\) sont tous égaux à \(\frac{\lambda}{n}\), quelle est la loi de \(S_{n}\) ? Quelle est alors la limite en loi de la suite \(\left(S_{n}\right)_{n \in \mathbb{N}^{*}}\) ?
  \item[9.] a) Montrer que \(\left[S_{n} \neq T_{n}\right] \subset \bigcup_{k=1}^{n}\left[X_{k} \neq Y_{k}\right]\).
  \begin{itemize}
    \item[b)] En déduire que pour tout \(n \in \mathbb{N}^{*}, \delta\left(S_{n}, T_{n}\right) \leqslant 4 \sum_{k=1}^{n} p_{k}^{2}\).
  \end{itemize}
\end{itemize}

\begin{itemize}
  \item[10.] Un cas particulier - On suppose dans cette question que tous les \(p_{k}\) sont égaux à \(\frac{\lambda}{n}\). Montrer que pour tout \(n \in \mathbb{N}^{*}, \sum_{k=0}^{n}\left|\binom{n}{k}\left(\frac{\lambda}{n}\right)^{k}\left(1-\frac{\lambda}{n}\right)^{n-k}-\frac{\lambda^{k}}{k!} \mathrm{e}^{-\lambda}\right| \leqslant \frac{4 \lambda^{2}}{n}\).
  \item[11.] Une application - Soit \(n\) un entier, \(n \geqslant 2\), on réalise \(n\) expériences de Bernoulli indépendantes avec les probabilités de succès respectives \(\frac{1}{n+1}, \frac{1}{n+2}, \ldots, \frac{1}{2 n}\). On note \(S_{n}\) le nombre de succès de l'expérience totale et \(s_{n}=\sum_{k=1}^{n} \frac{1}{n+k}\).
  \begin{itemize}
    \item[a)] Montrer que pour tout \(k \in \mathbb{N},\left|\mathbb{P}\left(S_{n}=k\right)-\frac{s_{n}^{k}}{k!} \mathrm{e}^{-s_{n}}\right| \leqslant \frac{4}{n}\).
    \item[b)] Établir que pour tout \(k \in \llbracket 1, n \rrbracket\),
\[
\int_{k}^{k+1} \frac{1}{n+t} d t \leqslant \frac{1}{n+k} \leqslant \int_{k-1}^{k} \frac{1}{n+t} d t
\]
    \item[c)] En déduire un encadrement de \(s_{n}\) puis que \(\lim _{n \rightarrow+\infty} s_{n}=\ln (2)\).
    \item[d)] En conclure que \(\left(S_{n}\right)_{n \geqslant 2}\) converge en loi vers une variable aléatoire \(S\) qui suit la loi de Poisson de paramètre \(\ln (2)\).
  \end{itemize}
\end{itemize}

\section*{Partie 3 - Étude du maximum d'une fonction}
Pour tout \(x \in\left[0,+\infty\left[\right.\right.\), on pose \(h(x)=\left\{\begin{array}{ll}\frac{\mathrm{e}^{x}-1}{x} & \text { si } x>0 \\ 1 & \text { si } x=0\end{array}\right.\).

\begin{itemize}
  \item[12.] Montrer que \(h\) est de classe \(C^{1}\) sur \(\mathbb{R}^{+}\)et préciser la valeur de \(h^{\prime}(0)\).
  \item[-] On définit alors la fonction \(g\) sur \(\mathbb{R}^{+}\)par \(g(x)=\mathrm{e}^{-x} \int_{0}^{x} h(t) d t\).
  \item[13.] a) Montrer que pour tout \(x \geqslant 0, g^{\prime}(x)=\mathrm{e}^{-x}\left(1-\int_{0}^{x}\left(h(t)-h^{\prime}(t)\right) d t\right)\).
  \begin{itemize}
    \item[b)] Montrer que pour tout \(t>0, h(t)-h^{\prime}(t)=\frac{\mathrm{e}^{t}-1-t}{t^{2}}\).
    \item[c)] En déduire que pour tout \(t \geqslant 0, h(t)-h^{\prime}(t)>0\) et \(\lim _{x \rightarrow+\infty} \int_{0}^{x}\left(h(t)-h^{\prime}(t)\right) d t=+\infty\).
    \item[d)] Dans un même tableau, représenter le signe de \(g^{\prime}\) et les variations de \(g\) en justifiant que \(g\) possède un maximum qui est atteint en un unique réel noté \(\gamma>0\) que l'on fera apparaitre dans ce tableau.
  \end{itemize}
  \item[14.] a) Montrer que pour tout \(t>0, \frac{1}{2} \leqslant \frac{\mathrm{e}^{t}-1-t}{t^{2}} \leqslant \frac{1}{2} \mathrm{e}^{t}\).
  \begin{itemize}
    \item[b)] En déduire que pour tout \(x \geqslant 0, \mathrm{e}^{-x} \frac{3-\mathrm{e}^{x}}{2} \leqslant g^{\prime}(x) \leqslant \mathrm{e}^{-x}\left(1-\frac{x}{2}\right)\).
    \item[c)] En conclure que \(\gamma \in[\ln (3), 2]\).
  \end{itemize}
  \item[15.] a) Soit \(n \in \mathbb{N}^{*}\) et \(x \in[0, n]\), montrer que :
\[
\sum_{k=0}^{n} \frac{x^{k}}{k!} \leqslant \mathrm{e}^{x} \leqslant \sum_{k=0}^{n} \frac{x^{k}}{k!}+\frac{x^{n}}{n!} \sum_{k=n+1}^{+\infty}\left(\frac{x}{n+1}\right)^{k-n} \leqslant \sum_{k=0}^{n} \frac{x^{k}}{k!}+\frac{x^{n+1}}{n!}
\]
et en déduire que \(\sum_{k=1}^{n} \frac{x^{k-1}}{k!} \leqslant h(x) \leqslant \sum_{k=1}^{n} \frac{x^{k-1}}{k!}+\frac{x^{n}}{n!}\).
\end{itemize}

\begin{itemize}
  \item[]
  \begin{itemize}
    \item[b)] En conclure que pour tout \(n \in \mathbb{N}^{*}\) et \(x \in[0, n]\),
\[
\mathrm{e}^{-x}\left(\sum_{k=1}^{n} \frac{x^{k}}{k!k}\right) \leqslant g(x) \leqslant \mathrm{e}^{-x}\left(\sum_{k=1}^{n} \frac{x^{k}}{k!k}\right)+\mathrm{e}^{-x} \frac{x^{n+1}}{(n+1)!}
\]
puis que pour tout \(x \geqslant 0, g(x)=\mathrm{e}^{-x} \sum_{k=1}^{+\infty} \frac{x^{k}}{k!k}\).
  \end{itemize}
  \item[16.] Compléter le programme suivant pour qu'il trace la partie de la courbe de \(g\) comprise entre les abscisses \(\ln (3)\) et 2 , les valeurs de \(g\) étant calculées à \(10^{-4}\) près :
\begin{verbatim}
X = np.linspace(np.log(3),2,100)
Y = []
for x in X:
    n = 2;s = x + x**2/4;d = x**3/6
    while d * np.exp(-x) > 0.0001:
        n = n + 1
        s = s + ...
        d = d*x/...
    Y.append(s*...)
plt.plot(X,Y)
plt.grid()
plt.show()
\end{verbatim}

On obtient le graphique suivant :\\
\includegraphics[max width=\textwidth, alt={}, center]{4ce37f55-f01f-4a29-b9a4-e14a084f677d-05_672_937_1428_579}
\end{itemize}

\section*{Partie 4 - Le problème du meilleur choix}
On considère \(\left(X_{k}\right)_{k \in \mathbb{N}^{*}}\) une suite de variables aléatoires indépendantes à densité.\\
Soit \(s \in \mathbb{R}\), pour tout \(n \in \mathbb{N}^{*}\) on définit des variables aléatoires, \(a_{n, s}, Y_{n, s}, Z_{n}\) et \(K_{n, s}\) par, pour tout \(\omega \in \Omega\) :

\[
a_{n, s}(\omega)=\min \left(\left\{k \in \llbracket 1, n \rrbracket / X_{k}(\omega)>s\right\} \cup\{n\}\right), Y_{n, s}(\omega)=X_{a_{n, s}(\omega)}(\omega), Z_{n}(\omega)=\max _{k \in \llbracket 1, n \rrbracket}\left(X_{k}(\omega)\right)
\]

et \(K_{n, s}(\omega)\) est égal au nombre d'indices \(k \in \llbracket 1, n \rrbracket\) tels que \(X_{k}(\omega)>s\).\\
On cherche à choisir \(s\) pour maximiser \(r_{n}=\mathbb{P}\left(Y_{n, s}=Z_{n}\right)\).\\
On pose pour tout \(k \in \llbracket 1, n \rrbracket, p_{k}=\mathbb{P}\left(X_{k}>s\right)\) et on suppose que \(p_{k} \neq 1\).

\begin{itemize}
  \item[17.] Une minoration dans le cas général
  \begin{itemize}
    \item[a)] Montrer que \(\mathbb{P}\left(Y_{n, s}=Z_{n}\right) \geqslant \mathbb{P}\left(K_{n, s}=1\right)\).
    \item[b)] On pose \(\theta=\mathbb{P}\left(Z_{n} \leqslant s\right)\). Montrer que \(\mathbb{P}\left(K_{n, s}=1\right)=\theta \sum_{k=1}^{n} \frac{p_{k}}{1-p_{k}}\).
    \item[c)] En déduire que \(\mathbb{P}\left(Y_{n, s}=Z_{n}\right) \geqslant-\theta \ln (\theta)\).
    \item[d)] En déduire l'existence d'au moins une valeur de \(s\), que l'on définira à l'aide de la fonction de répartition \(F_{n}\) de \(Z_{n}\), pour laquelle \(\mathbb{P}\left(Y_{n, s}=Z_{n}\right) \geqslant \frac{1}{\mathrm{e}}\).
  \end{itemize}
  \item[-] On suppose désormais que les \(X_{k}\) suivent la même loi donc que les \(p_{k}\) sont tous égaux, non nuls.\\
On note \(p\) cette valeur commune, \(F\) la fonction de répartition et \(f\) une densité communes aux \(X_{k}\).\\
On admet que si \(X\) et \(Y\) sont deux variables à densité indépendantes, de fonction de répartition \(F_{X}\) pour \(X\) et de densité \(f_{Y}\) pour \(Y\), on a alors \(\mathbb{P}(X \leqslant Y)=\int_{-\infty}^{+\infty} F_{X}(t) f_{Y}(t) d t\).\\
On rappelle que si \(A\) est un événement de probabilité non nulle, \(\left(\Omega, \mathcal{A}, \mathbb{P}_{A}\right)\) est un espace probabilisé admettant les mêmes variables aléatoires et ayant les mêmes propriétés que \((\Omega, \mathcal{A}, \mathbb{P})\).
  \item[18.] Une estimation - On suppose dans cette question les \(X_{k}\) suivent la loi uniforme sur \([0,1[\).
  \begin{itemize}
    \item[a)] Montrer que \(s=1-p\).
    \item[b)] Écrire une fonction Python simulCouple (n,p) qui renvoie une simulation du couple \(\left(Z_{n}, Y_{n, s}\right)\).
    \item[c)] Écrire un programme Python qui réalise et affiche une estimation de \(r_{10}\) dans ces conditions lorsque \(n=10\) et \(p=0.15\).
  \end{itemize}
  \item[19.] Une expression explicite de \(r_{n}\) - Soit \(k \in \llbracket 1, n \rrbracket\), on note \(I_{1}, \ldots, I_{\binom{n}{k}}\) les parties à \(k\) éléments de \(\llbracket 1, n \rrbracket\).\\
Pour tout \(j \in \llbracket 1,\binom{n}{k} \rrbracket\), on définit \(A_{j}\) l'événement \(\left(\bigcap_{i \in I_{j}}\left[X_{i}>s\right]\right) \bigcap\left(\bigcap_{i \notin I_{j}}\left[X_{i} \leqslant s\right]\right)\).
  \begin{itemize}
    \item[a)] Montrer que pour tout \(x\) réel et \(i \in I_{j}\),
\[
\mathbb{P}_{A_{j}}\left(X_{i} \leqslant x\right)= \begin{cases}\frac{F(x)-F(s)}{p} & \text { si } x>s \\ 0 & \text { sinon }\end{cases}
\]
et que les \(X_{i}\) pour \(i \in I_{j}\) sont indépendantes pour la probabilité \(\mathbb{P}_{A_{j}}\).
    \item[b)] En déduire que pour tout \(r \in I_{j}\) :
\[
\mathbb{P}_{A_{j}}\left(X_{r}=\max _{i \in I_{j}}\left(X_{i}\right)\right)=\frac{1}{p^{k}} \int_{s}^{+\infty}(F(t)-F(s))^{k-1} f(t) d t=\frac{1}{k}
\]
  \end{itemize}
\end{itemize}

\begin{itemize}
  \item[]
  \begin{itemize}
    \item[c)] Montrer que :
\[
\mathbb{P}\left(\left[Y_{n, s}=Z_{n}\right] \cap\left[K_{n, s}=k\right]\right)=\sum_{j=1}^{\binom{n}{k}} \mathbb{P}\left(\left[Y_{n, s}=Z_{n}\right] \cap A_{j}\right) \text { et que } \mathbb{P}_{A_{j}}\left(Y_{n, s}=Z_{n}\right)=\frac{1}{k}
\]
    \item[d)] Montrer de même que \(\mathbb{P}\left(\left[Y_{n, s}=Z_{n}\right] \cap\left[K_{n, s}=0\right]\right)=\frac{1}{n}(1-p)^{n}\).
    \item[e)] En conclure que \(r_{n}=\frac{1}{n}(1-p)^{n}+\sum_{k=1}^{n} \frac{1}{k}\binom{n}{k} p^{k}(1-p)^{n-k}\).
  \end{itemize}
  \item[20.] Comportement asymptotique - On suppose que \(p=\frac{\lambda}{n}, \lambda\) étant un réel strictement positif ne dépendant pas de \(n\).
  \begin{itemize}
    \item[a)] En utilisant un résultat de la partie 2, montrer que pour tout \(n \in \mathbb{N}^{*}\),
\[
\sum_{k=1}^{n} \frac{1}{k}\left|\binom{n}{k}\left(\frac{\lambda}{n}\right)^{k}\left(1-\frac{\lambda}{n}\right)^{n-k}-\frac{\lambda^{k}}{k!} \mathrm{e}^{-\lambda}\right| \leqslant \frac{4 \lambda^{2}}{n}
\]
    \item[b)] En déduire que \(\lim _{n \rightarrow+\infty} r_{n}=\left(\sum_{k=1}^{+\infty} \frac{\lambda^{k}}{k!k}\right) \mathrm{e}^{-\lambda}\).
  \end{itemize}
  \item[21.] On suppose que \(n\) est assez grand pour pouvoir considérer que \(r_{n}\) vaut \(\left(\sum_{k=1}^{+\infty} \frac{\lambda^{k}}{k!k}\right) \mathrm{e}^{-\lambda}\). Comment choisir \(s\) pour que cette probabilité soit maximale ?
\end{itemize}

\section*{Aide-mémoire Python}
Toutes les fonctions et instructions présentées ne sont pas utiles et il est possible d'utiliser d'autres fonctions ou instructions absentes de cet aide-mémoire.

Listes\\[0pt]
[] Créer une liste vide\\[0pt]
[a]\textit{n ou n}[a] Créer une liste avec \(n\) fois l'élément a\\
L. append (a) Ajoute l'élément a à la fin de la liste L\\
L1 + L2 Concatène les deux listes L1 et L2\\
L.count (a) Renvoie le nombre d'occurences de a dans la liste L\\
max(L) Renvoie le plus grand élément de la liste L\\
min(L) Renvoie le plus petit élément de la liste L\\
sum (L) Renvoie la somme de tous les éléments de la liste L\\
a in L Vaut True si a se trouve au moins une fois dans L et False sinon\\
Module mathématique numpy

\begin{verbatim}
import numpy as np
np.linspace(a,b,n) Crée un vecteur de n valeurs uniformément réparties entre a et b (inclus)
    np.zeros(n) Crée le vecteur nul de taille n
    np.ones(n) Crée le vecteur de taille (n dont tous les coefficients valent 1
    np.max(M) Renvoie le plus grand élément de M, matrice ou vecteur
    np.min(M) Renvoie le plus petit élément de M, matrice ou vecteur
np.arange(a,b,eps) Renvoie le vecteur des flottants de a à b de pas constant eps, b étant exclu
    np.sqrt(x) Renvoie 5x, 5i x 
    np.log(x) Renvoie ln(x) si x > 0
    np.exp(x) Renvoie ex
        np.e Renvoie e
\end{verbatim}

Sous module random de numpy pour la simulation probabiliste

\begin{verbatim}
import numpy.random as rd
        rd.random([r,s]) Simule une réalisation d'une matrice (r,s) dont les coefficients sont des variables
            aléatoires indépendantes qui suivent la loi uniforme U([0,1])
    rd.binomial(n,p,[r,s]) Simule une réalisation d'une matrice (r,s) dont les coefficients sont des variables
            aléatoires indépendantes qui suivent la loi binomiale B(n, p)
        rd.poisson(a,[r,s]) Simule une réalisation d'une matrice (r,s) dont les coefficients sont des variables
            aléatoires indépendantes qui suivent la loi de Poisson P(a)
\end{verbatim}

Si le paramètre [r,s] est remplacé par r, ces fonctions renvoient la réalisation d'un vecteur de longueur r correspondant à la loi en question, et si ce paramètre est omis, elles renvoient un seul coefficient suivant les mêmes contraintes.

Sous module graphique pyplot de matplolib

\begin{verbatim}
import matplotlib.pyplot as plt
    plt.plot(X,Y,options) Crée la courbe des points définis par les listes X, abscisses, et Y, ordonnées
            suivant les options graphiques définies par la chaîne de caractères
            facultative options
    plt.xlim(xmin,xmax) Fixe les bornes de l'axe des abscisses
    plt.ylim(ymin,ymax) Fixe les bornes de l'axe des ordonnées
        plt.show() Affiche le graphique
        plt.grid() Affiche un quadrillage
\end{verbatim}

FIN DE L'ÉNONCÉ


\end{document}