serie de taylor para cualquier funcion
Publicado por pous (4 intervenciones) el 11/06/2013 22:42:30
Estoy intentando escribir un algoritmo que resuelva la serie de Taylor para cualquier funcion. Sin embargo hasta el momento aun no me funciona. agradecere cualquier ayuda.
Lo que llevo hasta el momento es:
clear all
close all
clc
format longG
syms x
f = inline(input('f(x): ','s'));
x1 = input('X(i): ');
x2 = input('X(i+1): ');
n = input('Número de Orden: ');
h = x2 - x1;
v1 = f(x1);
v2 = f(x2);
et = abs((v2 - v1)/v2 * 100);
X = [0 v1 0 et];
for i = 1:n
d = diff(sym(f(x1),i));
v1 = v1 + (d*((h^i)/factorial(i)));
et = abs(((v2 - v1)/v2) * 100);
X(i+1,1) = i;
X(i+1,2) = v1;
ea = abs((X(i+1,2)-X(i,2))/X(i,2)*100);
X(i+1,3) = ea;
X(i+1,4) = et;
if ea <= 0
break
end
end
Lo que llevo hasta el momento es:
clear all
close all
clc
format longG
syms x
f = inline(input('f(x): ','s'));
x1 = input('X(i): ');
x2 = input('X(i+1): ');
n = input('Número de Orden: ');
h = x2 - x1;
v1 = f(x1);
v2 = f(x2);
et = abs((v2 - v1)/v2 * 100);
X = [0 v1 0 et];
for i = 1:n
d = diff(sym(f(x1),i));
v1 = v1 + (d*((h^i)/factorial(i)));
et = abs(((v2 - v1)/v2) * 100);
X(i+1,1) = i;
X(i+1,2) = v1;
ea = abs((X(i+1,2)-X(i,2))/X(i,2)*100);
X(i+1,3) = ea;
X(i+1,4) = et;
if ea <= 0
break
end
Valora esta pregunta
0