Como puedo graficar la función que ingreso el usuario, por favor
Publicado por JUAN ANGEL (17 intervenciones) el 26/11/2019 07:47:14
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
45
46
47
48
49
50
%Regla de sipson 1/3
clear all;close all; clc
input('Programa de integracion por el metodo de regla de simpson 1/3');
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)/(2*n);
sum=gx(a)+gx(b);
for i=1:n
eval=b+(2*i-1)*h;
sum=sum+4*gx(eval);
end
for i=1:n-1
eval=b+(2*i)*h;
sum=sum+2*gx(eval);
end
sum=(h/3)*sum;
x=sum;
fprintf('Resultado de la integral: %14.14f',rSimpsont(gx,a,b,e,x,n));
function y=rSimpsont(gx,a,b,e,x,n)
n=n+1;
h=(a-b)/(2*n);
aprox=gx(a)+gx(b);
for i=1:n
eval=b+(2*i-1)*h;
aprox=aprox+4*gx(eval);
end
for i=1:n-1
eval=b+(2*i)*h;
aprox=aprox+2*gx(eval);
end
aprox=(h/3)*aprox;
if abs(x-aprox)>e
y=rSimpsont(gx,a,b,e,aprox,n);
else
y=aprox;
fprintf('\nNumero de integrales: %d\n',n);
end
end
Valora esta pregunta
0