Problema de Caida Libre Metodo de Integracion00
Problema de Caida Libre Metodo de Integracion00
Problema de Caida Libre Metodo de Integracion00
Se desea saber que tan lejos ha caído un cuerpo después de 10 segundos. conociendo la
distancia exacta d=289.4315 m .Use el método de trapecio,Simpson 1/3 y Simpson 3/8
SOLUCION:
−c
g∗m ∗t
V (t)= ∗(1−e m )
c
−c
dx g∗m ∗t
= ∗(1−e m )
dt c
x 10 −c
∫ dx=∫ g∗m
c
∗t
(
∗ 1−e m ∗ dt )
0 0
10 −c
x=∫
g∗m
c
( ∗t
∗ 1−e m ∗dt )
0
function [h p] = trapezoide(f,a,b,n)
%el nombre de la funcion como string
%n:numero de segmentos
%a:limite inferior
%b:limite superior
%p:integracion
h=(b-a)/n;
n=n+1;
y=zeros(n,1);
x=zeros(n,1);
suma=0;
for i=1:n
x(i)=a+h*(i-1);
y(i)=feval(f,x(i));
end
for i=2:n-1
suma=suma+y(i);
end
p=0.5*h*(y(1)+2*suma+y(n));
end
Y la función:
function f = caidalibre( x )
g=9.8;
m=68.1;
c=12.5;
f= (g*m/c)*(1-exp(-(c/m)*x));
end
Programa principal:
clc
close all
disp('======================================')
disp(' Menu de Integracion ')
disp('======================================')
f=input('ingrese la funcion: ');
a=input('ingrese el limite inferior: ');
b=input('ingrese el limite superior: ');
es=input('ingrese el error fijado:');
d=input('ingrese la distancia verdadera:');
fprintf('Escoge el metodo de Integracion\n');
fprintf('1.Trapecio\n');
fprintf('2.Simpson 1/3\n');
fprintf('3.Simpson 3/8\n');
op=input('Escoger opcion:');
n=1;
et=100;
while (op~=1)&(op~=2)&(op~=3)
op=input(' Escoger nuevamente la opcion 1 , 2 o 3 : ');
end
fprintf('\nSegmentos Tamaño de segmento d[aprox.]
et[porcent.]\n');
while es<et
if op==1
[h p]=trapezoide(f,a,b,n);
et=abs((d-p)/d)*100;
fprintf('\n')
fprintf('%6.0f %20.6f %18.6f %15.6f\n \n',n,h,p,et);
n=n*2;
elseif op==2
[h p]=simpson3(f,a,b,n);
et=abs((d-p)/d)*100;
fprintf('\n')
fprintf('%6.0f %20.6f %18.6f %15.6f\n \n',n,h,p,et);
n=n*2;
elseif op==3
[h p]=simpson8(f,a,b,n);
et=abs((d-p)/d)*100;
fprintf('\n')
fprintf('%6.0f %20.6f %18.6f %15.6f\n \n',n,h,p,et);
n=n*3;
end
end