Apprendre Maple Index du Forum Apprendre Maple
Site dédié au logiciel de calcul formel Maple
 
  Page d'accueilPage d'accueil   FAQFAQ    RechercherRechercher    Liste des MembresListe des Membres    Groupes d'utilisateursGroupes d'utilisateurs 
S'enregistrerS'enregistrer    ProfilProfil   Se connecter pour vérifier ses messages privésSe connecter pour vérifier ses messages privés   ConnexionConnexion 

bissection

 
Poster un nouveau sujet   Répondre au sujet    Apprendre Maple Index du Forum -> Programmation
Voir le sujet précédent :: Voir le sujet suivant  
Auteur Message
lol



Inscrit le: 05 Nov 2006
Messages: 4

MessagePosté le: 05 Nov 2006 18:15    Sujet du message: bissection Répondre en citant

bonjour ,
j'ai fait une procédure qui approxime le 0 d'une fonction mais il y a un truc bizarre , lorsque je rentre f(x)=x-racine(2) la procédure me renvoit 1 c'est pas normal , en plus je vois pas vraiment l'erreur dans la procédure même : la tête de la procédure:
bisection:=proc(f,a,b,e);

a1:=a;

b1:=b;

while evalf(b1-a1)>=e do

if evalb(f(a1)*f((a1+b1)/2)<0)=true then b1:=evalf(((a1+b1)/2)) ;

elif evalb(f(a1)*f((a1+b1)/2)>0)=true then a1:=evalf(((a1+b1)/2));

else return(x=evalf((a1+b1)/2));

end if ;

end do;

return(x=evalf(((a1+b1)/2)))

end proc;
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé  
ALS



Inscrit le: 11 Sep 2006
Messages: 647

MessagePosté le: 06 Nov 2006 13:10    Sujet du message: Répondre en citant

Il fallait passer une fonction x->f(x) à la procédure.
Exemple: bissection(x->x-sqrt(2),0,2,10^(-6));


Code:


> bissection:=proc(f,a,b,e)
> local a1,b1,m;
> a1:=evalf(a);
> b1:=evalf(b);
> while evalf(abs(b1-a1))>=e do
>    m:=evalf(f((a1+b1)/2)):
>    if evalf(f(a1))*m<0 then b1:=evalf(((a1+b1)/2))
>     elif evalf(f(a1))*m>0 then a1:=evalf(((a1+b1)/2));
>    end if ;
>    print(`x entre `, a1,`  et `, b1);
> end do;
> return(x=evalf((a1+b1)/2))
> end proc:

> bissection(x->x-sqrt(2),0,2,10^(-6));

                   x entre , 1.000000000,   et , 2.


              x entre , 1.000000000,   et , 1.500000000


              x entre , 1.250000000,   et , 1.500000000


              x entre , 1.375000000,   et , 1.500000000


              x entre , 1.375000000,   et , 1.437500000


              x entre , 1.406250000,   et , 1.437500000


              x entre , 1.406250000,   et , 1.421875000


              x entre , 1.414062500,   et , 1.421875000


              x entre , 1.414062500,   et , 1.417968750


              x entre , 1.414062500,   et , 1.416015625


              x entre , 1.414062500,   et , 1.415039062


              x entre , 1.414062500,   et , 1.414550781


              x entre , 1.414062500,   et , 1.414306640


              x entre , 1.414184570,   et , 1.414306640


              x entre , 1.414184570,   et , 1.414245605


              x entre , 1.414184570,   et , 1.414215088


              x entre , 1.414199829,   et , 1.414215088


              x entre , 1.414207458,   et , 1.414215088


              x entre , 1.414211273,   et , 1.414215088


              x entre , 1.414213180,   et , 1.414215088


              x entre , 1.414213180,   et , 1.414214134


                           x = 1.414213657

Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé  
Montrer les messages depuis:   
Poster un nouveau sujet   Répondre au sujet    Apprendre Maple Index du Forum -> Programmation Toutes les heures sont au format GMT + 2 Heures
Page 1 sur 1

 
Sauter vers:  
Vous ne pouvez pas poster de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas voter dans les sondages de ce forum


phpBB

Développé par phpBB © 2001, 2006 phpBB Group
Traduction par : phpBB-fr.com


Apprendre Maple - © 22/06/2025 - Alain Le Stang - Navigation optimisée pour une résolution 1024 x 768.