Práctica 12-Integración Numérica MatLab

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 4

UNIVERSIDAD CENTRAL DE VENEZUELA.

FACULTAD DE INGENIERÍA-DEPARTAMENTO DE MATEMÁTICAS APLICADAS –CÁLCULO NUMÉRICO.


TEMA 5: DIFERENCIACIÓN E INTEGRACIÓN NUMÉRICA

INTEGRACIÓN NUMÉRICA UTILIZANDO MATLAB.


Introducción:
A continuación se presenta una sintaxis en Matlab para aplicar varios métodos de integración
numérica

Pasos para aproximar integrales definidas en Matlab.

PASO 1: Abra el script integración.m El mismo contiene el código en Matlab para aplicar
varios métodos de integración numérica. Este código es:

clear; format long


fprintf('\n\t::::::::::PROGRAMADO POR: GERARDO DOS RAMOS
::::::::::\n');
while 1
disp('[1] TRAPECIO SIMPLE (h>0)')
disp('[2] TRAPECIO COMPUESTO(h<0)')
disp('[3] FORMULA DE SIMPSON SIMPLE')
disp('[4] FORMULA DE LOS TRES OCTAVOS DE SIMPSON')
disp('[5] FORMULA DE SIMPSON COMPUESTO')
disp('[6] INTEGRACION DE ROMBERG')
disp('[7] ROMBERG MODIFICADO')
disp('[8] VOLVER')
elecc3 = input('ELIGA OPCION ');
switch elecc3
case 1
clc; clear;
fprintf('\t\tTRAPECIO SIMPLE\n')
funcion=input('ingrese la funcion \n f(x)=','s');
b=input('ingrese el limite superior de la integral\n');
a=input('ingrese el limite inferior de la integral\n');
h=b-a;
x=a;
f=eval(funcion); x=b;
f= (f+eval(funcion))*(h/2);
fprintf('El valor aproximado es: %10.15f\n\n',f)
case 2
clc; clear;
fprintf('\t\tTRAPECIO COMPUESTO\n')
funcion=input('ingrese la funcion \n f(x)=','s');
b=input('ingrese el limite superior de la integral\n');
a=input('ingrese el limite inferior de la integral\n');
n=input('ingrese el numero de intervalos\n');
h=(b-a)/n;
f=0;
for k=1:n-1
x=a+h*k;
f=f+eval(funcion);
end
f=2*f;
x=a; f=f+eval(funcion); x=b; f=f+eval(funcion);
f=(h/2)*(f);
fprintf('El valor aproximado es: %10.15f\n\n',f)
case 3
clc; clear;
fprintf('\t\tFORMULA DE SIMPSON SIMPLE\n')
funcion=input('ingrese la funcion \n f(x)=','s');
b=input('ingrese el limite superior de la funcion\n');
a=input('ingrese el limite inferior de la integral\n');
h=(b-a)/2;
x=a; f=eval(funcion); x=b; f=f+eval(funcion);
x=a+h; f=f+ 4*(eval(funcion));
f=(h/3)*f;
fprintf('El valor aproximado de la integral es: %10.15f\n\n',f)
case 4
clc; clear;
fprintf('\t\tFORMULA DE LOS TRES OCTAVOS DE SIMPSON\n')
funcion=input('ingrese la funcion \n f(x)=','s');
b=input('ingrese el limite superior de la funcion\n');
a=input('ingrese el limite inferior de la integral\n');
h=(b-a)/3;
x=a; 1
f=eval(funcion);x=a+h; f=f+3*(eval(funcion));
x=a+2*h; f=f+3*(eval(funcion)); x=b;
f=f+eval(funcion);
UNIVERSIDAD CENTRAL DE VENEZUELA.
FACULTAD DE INGENIERÍA-DEPARTAMENTO DE MATEMÁTICAS APLICADAS –CÁLCULO NUMÉRICO.
TEMA 5: DIFERENCIACIÓN E INTEGRACIÓN NUMÉRICA

PASO 2: Seleccione en la interfaz de Matlab la opción RUN. Aparecerá en el área de trabajo un


mensaje como se muestra a continuación:

Proceda a seleccionar el método que desee aplicar, tomando el número que se encuentra a la
izquierda de la descripción del método y presionando enter.

PASO 3: Deberá introducir la siguiente información:

 ingrese la funcion
 ingrese el limite superior de la integral
 ingrese el limite inferior de la integral

Matlab arrojará el valor numérico de la integral según el método seleccionado inicialmente.

Ejemplo 1: Utilice Matlab para estimar el valor numérico de la integral usando la regla de
Simpson 3/8
2
1

1 x
dx

Solución: Tomando en cuenta la secuencia de pasos descrita, tenemos que:

2
UNIVERSIDAD CENTRAL DE VENEZUELA.
FACULTAD DE INGENIERÍA-DEPARTAMENTO DE MATEMÁTICAS APLICADAS –CÁLCULO NUMÉRICO.
TEMA 5: DIFERENCIACIÓN E INTEGRACIÓN NUMÉRICA

PASO 1:

PASO 2:

3
UNIVERSIDAD CENTRAL DE VENEZUELA.
FACULTAD DE INGENIERÍA-DEPARTAMENTO DE MATEMÁTICAS APLICADAS –CÁLCULO NUMÉRICO.
TEMA 5: DIFERENCIACIÓN E INTEGRACIÓN NUMÉRICA

PASO 3:

También podría gustarte