Area debajo de una curva

Area debajo de una curva


T=imput('insertar cantidad de trapecios')
F=imput('insertar la funcion')
M=imput('insertar valor extremo izq')
P=imput('insertar valor extremo der')
H=(b-a)/n
for k=1:T
a=a+h
F(a+h)=eval


Biseccion
function [ xr ] = biseccion( y )
%programa realizado por Julio Echeverri
%la funcion recibe una funcion y que depende de x
%este programa retorna el valor de la raizde una funcion
%con un error aproximado del 0.000001%
%recibe la funcion simbolica dependiente de x

%se hace un analisis, graficando lafuncion, para poder elegir
%bien el intervalo
t=linspace(-50,50,10000);
q=subs(y,t);
plot(t,q)
grid on

%se pide el ingreso de los limites
xl=input('Ingrese el limite inferior: ');
xu=input('Ingrese el limite superior: ');
xr=0;error=100;

if ((subs(y,xu)*subs(y,xl))>0)%sentencia de control para descartar el intervalo
disp(' ')
disp('Intervalo no valido para hallar la raiz')
error=0.0001;
xr=NaN;
end
%se inicia el metodo de biseccion, con un ciclo en este caso while y
%error=0.000001%
while (abs(error)>0.000001)
aux=xr;
xr=(xl+xu)/2;
fxl=subs(y,xl); %se evalua la funcion en XL
fxr=subs(y,xr); %seevalua la funcion en Xr
h=fxl*fxr; %se multiplican los reultados para f(Xr)*f(XL)
if (h<0) %sentencias condicionales para definir los nuevos limites.
xu=xr;
elseif(h>0)
xl=xr;
else
error=0.000001; %en caso de que h=0 se terminan las iteraciones
end
error=((xr-aux)/xr)*100; %se calcula el error relativo aproximado
end
end

biseccion 2clear
format short;
a=input('Introduzca el valor de a: ');
b=input('Introduzca el valor de b: ');
cont=input('Introduzca el número de iteraciones: ');
for k=1:1:cont
c=(a+b)/2;
e=abs((b-a)/2);
A=[k a b c f(c) e];
if f(a)*f(c)<0
b=c;

else
a=c;

end

end
fprintf('\n k a b c f(c) error \n')
fprintf('Solución: c='), disp(c)
fprintf('f(c)='), disp(f(c))
fprintf('error'), disp(e)

fisica
clc
cleart=0:2:100;
%altura=h=60+2.13*t^2-0.0013*t^4+0.000034t^4.751;
h=60+2.13*t.^2-0.0013*t.^4+0.000034*t.^4.751;
x= find(h<0);
N=x(1);

for (K=1:N);
t1(K)=t(K);
h1(K)=h(K);
end
R=[t1,h1];
fprintf ('\n Tiempo Altura\n')
fprintf ('%3.0f %.3f\n',R)
plot (t1,h1),title('trayectoria del cohete'),xlabel('tiempo(seg)'), ylabel('Altura (pies)'),grid


fibonacci
clear
clc
while 1N=input('Indique la cantidad de terminos de la sucesion: ');
if N<3
break
end
fib(1)=1;
fib(2)=1;
k=3;
while k <= N
fib(k)=fib(k-2)+fib(k-1);
k=k+1;
end
fprintf('La secuencia de Fibonacci a %d términos es\n',N);
fprintf('%g ',fib);
fprintf('\n ');
end

biseccion
%metodo de la biseccion.
clear, clc
h=input ('ingrese la funcion a analizar: ');
f=inline(h);
a=input('ingrese limite inferiordel intervalo: ');
b=input('ingrese limite superior del intervalo: ');
tol=input('ingrese la tolerancia deseada: ');

c=0; n=0; MEP =(b-a)/2;
fprintf('\t n \t\ta \t\tc \t\tb \t\tMEP \n')

while(MEP>tol)
c=(a+b)/2;
disp([n,a,c,b,MEP])
if (f(a)*f(c)<0 )
b=c;
else a=c;
end
MEP = (b-a)/2;
n = n+1;
end
fprintf('Raiz encontrada con una tolerancia de%f:\n\t%f\n', tol,c)

practica hecha por mi
clear
clc
A= [2; 2; 4; 9; 30];
B= [43; 32; 54; 132; 45];
C= [32; 64; 922; 43; 45];
Z= [12; 431; 54; 45; 12];
R= [21; 65; 43; 566; 23;];
Q= [A, B, C, Z, R ];
disp ('las notas de los ultimos 5 meses')
disp(Q)
P= [max(max(Q)); min(min(Q))];
disp('la calificacion mas alta y minima son las siguientes')
disp(P)

Ejercicio de quimica
clc
clear
r=0.08205746;M=menu('Seleccione la variable', 'presion', 'volumen', 'No. moles', 'Temp');
if M==1;
n=input('ingrese el No. moles ');
t=input('Ingrese temp ');
v=input('ingrese Vol. ');
p=n*r*t/v;
disp(p)
end
if M==2
n=input('ingrese el No. moles ');
t=input('Ingrese temp ');
p=input('ingrese presion ');

v=n*t*r/p;
disp(v)
end
if M==3
v=input('ingrese volumen...