-
Notifications
You must be signed in to change notification settings - Fork 0
/
MN_newtonRapshon.m
55 lines (38 loc) · 1.13 KB
/
MN_newtonRapshon.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
clc;clear;close all
%Metodo numerico de NewtonRaphson derivada analitica
%Autor= Juan Felipe Hernandez Arango. Ing.Quimico
%Parametros iniciales del metodo
Error=1;
tol=1e-6;
contador=0;
% valores incial para el metodo abierto
xr=12;
[solucion,Error,contador]=NewtonRapshon(xr,tol,Error,contador)
% Se debe colocar a continuacion la funcion a solucionar igualada a 0
% junto con los parametros necesarios para solucionar la funcion
function [y]= Ecuacion (x)
g=9.8;
m=68.1;
t=10;
v=40;
y=((g*m)/x)*(1-exp(-(x/m)*t))-v;
end
% En esta funcion debe ir la derivada analitica junto con los parametros
% necesarios
function [y]= derivada (x)
g=9.8;
m=68.1;
t=10;
v=40;
y=(-g*m*(m*exp(x*t/m)-x*t-m)*exp(-x*t/m))/(m*x^2);
end
% Metodo de Newton Rapshon con derivada analitica
function [xr,Error,contador] = NewtonRapshon(xr,tol,Error,contador)
while Error>tol
% se aplica la formula de Newton Rapshon para aproximarce a la raiz
xsol=xr-(Ecuacion(xr)/derivada(xr));
Error=abs((xsol-xr)/xsol); % se calcula el error relativo
contador=1+contador;
xr=xsol;
end
end