La méthode de Newton pour la langage MATLAB

Réponses à toutes vos questions après le Bac (Fac, Prépa, etc.)
bilalmaths
Messages: 1
Enregistré le: 03 Jan 2014, 18:35

La méthode de Newton pour la langage MATLAB

par bilalmaths » 03 Jan 2014, 18:49

Bonjour tout le monde,
Au dessous mon essai...
[CENTER]J'espère que vous l'aimez :king2: [/CENTER]
clc,clear all
a = 1;
b = 2;
psid = 10^(-8);
n = 0;
f = inline('3*x+sin(x)-exp(x)')
syms x;
f1 = inline('3+cos(x)-exp(x)');
f2 = inline('-sin(x)-exp(x)');
M2 = max(abs(f2(a)),abs(f2(b)));
m1 = min(abs(f1(a)),abs(f1(b)));
tic
if f(a)*f2(a) > 0
xn = a;
xn1 = xn-f(xn)/f1(xn);
psi = (M2/2*m1)*((xn1-xn)^2);
while psi > psid
xn = xn1;
xn1 = xn-f(xn)/f1(xn);
psi = (M2/2*m1)*((xn1-xn)^2);
n = n+1;
end
end
if f(b)*f2(b) > 0
xn = b;
xn1 = xn-f(xn)/f1(xn);
psi = (M2/2*m1)*((xn1-xn)^2);
while psi > psid
xn = xn1;
xn1 = xn-f(xn)/f1(xn);
psi = (M2/2*m1)*((xn1-xn)^2);
n = n+1;
end
end
toc
sprintf('La racine r = %.18f, l''erreur e = %.18f et le nombre d''itération n = %d',xn1,psi,n)
fplot(f,[-0.05,3])
xlabel('x')
ylabel('y')
title('Méthode de Newton : Courbe y = f(x)')
hold on
plot(xn1,f(xn),'*r')
grid on




 

Retourner vers ✯✎ Supérieur

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 67 invités

Tu pars déja ?



Fais toi aider gratuitement sur Maths-forum !

Créé un compte en 1 minute et pose ta question dans le forum ;-)
Inscription gratuite

Identification

Pas encore inscrit ?

Ou identifiez-vous :

Inscription gratuite