Commit 2abe74fe authored by Thibault Chassat's avatar Thibault Chassat
Browse files

Téléverser un nouveau fichier

parent 7e9665d9
function [y]=f(x)
y=exp(x)./(cos(x).^3+sin(x).^3)
end
function d =fprime(x)
d = (exp(x).* ((cos(x).^3 + sin(x).^3) + 3*sin(x)*cos(x).^2 - 3*sin(x).^2 *cos(x)))/(cos(x).^3 + sin(x).^3).^2
end
//a=linspace(-%pi/5,%pi/2,100);
//b=f(a);
//disp(y)
//clf
//plot(a,f(a))
h=10.^(-20:0.1:-1) //vecteur de valeure de h
x0=%pi/4;
D1 = (f(x0+h)-f(x0))./h; //approximation d'ordre un de la dérivé
e1=D1-fprime(x0) //erreur
D2 = (f(x0+h)-f(x0-h))./(2*h); //approximation d'ordre un de la dérivé
e2=D2-fprime(x0) //erreur
clf
plot(h,abs(e1),h,abs(e2))
legend("ordre1","ordre2")
gca().log_flags="ll"
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment