Interpolation de Lagrange

Tout ce qui concerne les différents logiciels mathématiques (Matlab, Maxima, Scilab, R, Maple, Mathematica, ...) et les calculatrices scientifiques.
Règles du forum
Merci d'éviter le style SMS dans vos messages et de penser à utiliser la fonction Recherche avant de poster un message. Pour joindre des fichiers à vos messages, consulter ce sujet.
> Penser à utiliser les balises Code pour poster du code.
> Préciser le nom du logiciel concerné dans le titre de votre message.

Interpolation de Lagrange

Messagepar rouxc » Jeudi 27 Novembre 2008, 16:04

Bonjour à tous,

je cherche à écrire l'algorithme d'interpolation lagrangienne sous Maple (version 10) : http://fr.wikipedia.org/wiki/Interpolation_lagrangienne

En fait, j'ai un petit problème d'entrée, pour calculer les Lk. Voici ce que je rentre :

Code: Tout sélectionner
> restart;
> with(linalg):
Warning, the protected names norm and trace have been redefined and unprotected

> Lk:=proc(x,n)local S1,S2,i,k,L;
>   L:=vector(n,0):
>   for k from 1 by 1 to n do
>        S1:=1:
>        S2:=1:
>        print('x[k]'=x[k]);
>        for i from 1 by 1 to n do
>              if (i!=k) then
>                 print('x[i]'=x[i]);
>                 S1:=S1*(x-x[i]):
>                 print('S1'=S1);
>                 S2:=S2*(x[k]-x[i]):
>                 print('S2'=S2);
>              else end if:   
>        end do:
>   end do:
> end proc:
>
> n:=3:
> x:=vector(n,[1,2,3]);
                            x := [1, 2, 3]
> Lk(x,n);
                               x[k] = 1
                               x[i] = 1
                              S1 = x - 1
                                S2 = 0
                               x[k] = 2
                               x[i] = 2
                              S1 = x - 2
                                S2 = 0
                              x[k] = 3


Je ne comprends pas du tout pourquoi Maple ne fais pas de boucle sur i pour chacun des k allant de 1 à n. Du coup je ne peux pas terminer l'algo en faisant le quotient de S1 par S2 (division par 0) ... Quelqu'un peux m'en dire un peu plus ?

Merci d'avance :wink: .
Cordialement.

Clément.
rouxc
Utilisateur
 
Messages: 1
Inscription: Jeudi 27 Novembre 2008, 15:51
Statut actuel: Post-bac | Licence

Publicité

Re: Interpolation de Lagrange.

Messagepar OG » Jeudi 27 Novembre 2008, 16:31

Code: Tout sélectionner
(i<>k)


O.G.
OG
Modérateur
 
Messages: 2286
Inscription: Lundi 12 Mars 2007, 11:20
Localisation: Rouen
Statut actuel: Actif et salarié | Maître de conférence


Retourner vers Logiciels mathématiques

 


  • Articles en relation
    Réponses
    Vus
    Dernier message

Qui est en ligne

Utilisateurs parcourant ce forum: Aucun utilisateur enregistré et 1 invité