Integracion con Matlab



























































TRABAJO ENCARGADO
1. Codigo Matlab:

clc
close all
clear all

format longpromptdd={'Ingrese su funcion'};
namedd='Ingreso de la funcion';
numlinesdd=1;
defaultanswer1dd={'0'};
options.Resize='on';respuesta1dd=inputdlg(promptdd,namedd,numlinesdd,defaultanswer1dd,options);

if isempty(respuesta1dd) %Si se aprieta el boton cancelar se sale de la rutina

return;

else

end%888888888888888888888888888888888
%METODO TRAPEZOIDAL

prompt={'Ingrese el limite inferior:','Ingrese el limite superior','numero de intervalos N:'};
name='Ingrese valores del grafico y el número';
numlines=1;defaultanswer={'0','0','0'};
options.Resize='on';

respuesta=inputdlg(prompt,name,numlines,defaultanswer,options);

if isempty(respuesta) %Si se aprieta el boton cancelar se sale de la rutinareturn;
else
end

syms x;

iniciodd=respuesta1dd

yu=char(iniciodd)

%derfunddd=diff(yu,x)
% iniciodd1 = char(iniciodd)

f = inline(char(iniciodd));
% f1 = inline('0*x');
%f2 =inline(char(derfunddd));

Valoreslimite = str2double(respuesta);

a1=Valoreslimite(1,1);
b1=Valoreslimite(2,1);
N=Valoreslimite(3,1);
% myfun1(x,[a1 b1]);
h=(b1-a1)/N
% Otraforma de generar el gráfico
x=[a1:h:b1]

x=x';

fa=f(a1)
fb=f(b1)

cont=0

for i=2:N
dd=f(x(i))
cont=cont+dd
endAreaTotal=(h/2)*(fa+fb)+h*cont


I(1,1)=AreaTotal

xlswrite('AREA',I,'Integracion','A2');

% Mostrando resultados
IteracTrapezoide=1
Areatotal=I(1,1)







%----------------------------------------------------





%888888888888888888888888888888888

%METODO SIMPSONS3


prompt={'Ingrese el limite inferior:','Ingrese el limite superior','numero...