INFORME N°05

INTRODUCCION

Hidráulica de tuberías representa la diferencia más importante, ya que cubre

aspectos que usualmente son toma de cursos o aun carreras diferentes. Dicho

capítulo aborda las redes de riego a presión, incluyendo las redes de sistemas

de riego localizado de alta frecuencia. El diseño de este tipo de redes no había

'sido tratado en forma coherente en un texto de hidráulíca. a pesar de que tiene

muchos elementos comunes con los sistemas de tuberías abiertas y con las

redes de distribución de agua potable.En el caso del flujo en tuberías actúan tres

fuerzas: fuerzas de presión, fuerzas gravitacionales y fuerzas de fricción. Las dos

primeras tratan de acelerar el flujo y las últimas tratan de frenarlo.

En el caso del flujo uniforme existe un equilibrio entre las fuerzas de fricción, por

un lado, y las fuerzas gravitacionales y de presión, por el otro. Dada la

importancia de las fuerzas de fricción en el problema del flujo uniforme, el cual

es básico para el diseño de sistemas de tuberías, el resto de este capítulo se

dedica a estudiarlas
II. OBJETIVOS

2.1. Objetivo general.

 Realizar la Práctica N°5 de Laboratorio con un programa de simulación llamado Matlab.

2.2. Objetivos Específicos.

 Importar los datos de Excel al Matlab.

 Reproducir en Hoja de Cálculo Excel
 Calculo el gasto sistema mostrado en la figura.
 Calculo el gasto sistema mostrado en la figura
 Cálculo y análisis de redes:método de newton-raphson para solucionar redes de
tuberías
3. MATLAB.
3.1 Ejemplo 4.3. Calculo el diámetro que debe tener una tubería nueva, de cemento
enlucido (K=0.0004 m) para conducir 2 m3/s. la viscosidad del agua es de 1.2 x
10-6 m2/s. la longitud de la tubería dela tubería es de 1000 m. la perdida de carga
admisible. Es de 25 m.
CODIGO

clc;clear;format short, format compact

%ASIGNAR EL VALOR f hasta igualarse.
f=0.0171; %coeficiente de Darcy - weisbach
%DATOS DE ENTRADA
H=25;%perdida de carga
g=9.81; %aceleracion de la gravedad
Ks=0.0004; %Rugosidad absoluta PVC
Q=2; %Caudal
L=1000;%Longitud
v=1.2e-6;%Viscosidad cinematica
%CALCULOS PRELIMIARES
S=H/L;%pendiente
d=(0.0827*f*Q^2/S)^(1/5);
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=Q/A; %velocidad flujo
Re=V*d/v; %Número de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
[f f1 d]
disp(' f f d')
MATLAB
RESULTADO

1.1. Ejemplo 4.4. Que presión se requiere impulsar 20 lps a lo largo de una tubería lisa, horizontal,
de 2” de diámetro. La longitud del tramo es 300 m. la viscosidad de agua es de 10 -6 m2/s.
CODIGO

clc;clear;format short, format compact

%ASIGNAR EL VALOR f hasta igualarse.
f=0.0136; %coeficiente de Darcy - weisbach
%DATOS DE ENTRADA
d=2*2.54/100; %Diametro tuberia
Ks=1.5e-6; %Rugosidad absoluta PVC
Q=20e-3; %Caudal
g=9.81;%Aceleracion de la gravedad
v=1e-6;%Viscosidad cinematica
L=300;%longitud
%CALCULOS PRELIMIARES
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=Q/A; %velocidad flujo
Re=V*d/v; %Número de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-
(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))
-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
hf=f1*L*V^2/(2*g*d);
[f f1 hf]
disp(' f f hf')

MATLAB
RESULTADO
1.2. Ejemplo 4.5. Calculo el gasto sistema mostrado en la figura. La viscosidad del agua es 1.2x10 -
6 m2/s. la tubería es lisa. Considerar únicamente las pérdidas de carga continuas. El diámetro
de ls tubería de descarga es de 2 cm.

CODIGO

clc;clear;format short, format compact

%ASIGNAR EL VALOR f hasta igualarse.
f=0.0186; %coeficiente de Darcy - weisbach
%DATOS DE ENTRADA
hf=5;%Perdida de carga continua
d=2/100; %Diametro tuberia
Ks=0; %Rugosidad absoluta PVC
g=9.81;%Aceleracion de la gravedad
v=1.2e-6;%Viscosidad cinematica
L=4;%longitud
%CALCULOS PRELIMIARES
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=((2*g*d*hf)/(f*L))^0.5; %velocidad flujo
Re=V*d/v; %Número de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
Q=V*A;
[f f1 Q]
disp(' f f Q')
MATLAB

RESULTADO
II. EXCEL.
2.1. Ejemplo 4.3. Calculo el diámetro que debe tener una tubería nueva, de cemento enlucido
(K=0.0004 m) para conducir 2 m3/s. la viscosidad del agua es de 1.2 x 10^-6 m2/s. la longitud
de la tubería de la tubería es de 1000 m. la perdida de carga admisible. Es de 25 m.

Solución

Perdida de carga. H= 25 m
Aceleración de la gravedad. g = 9.81 m2/s
Caudal. Q= 2 m3/s
Longitud. L= 1000 m
Viscosidad cinemática. v= 1.20E-06 m2/s
Rugosidad del material. Ks = 0.0004 m
Coeficiente de Darcy -
Weisbach. f= 0.0171

Pendiente S= 0.025
Diámetro

D= 0.74 m

Área. A= 0.433 m2
Velocidad del flujo V= 4.614 m/s
Número de Reynolds.
Re = 2.86E+06

i xi f(xi) F'(xi) xi+1 er

1 7.647191129 -0.0100270 -1.01E+00 7.63721432 0.130634
2 7.637214317 0.0000000 -1.01E+00 7.63721432 0.000000

X= 7.637
f = 1/x^2 f= 0.0171

D= 0.74 m
CODIGO

%ASIGNAR EL VALOR f hasta igualarse.

f=0.0171; %coeficiente de Darcy - weisbach
%DATOS DE ENTRADA
H=25;%perdida de carga
g=9.81; %aceleracion de la gravedad
Ks=0.0004; %Rugosidad absoluta PVC
Q=2; %Caudal
L=1000;%Longitud
v=1.2e-6;%Viscosidad cinematica
%CALCULOS PRELIMIARES
S=H/L;%pendiente
d=(0.0827*f*Q^2/S)^(1/5);
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=Q/A; %velocidad flujo
Re=V*d/v; %Número de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
[f f1 d]
disp(' f f d')

MATLAB
 RESULTADO

2.2. Ejemplo 4.4. Que presión se requiere impulsar 20 lps a lo largo de una tubería lisa, horizontal,
de 2” de diámetro. La longitud del tramo es 300 m. la viscosidad de agua es de 10 -6 m2/s.

Solución

Diámetro tubería d= 0.0508 m

Aceleración de la gravedad. g = 9.81 m2/s
Caudal. Q= 0.02 m3/s
Longitud. L= 300 m
Viscosidad cinemática. v= 1.00E-06 m2/s
Rugosidad del material. Ks = 1.50E-06 m
Coeficiente de Darcy -
Weisbach. f= 0.0136

Área. A= 0.002 m2
Velocidad del flujo V= 9.868 m/s
Número de Reynolds. Re = 5.01E+05

i xi f(xi) F'(xi) xi+1 er

1 8.574929257 0.0117114 -1.09E+00 8.59E+00 0.125667
2 8.59E+00 0.0000005 -1.09E+00 8.58571914 0.000005
3 8.59E+00 0.0000000 -1.09E+00 8.58571914 -

X= 8.586
f = 1/x^2 f= 0.0136
hf = f * L * V^2 / (2 * g * d) hf = 397.59 m
CODIGO

%ASIGNAR EL VALOR f hasta igualarse.

f=0.0136; %coeficiente de Darcy - weisbach
%DATOS DE ENTRADA
d=2*2.54/100; %Diametro tuberia
Ks=1.5e-6; %Rugosidad absoluta PVC
Q=20e-3; %Caudal
g=9.81;%Aceleracion de la gravedad
v=1e-6;%Viscosidad cinematica
L=300;%longitud
%CALCULOS PRELIMIARES
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=Q/A; %velocidad flujo
Re=V*d/v; %Número de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
hf=f1*L*V^2/(2*g*d);
[f f1 hf]
disp(' f f hf')

MATLAB
 RESULTADO

2.3. Ejemplo 4.5. Calculo el gasto sistema mostrado en la figura. La viscosidad del agua es
1.2x10-6 m2/s. la tubería es lisa. Considerar únicamente las pérdidas de carga
continuas. El diámetro de la tubería de descarga es de 2 cm.

Solución

Diámetro tubería d= 0.02 m

Aceleración de la
gravedad. g = 9.81 m2/s
Perdida de carga por
dif. Cota. hf = 5 m
Longitud. L= 4 m
Viscosidad cinemática. v= 1.20E-06 m2/s
Rugosidad del material. Ks = 0 m
Coeficiente de Darcy -
Weisbach. f= 0.0186

Area. A= 0.000314 m2
Velocidad del flujo V= 5.135 m/s
Número de Reynolds. Re = 8.56E+04
 i xi f(xi) F'(xi) xi+1 er
7.33235575 0.002633
1 1 3 -1.12E+00 7.33E+00 0.032100
0.000000 7.334710
2 7.33E+00 0 -1.12E+00 2 0.000001
0.000000 7.334710
3 7.33E+00 0 -1.12E+00 2 -

X= 7.335
f = 1/x^2 f= 0.0186
Q = V*A Q= 0.00161 m

CODIGO

hf=5;%Perdida de carga continua

d=2/100; %Diametro tuberia
Ks=0; %Rugosidad absoluta PVC
g=9.81;%Aceleracion de la gravedad
v=1.2e-6;%Viscosidad cinematica
L=4;%longitud
%CALCULOS PRELIMIARES
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=((2*g*d*hf)/(f*L))^0.5; %velocidad flujo
Re=V*d/v; %Número de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
Q=V*A;
[f f1 Q]
disp(' f f Q')
MATLAB

RESULTADO
 PRÁCTICA V-(b)
“HIDRAULICA DE TUBERÍAS – CÁLCULO Y ANÁLISIS DE REDES:
MÉTODO DE NEWTON-RAPHSON PARA SOLUCIONAR REDES DE
TUBERÍAS”

1. Ejemplo 3: La red mostrada en la figura tiene una válvula en la tubería 2-3, la cual se encuentra
parcialmente cerrada y produce una perdida menor local de 10𝑉23 2
/2𝑔, la presión en el nodo
1 es de 100 mca. Analizar los caudales y presiones en la red. Los diámetros (en milímetros) y
las longitudes (en metros) para cada una de las tuberías son los indicados en la figura. Los
caudales de demanda en los nodos están dados en l/s., la viscosidad cinemática es 𝑣 =
1.141𝑥10−6 𝑚2/𝑠 y Rugosidad absoluta en todas la tuberías: Ks=0.00005.
METODO DE NEWTON-RAPHSON PARA SOLUCIONAR REDES
DE TUBERIAS
PRIMERA ITERACIÓN
Rugosidad Ks (m): 6.00E-05
Viscocidad V (m2/s): 1.14E-06
Gravedad g (m/s2): 9.81

Nodo 1 2 3 4 5 6
Presión (m) 100 92 80 75 90 95

Primera iteración Ec. 7.36 Ec. 7.37 Ec. 7.38 Ec. 7.33
∑

Tubería Hi (m) Hj (m) hf (m) L (m) d (m) A (m2) Q (m3/s) hm (m) km+f*l/d df/dx
12 100 92 8.0 500 0.25 0.0491 0.1093 31.6333 0.006834
23 92 80 12.0 400 0.15 0.0177 0.0395 2.5496 47.0666 0.001647
9.450 400 0.15 0.0177 0.0349 1.9865 60.4068 0.001454
10.013 400 0.15 0.0177 0.0360 2.1105 56.8583 0.001498
9.889 400 0.15 0.0177 0.0357 2.0832 57.6039 0.001489
9.917 400 0.15 0.0177 0.0358 2.0892 57.4380 0.001491
9.911 400 0.15 0.0177 0.0358 2.0879 57.4745 0.001490
9.912 400 0.15 0.0177 0.0358 2.0882 57.4664 0.001490
34 80 75 5.0 200 0.10 0.0079 0.0123 40.0030 0.001230
45 75 90 15.0 400 0.15 0.0177 0.0444 46.6334 0.001480
25 92 90 2.0 200 0.10 0.0079 0.0076 42.3564 0.001890
56 90 95 5.0 600 0.20 0.0314 0.0431 52.0718 0.004312
16 100 95 5.0 300 0.25 0.0491 0.1117 18.9475 0.011169

Ec. 7.35 Ec. 7.29

∑

Nodo N° Tubería N° Hi (m) Hj (m) hf (m) Q (m3/s) df/dx df/dH f= Q-Qd

2 -0.0600
21 92 100 8 0.1093 0.006834
23 92 80 -12 -0.0358 0.001647
-
25 92 90 -2 -0.0076 0.001890 0.01037065 0.0060148
3 -0.0400
32 80 92 12 0.0358 0.001647
-
34 80 75 -5 -0.0123 0.001230 0.00287674 -0.0165303
4 -0.0300
43 75 80 5 0.0123 0.001230
-
45 75 90 15 0.0444 0.001480 0.00270970 0.0266927
5 -0.0300
52 90 92 2 0.0076 0.001890
54 90 75 -15 -0.0444 0.001480
-
56 90 95 5 0.0431 0.004312 0.00768171 -0.0237135
6 -0.0400
61 95 100 5 0.1117 0.011169
-
65 95 90 -5 -0.0431 0.004312 0.01548140 0.0285732

En forma matricial [Df/Dx][Dx]=[-F]

Nodos 2 3 4 5 6
2 -0.010371 0.001647 0.00 0.001890 0.00 dH2 -0.0060148
3 0.001647 -0.002877 0.001230 0.00 0.00 dH3 0.0165303
 4 0.00 0.001230 -0.002710 0.001480 0.00 dH4 = -0.0266927
5 0.001890 0.00 0.001480 -0.007682 0.004312 dH5 0.0237135
6 0.00 0.00 0.00 0.004312 -0.015481 dH6 -0.0285732

Solución de Dx=[Df/Dx]-1*[-F] por Inversa

-
-122.276563 -102.159060 -75.222555 -52.834235 14.715942 0.222237475
-
-102.159060 -531.988798 -294.453905 -97.026387 27.024802 -1.848374766
-
[Df/Dx]-1= -75.222555 -294.453905 -587.995606 -156.197849 43.505855 = 8.819306355
-
-52.834235 -97.026387 -156.197849 -205.378307 57.204109 -0.352486556
-
-14.715942 -27.024802 -43.505855 -57.204109 80.526720 1.747467683

Los resultados para el vector dHi son: Luego las nuevas alturas piezométricas son:

dH2= 0.222237 H2= 92.2222

dH3= -1.848375 H3= 78.1516
dH4= 8.819306 H4= 83.8193
dH5= -0.352487 H5= 89.6475
dH6= 1.747468 H6= 96.7475
METODO DE NEWTON-RAPHSON PARA SOLUCIONAR REDES
DE TUBERIAS
SEGUNDA ITERACIÓN
Rugosidad Ks (m): 6.00E-05
Viscocidad V (m2/s): 1.14E-06
Gravedad g (m/s2): 9.81

Nodo 1 2 3 4 5 6
Presión
(m) 100 92.222 78.152 83.819 89.648 96.747

Primera iteración Ec. 7.36 Ec. 7.37 Ec. 7.38 Ec. 7.33
∑

Tubería Hi (m) Hj (m) hf (m) L (m) d (m) A (m2) Q (m3/s) hm (m) km+f*l/d df/dx
12 100 92.222 7.8 500 0.25 0.0491 0.1077 31.6712 0.006927
23 92.222 78.152 14.1 400 0.15 0.0177 0.0429 3.0095 46.7538 0.001526
11.061 400 0.15 0.0177 0.0379 2.3418 60.0856 0.001346
11.729 400 0.15 0.0177 0.0391 2.4895 56.5196 0.001388
11.581 400 0.15 0.0177 0.0388 2.4568 57.2722 0.001379
11.614 400 0.15 0.0177 0.0389 2.4640 57.1038 0.001381
11.607 400 0.15 0.0177 0.0388 2.4624 57.1410 0.001380
11.608 400 0.15 0.0177 0.0388 2.4628 57.1328 0.001380
34 78.152 83.819 5.7 200 0.10 0.0079 0.0131 39.7359 0.001159
45 83.819 89.648 5.8 400 0.15 0.0177 0.0271 48.7421 0.002322
25 92.222 89.648 2.6 200 0.10 0.0079 0.0087 41.6328 0.001680
56 89.648 96.747 7.1 600 0.20 0.0314 0.0519 51.0990 0.003653
16 100 96.747 3.3 300 0.25 0.0491 0.0892 19.3152 0.013716

Ec. 7.35 Ec. 7.29

∑

Nodo N° Tubería N° Hi (m) Hj (m) hf (m) Q (m3/s) df/dx df/dH f= Q-Qd

2 -0.0600
21 92.22224 100 7.77776 0.1077 0.006927
23 92.22224 78.15163 -14.07061 -0.0388 0.001526
-
25 92.22224 89.64751 -2.57472 -0.0087 0.001680 0.01013273 0.0002529
3 -0.0400
32 78.15163 92.22224 14.07061 0.0388 0.001526
-
34 78.15163 83.81931 5.66768 0.0131 0.001159 0.00268499 0.0119837
4 -0.0300
43 83.81931 78.15163 -5.66768 -0.0131 0.001159
-
45 83.81931 89.64751 -5.82821 0.0271 0.002322 0.00348114 -0.0160718
5 -0.0300
52 89.64751 92.22224 2.57472 0.0087 0.001680
54 89.64751 83.81931 -5.82821 -0.0271 0.002322
-
56 89.64751 96.74747 7.09995 0.0519 0.003653 0.00765500 0.0034552
6 -0.0400
61 96.74747 100 3.25253 0.0892 0.013716
-
65 96.74747 89.64751 -7.09995 -0.0519 0.003653 0.01736892 -0.0026469

En forma matricial [Df/Dx][Dx]=[-F]

Nodos 2 3 4 5 6
 2 -0.010133 0.001526 0.00 0.001680 0.00 dH2 -0.0002529
3 0.001526 -0.002685 0.001159 0.00 0.00 dH3 -0.0119837
4 0.00 0.001159 -0.003481 0.002322 0.00 dH4 = 0.0160718
5 0.001680 0.00 0.002322 -0.007655 0.003653 dH5 -0.0034552
6 0.00 0.00 0.00 0.003653 -0.017369 dH6 0.0026469

Solución de Dx=[Df/Dx]-1*[-F] por Inversa

- - - -
-122.358802 98.98319839 68.20987449 52.84894262 11.11473234 0.274054911
- - - -
98.98319839 537.3059217 251.5959087 108.9797891 -22.9196863 2.736219058
- - - - -
[Df/Dx]-1= 68.20987449 251.5959087 493.0184588 182.8744679 38.46057577 = -4.361348604
- - - - -
52.84894262 108.9797891 182.8744679 219.7600227 46.21802651 -0.982806526
- - - -
11.11473234 -22.9196863 38.46057577 46.21802651 67.29429214 -0.359086195

Los resultados para el vector dHi son: Luego las nuevas alturas piezométricas son:

dH2= 0.274055 H2= 92.4963

dH3= 2.736219 H3= 80.8878
dH4= -4.361349 H4= 79.4580
dH5= -0.982807 H5= 88.6647
dH6= -0.359086 H6= 96.3884

CODIGO

clc;clear;format short;format compact

%DATOS DEL SISTEMA

H=43.5;
MATLAB %Carga Hidraulica
L=730; %Longitud de tuberia
d=0.293;
RESULTADO %Diámetro
Ks=1.5e-6; %Rugosidad
g=9.81; %Aceleracion de gravedad
v=1.007e-6; %Viscosidad Cinematica
Km=11.80; %Coeficiente de perdida menor

%CALCULOS PRELIMINARES
A=pi*d^2/4; %Area de tuberia
i=1;hfs(i)=0;
%CONDICIONES ASUMIDAS
tol=0.1;er=100;
while er(i)>=tol
hf(i)=H-hfs(i);
Q(i)=(pi*d^2/4)*(-2*sqrt(2*g*d*hf(i)/L))*...
log10(Ks/(3.71*d)+2.51*v/(d*sqrt(2*g*d*hf(i)/L)));
V(i)=Q(i)/A; Re(i)=V(i)*d/v;
hfs1(i)=hfs(i);
hfs(i+1)=Km*V(i)^2/(2*g);
er(i+1)=abs((hfs(i+1)-hfs1(i))/hfs(i+1))*100;
i=i+1;
end

TABLA =[hf' Q' V' hfs1']

MATLAB

RESULTADO
IV. CONCLUCIONES YRECOMENDACIONES

4.1 CONCLUCIONES
 Se estableció la Práctica N°5 de Laboratorio con un programa de simulación
llamado Matlab.
 Se importaron los datos de Excel al Matlab.
 Se ingresó los datos del archivo Excel con el nombre de Newton Rapson.
4.2 RECOMENDACIONES

 Crear una carpeta para guardar los archivos del Matlab, de esa manera se
trabajará ordenadamente.