Metodos Numericos

Metodos Numericos (MATLAB)- VIDEOTUTORIAL
Posted by Rosemberg Benavides on 18:03
[pic]

PARA VER EL VIDEOTUTORIAL DAR CLIC EN LA IMAGEN 
|[pic]|
|[pic][pic][pic] |
|[pic]|
|DAR CLIC EN LA IMAGEN PARA REPRODUCIR EN HD |
|Si no quieresesperar presiona (Cerrar X). |

METODO DE EULER
clear all
disp('METODO DE EULER')
clc
syms x
syms y
f=inline(input('ingrese laderivada:','s'));
x=input('ingrese el valor de x:');
y=input('ingrese el valor de y:');
h=input('ingrese el valor de h:');
n=input('ingrese numero de iteraciones:');
clc
disp('x(n) y(n) y´(n)hy´(n)');
for i=1:n
y1=feval(f,x,y);
hy1=h*y1;
fprintf('\n%0.1f %0.4f %0.4f %0.4f ',x,y,y1,hy1);
y=y+hy1;
x=x+h;
end
METODO DE NEWTON-COTES
clc
clear all
disp('METODO DE NEWTON-COTES')
syms xf=inline(input('ingrese la funcion:','s'));
a=input('ingrese intervalo inferior:');
b=input('ingrese intervalo superior:');
n=1
h=(b-a)/3
x0=a
x1=x0+h
x2=x1+h
x3=x2+h
x4=x3+h
f0=feval(f,x0)f1=feval(f,x1)
f2=feval(f,x2)
f3=feval(f,x3)
f4=feval(f,x4)
I1=(h/2)*(f0+f1)
error=(-1/12)*(h^5)*((f2-(2*f1)+f0)/(h^2))
I1=I1+error
n=2
I2=(h/3)*(f0+4*f1+f2)error=(-1/90)*(h^5)*((f4-(4*f3)+(6*f2)-(4*f1)+f0)/(h^4))
I2=I2+error
METODO DE EULER MODIFICADO
clear all
disp('METODO DE EULER MODIFICADO')
clc
syms x
syms y
f=inline(input('ingrese la derivada:','s'));
x=input('ingrese elvalor de x:');
y=input('ingrese el valor de y:');
h=input('ingrese el valor de h:');
n=input('ingrese numero de iteraciones:');
clc
disp('x(n) y´(n) hy´(n) y(n+1),p hy´(n+1),p y(n+1),c');
for...