Apprendre Maple

[manonlilou]

Bonjour tout le monde

Je suis toute nouvelle en maple mais je rencontre déjà un souci
En effet, j'ai du programmer une fonction qui regarde si une famille de vecteurs est libre. Si elle est libre je dois la compléter en une base de R^n et si elle n'est pas libre je dois renvoyer 'erreur' . Voici ce que j'ai fais, et le message d'erreur que maple renvoie :


> complete:=proc(L)
> local S;
> convert(vector(L),matrix);
> S:=nrow(L);
> if nrow(L)<> ncol(L)then delrows(L,S..S);
> if det(L)= 0 then print(erreur);
> else base(L);
> fi:
> fi:
> end:
>
> base:=proc(L)
> local V,U,X;
> X:= transpose(L);
> V:=kernel(X);
> U:= L,V;
> end:
>
> v1:=randvector(3);

v1 := [-54, 10, -77]

> v2:=randvector(3);

v2 := [-63, -90, 61]

> v3:=randvector(3);

v3 := [-3, -82, 16]

> L:= [v1,v2,v1+v2];

L := [v1, v2, v1 + v2]

> complete(L);
Error, (in delrows) First argument must be a matrix

Je ne comprend pas pourquoi il me dit cela puisque j'utilise un convert qui devrait justement transformer ma famille de vecteurs en matrice

Je vous remercie par avance pour votre aide

[ALS]

Bonjour,
L était une liste de vecteurs, donc S:=convert(L,matrix); permet la conversion en matrice.
J'utilise ensuite la fonction rank qui calcule le rang de S: si ce rang n'est pas égal au nb de vecteurs de L, la famille est liée: donc le résultat est erreur sinon on utilise base.

Revoyez votre procédure base car la méthode employée en adjoignant un vecteur du noyau ne donnera pas forcément une base de R^n, car certains des premiers vecteurs peuvent être dans le noyau.

[manonlilou]

Ha merci beaucoup, j'avais vu qu'on pouvait utiliser nops mais seulement je n'avais pas compris ce que cette fonction faisait.

Lorsque je rentre ce que vous m'avez indiqué maple me dit :
Error, (in convert/matrix) expecting array or list , pourtant L est bien une liste!
Je suis désolée, à vrai dire j'ai encore beaucoup de mal avec maple

[manonlilou]

Je viens de voir la fonction nullspace() .
Peut-elle être utile pour résoudre mon problème?

[ALS]

Oui, il semble avoir qqs problèmes, je travaille avec une version récente, Maple 15.

J'ai écrit la procédure base: je complète la famille libre avec des vecteurs nuls de façon à obtenir une matrice M de déterminant 0.
Ensuite je remplace les vecteurs nuls ainsi ajoutés par des vecteurs aléatoires jusqu'à obtenir det(M) non nul, ce qui fait sortir de la boucle while.
L'opérateur $ permet de répéter plusieurs fois: exemple x$3 équivaut à x,x,x
seq crée une séquence (ici de vecteurs aléatoires).