1. Définition de fonctions, lecture de données et affichage

1.1. Définition d'une fonction simple, affichage

Nous allons nous intéresser à la fonction {$f(x) = sin(x)+cos(\sqrt{2}*x)$} ; nous souhaitons plus particulièrement trouver ses extrema. La voici sous forme graphique sur l'intervalle {$ [-10,10] $} :

• Définir une fonction de signature

  double sin_x_plus_cos_sqrt2_times_x( double x );
  

  qui calcule et retourne la valeur de {$f$} en {$x$}.

• Définir une fonction test_11() qui appelle la précédente avec 1 puis avec -4.5 comme argument, et qui affiche les résultats retournés.
  
  Appeler cette fonction dans main(), le résultat attendu est :

1.2. Lecture

• Définir une fonction test_12() qui demande une valeur à l'utilisateur, appelle la fonction sin_x_plus_cos_sqrt2_times_x() avec la valeur fournie par l'utilisateur comme argument et affiche le résultat retourné.
  
  Appeler cette fonction dans main(), voici un exemple du résultat attendu (en vert, ce qui est fourni par l'utilisateur) :

2. Structures de contrôle

2.1. Boucle for version 1

• Définir une fonction test_21() qui demande une valeur à l'utilisateur, et affiche le résultat de l'appel de la fonction sin_x_plus_cos_sqrt2_times_x() pour 10 valeurs espacées de 1 en commençant avec la valeur fournie par l'utilisateur.
  
  Vous utiliserez une boucle for avec une variable de boucle commençant à 0.
  
  Appeler cette fonction dans main(), voici un exemple du résultat attendu (en vert, ce qui est fourni par l'utilisateur) :

2.2. Boucle for version 2

• Définir une fonction de signature

void print_sin_x_plus_cos_sqrt2_times_x( double begin, double end, double step );

qui affiche le résultat de la fonction sin_x_plus_cos_sqrt2_times_x() pour toutes les valeurs entre begin et end avec un pas de step.

Vous utiliserez une boucle for avec une variable de boucle commençant à begin.

• Définir une fonction test_22() qui appelle la précédente avec -10, 10 et 2 comme arguments. Appeler cette fonction dans main(), le résultat attendu est :

2.3. Conditionnelle, boucle while

• Définir une fonction test_23() qui demande à l'utilisateur une borne basse, une borne haute (qui devra être supérieure à la borne basse) et un nombre de valeurs à afficher, puis appelle la fonction print_sin_x_plus_cos_sqrt2_times_x() pour afficher les valeurs demandées.
  
  Vous utiliserez une instruction if pour vérifier que la borne haute donnée par l'utilisateur est strictement supérieure à la borne basse, et une instruction while pour refaire une demande en cas de test négatif ; vous ferez de même pour vérifier qu'au moins 2 valeurs sont demandées.
  
  Appeler cette fonction dans main(), voici un exemple du résultat attendu (en vert, ce qui est fourni par l'utilisateur) :

3. Fonction comme argument

3.1. Calcul de dérivée

Pour trouver les extrema de notre fonction, nous allons chercher les points où sa dérivée s'annule. Nous souhaitons que notre recherche d'extrema soit générique et indépendante de la fonction étudiée, ce qui signifie que nous ne connaissons pas sa dérivée formelle. Nous allons donc utiliser un calcul approché de cette dérivée en utilisant la formule {$f'(x) ≃ \frac{f(x + ε) - f(x)}{ε}$}.

Une signature possible de notre fonction, calculant la valeur estimée de la dérivée de la fonction func en x, pourrait être :

double compute_derivative( double func( double ), double x, double epsilon );

Et nous pourrions alors l'utiliser ainsi :

double estimate = compute_derivative( sin_x_plus_cos_sqrt2_times_x, 1, 1e-5 );

Pour recevoir la fonction avec l'argument func, nous allons plutôt utiliser le type std::function de la bibliothèque standard qui est beaucoup plus général, ce qui sera nécessaire pour la suite ; cela ne changera rien à son utilisation.

• Définir une fonction de signature

double compute_derivative( std::function< double( double ) > func, double x, double epsilon );

qui retourne la valeur estimée de la dérivée de la fonction func en x.

• Définir une fonction test_31() qui affiche l'estimation de la dérivée de sin_x_plus_cos_sqrt2_times_x() dans l'intervalle [-4.6, -4.5] par pas de 0.01 avec une valeur d'epsilon égale à 10-5.
  
  Appeler cette fonction dans main(), le résultat attendu est :

3.2. Recherche d'un zéro

• Définir une fonction de signature

  double find_zero( std::function< double( double ) > func, double begin, double end, double precision );
  

  qui retourne un nombre compris dans l'intervalle [begin, end] pour lequel le résultat de l'appel de func sur ce nombre est inférieur en valeur absolue à precision.
  Cette fonction a pour précondition que func( begin ) et func( end ) sont de signes opposés ; si cette précondition n'est pas vérifiée, la fonction retournera NAN.
  Vous procéderez par dichotomie.

• Définir une fonction test_32() qui cherche un zéro de la fonction sin_x_plus_cos_sqrt2_times_x() dans l'intervalle [-2, 0] avec une précision de 10-5.
  
  Appeler cette fonction dans main(), le résultat attendu est :

4. Recherche des extrema

4.1. Recherche des zéros

• Définir une fonction de signature

  int find_all_zeros( std::function< double( double ) > func, double begin, double end, double width, 
                      double precision, double results[], int max_number_of_results );
  

  qui cherche dans chaque intervalle de largeur width (de la forme [begin + n * width, begin + (n + 1) * width] où n est un naturel quelconque), sans dépasser end comme borne haute, un nombre pour lequel le résultat de l'appel de func est inférieur en valeur absolue à precision ; au maximum, max_number_of_results seront retournés dans le tableau results ; la fonction a pour résultat le nombre exact de valeurs présentes dans results.
• Définir une fonction test_41() qui cherche les zéros de la fonction sin_x_plus_cos_sqrt2_times_x() dans l'intervalle [-10, 10], avec une largeur de 0.5, une précision de 10-5 et un maximum de 10 zéros retournés.
  
  Appeler cette fonction dans main(), le résultat attendu est :

4.2. Recherche des extrema

• Définir une fonction de signature

  int find_all_extrema( std::function< double( double ) > func, double begin, double end, double width, 
                        double precision, double epsilon, double results[], int max_number_of_results );
  

  qui cherche dans chaque intervalle de largeur width (de la forme [begin + n * width, begin + (n + 1) * width] où n est un naturel quelconque), sans dépasser end comme borne haute, un extrema de func, un extrema étant défini comme un nombre pour lequel la valeur estimée (calculée en utilisant epsilon) de la dérivée de func en ce nombre est inférieur en valeur absolue à precision ; au maximum, max_number_of_results seront retournés dans le tableau results ; la fonction a pour résultat le nombre exact de valeurs présentes dans results.

• Définir une fonction test_42() qui cherche les extrema de la fonction sin_x_plus_cos_sqrt2_times_x() dans l'intervalle [-10, 10], avec une largeur de 0.5, une précision de 10-5, une dérivée estimée en utilisant 10-5 pour epsilon et un maximum de 10 extrema retournés.
  
  Appeler cette fonction dans main(), le résultat attendu est :