Bien que CaRMetal ne soit ni un logiciel de calcul numérique ni un logiciel de calcul formel, il possède des macros sur les nombres complexes qui font l’affaire :

Création des macros

Dans une figure CaRMetal, on commence par créer les points d’affixes respectifs $0$, $i$, $-i$ et $1$ dont on aura besoin par la suite.

On va utiliser deux transformations différentes selon la parité de l’indice (M est un point quelconque dans le plan, d’affixe $z$) :

• Si l’indice est pair, c’est $z\mapsto 1+iz$ qu’on matérialise avec la macro "produit de complexes", en sélectionnant les points d’affixes $i$ et $z$, puis la macro "somme de complexes", en sélectionnant le point nouvellement créé (rendu invisible ci-dessous) et le point d’affixe $1$. Le point obtenu est nommé "pair" ci-dessous.
• Si l’indice est impair, c’est $z\mapsto i-iz$ qu’on matérialise avec la macro "produit de complexes", en sélectionnant les points d’affixes $-i$ et $z$, puis la macro "somme de complexes", en sélectionnant le point nouvellement créé (rendu invisible ci-dessous) et le point d’affixe $i$. Le point obtenu est nommé "impair" ci-dessous.

C’est M qu’on peut bouger ci-dessus. On constate que les points d’affixes respectives $0$ et $i$ sont nommés u0 et u1, c’est pour la cohérence avec la suite du TP.

Ensuite (ça ne se voit pas sur la figure, il faut la télécharger pour s’en rendre compte) deux macros ont été créées :

• "pair" ayant pour objet final "pair" et pour objets initiaux P1 (d’affixe $1$), u1 et M
• "impair" ayant pour objet final "impair" et pour objets initiaux u1, P2 (d’affixe $-i$) et M.

Utilisation des macros

Le nombre $n$ tel que l’indice est égal, soit à $2n$, soit à $2n+1$, est le quotient euclidien de l’indice par 2 :

Math.floor(indice/2)

Le reste de la division euclidienne de l’indice par 2 détermine laquelle des macros on applique au point d’indice $n$ :

/*Programme tp 112a
*/
for(indice=2;indice<32;indice=indice+1){
	n=Math.floor(indice/2);
	if(indice%2==0){
		u=ExecuteMacro("u"+indice,"pair","P1,u1,u"+n);
	} else {
		u=ExecuteMacro("u"+indice,"impair","u1,P2,u"+n);
	}
	SetColor(u,"red");
	SetShowName(u,true);
	SetAlias(u,u);
	Println(u+"="+X(u)+"+i*"+Y(u));
}

Ce script est censé être appliqué à la figure dite de base, téléchargeable ci-dessous, où on a supprimé le point M et recadré. Les macros ont donné les alias "pair" et "impair" aux points créés, il a fallu les remplacer par le nom du point u. Comme les points se superposent, leurs coordonnées ont été récupérées par X(u) et Y(u) puis affichées dans une fenêtre texte.

Somme des termes de la suite

La macro "somme de nombres complexes" permet d’additionner les termes de la suite $u$. La modification suivante du script

/*Programme tp 112b
*/
somme="u1";
for(indice=2;indice<32;indice=indice+1){
	n=Math.floor(indice/2);
	if(indice%2==0){
		u=ExecuteMacro("u"+indice,"pair","P1,u1,u"+n);
	} else {
		u=ExecuteMacro("u"+indice,"impair","u1,P2,u"+n);
	}
	somme=ExecuteMacro("s"+indice,"Opérations sur les complexes/Somme","_somme,_u");
	SetColor(somme,"red");
	SetShowName(somme,true);
	SetHide(u,true);
	Println(somme+"="+X(somme)+"+i*"+Y(somme));
}

appliquée à la figure de base, produit la figure suivante :

Pour avoir les valeurs de $S_{100}$ etc. il suffit de modifier la fin de la boucle ci-dessus.

Question B.6

Un test permet de vérifier si oui ou non, $S_{2n}$ est égal à $n+ni$. Mais il est multiple : $n$ doit être pair, X(S) doit être égal à $n$ et Y(S) aussi :

/*Programme tp 112c
*/
somme="u1";
for(indice=2;indice<102;indice=indice+1){
	n=Math.floor(indice/2);
	if(indice%2==0){
		u=ExecuteMacro("u"+indice,"pair","P1,u1,u"+n);
	} else {
		u=ExecuteMacro("u"+indice,"impair","u1,P2,u"+n);
	}
	somme=ExecuteMacro("s"+indice,"Opérations sur les complexes/Somme","_somme,_u");
	SetColor(somme,"red");
	SetShowName(somme,true);
	SetHide(u,true);
	if(indice%2==0){
		if((X(somme)==n)&&(Y(somme)==n)){
			Println(somme+"="+X(somme)+"+i*"+Y(somme));
		}
	}
}

Son exécution sur la figure de base produit le résultat suivant :

qui mène à la recherche de conjectures.