ANEXO 1:


clear all
N=20; n = [0:1:N]; h = 10.^-n; % Genero vector de h
x=1; der_f_exacto = cos(x); % Verdadera solución
der_f_aprox = (sin(x+h)-sin(x))./h; % Estimaciones derivada para distintas h's
E_rel = abs(der_f_exacto-der_f_aprox)./abs(der_f_exacto); % Error relativo
E_rel_trunc = h*sin(x)/(2*cos(x));
E_rel_maq = eps*sin(x)./(h*cos(x));
E_rel_cota=E_rel_trunc+E_rel_maq;
subplot(211),semilogx(h,der_f_exacto ,'bo',h,der_f_aprox,'r*'),title('Grafica del valor "exacto" (azul) y aproximado (rojo), respecto de h')
subplot(212),loglog(h,E_rel,'bo'),title('Grafica del Error relativo de la formula (1), respecto de h')


[Cerrar]