Error en metodo del trapecio
Publicado por JUAN ANGEL (17 intervenciones) el 10/11/2019 02:13:10
Programo la regla del trapecio y debe terminar con un error mínimo ingresado, lo intento hacer por función re cursiva pero tengo este error. Saben por que?
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
clear
input('Programa de interpolacion por el metodo del Trapecio');
f=input('Ingrese la funcion a evaluar');
a=input('Ingrese el intervalo superior');
b=input('Ingrese el intervalo inferior');
e=input('Ingrese el error minimo');
n=2;
gx=str2func(strcat('@(x) ',f));
h=(a-b)/n;
sum=gx(a)+gx(b);
for i=1:n-1
eval=b+i*h;
sum=sum+2*gx(eval);
end
sum=(h/2)*sum;
x=sum;
%disp(sum);
rTrapecio(f,a,b,e,x,n)
function y=rTrapecio(gx,a,b,e,x,n)
n=n+1;
h=(a-b)/n;
aprox=gx(a)+gx(b);
for i=1:n-1
eval=b+i*h;
aprox=aprox+2*gx(eval);
end
aprox=(h/2)*aprox;
if abs(x-aprox)>e
y=rTrapecio(gx,a,b,e,aprox,n);
else
y=aprox;
end
end
Valora esta pregunta
0