1

FACULTAD DE INGENIERIA



ESCUELA ACADEMICO PROFESIONAL DE
INGENIERA CIVIL


ASIGNATURA: INGENIERA MATEMTICA II



GUA DE LABORATORIO


AUTOR: Dr. LUIS PAIHUA MONTES





SEMESTRE ACADMICO 2013-0


2
















Este material de enseanza se hace en concordancia con lo dispuesto por la legislacin
sobre derechos de autor:

Ley 13714
Art 69.- Pueden ser reproducidos y difundidos breves fragmentos de obras literarias,
cientficas y artsticas, y an la obra entera, si su breve extensin y naturaleza lo
justifican siempre que la reproduccin se haga con fines culturales y no comerciales, y
que ella no entrae competencia desleal para el autor en cuanto a aprovechamiento
pecuniario de la obra, debiendo indicarse, en todo caso, el nombre del autor, el ttulo de
la obra y la fuente de donde se hubieren tomado.






















3

INGENIERA MATEMTICA II
Laboratorio N 1
Introduccin

OBJETIVO GENERAL:

Familiarizarse con la forma de trabajo en el Mathcad.
Verificar la facilidad de trabajo y potencia del asistente.
Trabajar en el plano complejo con el Mathcad.
Presentacin en 2D.
Presentacin en 3D

COMPETENCIAS:

Edita textos, ecuaciones.
Comprende y desarrolla el clculo aritmtico y formal.
Define objetos en el Mathcad y opera con ellos.
Edita nmeros complejos y los opera.
Grafica en 2D.
Grafica curvas parametrizadas con dominio predeterminado.
Grafica en 3D, en sus diferentes presentaciones.
Parametriza superficies y grafica en 3D.

INTRODUCCIN:

Con frecuencia se har mencin de las herramientas de trabajo que estn agrupadas en:
CALCULADORA, GRFICO, MATRIZ, EVALUACIN, CLCULO, BOOLEANO,
PROGRAMACIN, GRIEGO, SIMBLICO.

Inspeccionar estas barras de herramienta.

Para escribir un texto: Posicione el cursor en cualquier lugar de la pantalla, luego haga
MAYSCULA + 2 ([ ] doble comilla), enseguida escriba la frase:

Ingeniera Matemtica II con Mathcad
Ciclo 2013-0

Luego subraye y resalte para que obtener:

Ingeniera Matemtica II con Mathcad
Ciclo 2013-0

Despus desplace arrastrando el ratn de manera que est centrado en la parte superior.

Posicionar el cursor en forma adecuada y escriba lo indicado sin dejar espacio libre (de corrido),
observe el resultado.


1 1
7 4
+ = (resultado aritmtico, en decimal)

1 1
7 4
+ (resultado formal, exacto)

3
1 2
2
3
4 1
cosh 1
2 1
1 ln 3 1
1
x
e x
x x
x x

+
| |

|

| |
\ .
+
|
+
\ .
(una expresin algebraica)

4

3 2 3
3 z x y y + = 2 (una ecuacin, el signo igual es el lgico)


Para definir los objetos en el Mathcad: variables, funciones, matrices, constantes, etc se emplea la
siguiente sintaxis:
Objeto:=valor

Como aplicacin defina las matrices y las funciones dadas:

Para crear matrices emplear las herramientas en el grupo MATRIZ , no olvidar que la
enumeracin por defecto empieza en cero para modificar haga ORIGIN:=1
1 2 3 4
5 6 7 8
:
9 0 1 2
3 4 5 6
A
(
(

(
=
(
(

1 3 5 7
9 0 2 4
:
6 8 1 2
3 4 5 6
B
(
(
(
=
(
(

. Haga las operaciones bsicas, y visualizar
el elemento de la fila 2, columna 3 haciendo A
2,3
.

2
( ): 3 5 f x x x = (funcin real de variables real)
Para evaluar la funcin definida en un valor determinado por ejemplo x=1 es suficiente hacer
f(1)= (ver lo que se obtiene)

2 3
( , ): 5 2 f x y x y = (funcin real de dos variables
Para evaluar de manera similar al caso anterior escribir el valor correspondiente a cada variable,
por ejemplo hacer f(3,2)= (verificar el valor hallado)

Escriba lo siguiente: z1:=4-3i z2:=-5+4i, (emplear el i que se encuentra en la calculadora) luego
realice la operaciones.

a. Suma z1+z2
b. Multiplicacin z1.z2
c. Divisin
1
2
z
z

d. Conjugado 1 z (emplear MAY+2)
e. |z2|

Para graficar una funcin real de variable real emplear el cono que se encuentra en
GRAFICO, como aplicacin graficar ambas funciones
2
2
( ) 3 , ( ) 1
9
x
f x g x x = = en el
intervalo [-3 , 4] con una presentacin adecuada.

Para graficar una curva determine su ecuacin paramtrica, es decir escribir
0 0
( )
( ) , [ , ]
x f t
y g t t a b
=

= e

y hacer:









Aplicacin: Parte de la elipse
2 2
9 4 36 x y + = que se encuentra en el 2do y 3er cuadrantes.
a:=a
0
b:=b
0
n:=200 :
b a
h
n

= k:=0 .. n :
k
t a k h = +
: ( )
: ( )
k k
k k
x f t
y g t
=
=

A continuacin graficar, este par escribiendo x
k
en el dominio, y
k
en el rango
5

Para graficar una superficie que es grfica de una funcin real de dos variables emplear , que
se encuentra en GRAFICO, aplicar a las funciones.
















Y proceder de acuerdo a la indicacin del profesor, haga varias presentaciones

Si la superficie est parametrizada
( , ):
( , ):
( , ):
x r t
y r t
z r t
=

con el mismo comando anterior en la parte del

rectngulo rojo escribirla terna (x,y,z), enseguida escoger una buena presentacin con el cuadro de
dilogo de formato 3D. Aplicacin

2 2 2 2
1 2
: , : 8 S z x y S z x y = + = .





























f x y , ( ) x
2
y
2
+ :=

g x y , ( ) 8 x
2
y
2
:=



x1r t , ( ) r cos t ( ) :=

y1 r t , ( ) r sin t ( ) :=

z1r t , ( ) r
2
:=


x2r t , ( ) r cos t ( ) :=

y2 r t , ( ) r sin t ( ) :=

z2r t , ( ) 8 r
2
:=




x1y1 , z1 , ( ) x2y2 , z2 , ( ) ,


6


INGENIERA MATEMTICA II
Laboratorio N 2
Mapeo y Flujo

OBJETIVO GENERAL:

Mapeo en forma grfica.
Aplicar la funcin compleja en explicar el comportamiento del flujo de fluido.
Explicar el potencial complejo de flujo

COMPETENCIAS:

Determina la imagen de curvas.
Determina la imagen de regiones limitadas por curvas.
Descompone una funcin compleja en parte real e imaginaria.
Halla la funcin corriente para una funcin compleja
Graficar las lneas de corriente e interpretar el flujo.
Identifica la funcin corriente y potencial en el potencial complejo de flujo.
Grafica las Lneas de Corriente y curvas equipotenciales e interpretar el flujo.
Calcula el vector velocidad.


1. Para descomponer en parte real e imaginaria una funcin compleja en el Mathcad proceder como
se indica en el cuadro siguiente, es importante las tareas que se encuentran en SIMBOLICO para
desarrollar esta parte. La unidad imaginaria es el que se encuentra el CALCULADORA.
















2. Para hallar la imagen de una curva proceda como sigue::

Hallar la ecuacin paramtrica de la curva
1( ): ( )
1( ): ( )
x t f t
y t g t
=

, y a continuacin de lo
desarrollado anteriormente agregar:








Graficar poniendo ( 1 , 1 ) ( 1 , 1 )
k k k k
u x y v x y en los ejes correspondientes.
f z ( ) z
2
:=
Desarrollo de la funcin
g x y , ( ) f z ( )
reemplazar z x y i + ,
expandir
simplificar
x
2
y
2
2i xy + :=
Parte real y parte imaginaria
u x y , ( )
g x y , ( ) g x y , ( )

+
2
complejos x
2
y
2
:=
v x y , ( )
g x y , ( ) g x y , ( )

2i
complejos 2x y :=

a:=a
0
b:=b
0
n:=200 :
b a
h
n

= k:=0 .. n :
k
t a k h = +
1 : ( ) 1 : ( )
k k k k
x f t y g t = =

7
Si hay ms curvas proceda como en el caso anterior pero denotar por
2( ):
2( ):
x t
y t
=

,
3( ):
3( ):
x t
y t
=

las otras curvas (respectivamente), para dibujar varias imgenes por ejemplo si fueran tres
curvas en el cuadro del dibujo escriba:

En el dominio ( 1 , 1 ), ( 2 , 2 ), ( 3 , 3 )
k k k k k k
u x y u x y u x y
En el rango ( 1 , 1 ), ( 2 , 2 ), ( 3 , 3 )
k k k k k k
v x y v x y v x y

Como aplicacin hallar la imagen de la figura cuando aplicamos la funcin compleja f(z)=z
2
-z





















3. Dada la funcin compleja
1
( ) f z
z
= describa el flujo asociado a esta funcin.
4. Dado el potencial complejo
2
( ) F z z z = + explicar la naturaleza del flujo.




























2 2
4 16 x y + =
2 2
( 1) ( 1) 16 x y + =






2 2
16 x y + =


8


INGENIERA MATEMTICA II
Laboratorio N 5
Proceso estable inestable y Propagacin de error

OBJETIVO GENERAL:

Identificar si el proceso es estable o inestable
Estimar la propagacin de error en cierta situacin.

COMPETENCIAS:

Distingue cuando hay estabilidad o inestabilidad.
Modela ciertas situaciones y analiza la propagacin de error.

PROCESO ESTABLE-INESTABLE
1. Haga el grfico simultneo de las siguientes funciones polinmicas.
5
6
0
( ): ( ) , ( ): ( ) 0.001
j
P x x j Q x P x x
=
= =
[

Luego tramo a tramo analizar la diferencia que existe entre ellos y comentar a que se debe esa
diferencia.

PROPAGACIN DE ERROR
2. La frmula de Manning para un canal rectangular se escribe como
2 / 1
3 / 2
3 / 5
) 2 (
) ( 1
S
H B
BH
n
Q
+
= ,
donde Q=flujo (m
3
/s), n=coeficiente de rugosidad, B=ancho (m), H=profundidad (m) y
S=pendiente. Estimar Q si tenemos un arroyo con B=20 m H=0.3 m , n=0.03 10% ,
S=0.0003 10%. Cul de estos valores con error recomienda sea medido con ms
cuidado?




Evaluacin de Laboratorio N 1

Las preguntas medirn las siguientes competencias

Grafica en 2D y 3D
Mapeo de curvas y regiones del plano.
Descompone una funcin compleja en parte real e imaginaria
















9



INGENIERA MATEMTICA II
Laboratorio N 4
Solucin numrica de una EDO

OBJETIVO GENERAL:

Resolver numricamente un EDO de orden 1 o ms.
Resolver una EDO con condiciones de frontera, empleando la diferencia finita

COMPETENCIAS:

Aplica el mtodo de Euler y resuelve una EDO.
Emplea el mtodo de RK$ y resuelve una EDO.
Plantea una EDO de orden n>1 como un sistema de EDO de orden 1.
Resuelve un sistema de EDO.


1. Prepara la siguiente plantilla que corresponde al mtodo de Euler
































Resuelva el mismo ejercicio con Runge-Kutta de orden 4, archivo del aula virtual, estime el error
segn las indicaciones del profesor.

Para resolver una EDO de orden superior y/o sistema de EDO, seguir las indicaciones del profesor
y emplear los archivos del aula virtual.

Mtodo de Euler
f x y , ( )
x
y
y
x
+ :=
a 1 :=

b 3 :=

n 10 :=

k 1 n .. :=
x
0
a :=

y
0
0.5 :=
h
b a
n
:=

x
k
a k h + :=

y
k
y
k 1
h f x
k 1
y
k 1
,
( )
+ :=

SOL augment x y , ( ) :=


SOL
0 1
0
1
2
3
4
5
6
7
8
9
10
1 0.5
1.2 1
1.4 1.406666667
1.6 1.80667118
1.8 2.209626407
2 2.618063921
2.2 3.032654971
2.4 3.453438269
2.6 3.880216707
2.8 4.312708049
3 4.750607432
=

10
2. EDO con condiciones de frontera, aplicacin de la diferencia finita.


3. La ecuacin diferencial que explica el movimiento de un resorte rgido est dad por la ecuacin
3
'' y ky py = , donde k y p son constantes. Resuelva la ED para k=4 p=1, teniendo las
condiciones iniciales (0) 0 , '(0) 1 y y = = (haga lo mismo para '(0) 0, 2, 3 y = )

4. La siguiente EDO modela el comportamiento de un resorte envejecido
/100
'' 0.5 0
(0) 2 , (2) 1
t
y e y
y y

+ =

= =

, resolver y explicar el significado cuando t es grande.

Evaluacin de Laboratorio N 2

Las preguntas medirn las siguientes competencias

Numricamente explica el proceso estable inestable.
Analiza la propagacin de error.
Explica el comportamiento del flujo de un fluido modelado por una funcin compleja.




































11



INGENIERA MATEMTICA II
Laboratorio N 5
Ecuacin no lineal


OBJETIVO GENERAL:

Plantear procesos iterativos convergentes.
Aplicar correctamente los mtodos que resuelvan una ecuacin no lineal.
Resolver numricamente una ecuacin no lineal.

COMPETENCIAS:

Localiza cada solucin de una ecuacin.
Distingue si el proceso iterativo es convergente o no.
Verifica las condiciones favorables de cada mtodo al aplicarla.
Calcula la solucin de una ecuacin no lineal con cierta exigencia.

1. En cada caso identificar la funcin a partir de la ecuacin ( ) 0 f x = , luego hacer el grfico como
se indica a continuacin para localizar las soluciones deseadas.
Para ello hacer:


















2. En cada caso, aplicar un mtodo segn indicacin del docente y determinar la solucin para la
localizacin hallada.

a. La ecuacin que permite determinar la frecuencia de la viga empotrada en un extremo y libre en
el otro extremo, est dada por cos( )cosh( ) 1 0 ul ul + = , donde:
2
p
u
a
= ,
2
EIg
a
A
=
p=Frecuencia circular natural de la viga, radianes /s.
l=300cm (longitud de la viga)
I=7000cm
4
(momento de inercia del rea del material de la viga)
E=2 (10)
5
kg/cm
2
(mdulo de elasticidad del material de la viga)
0.002 = kg/cm
3
(densidad del material de la viga)
A=200 cm
2
(rea de la seccin transversal de la viga)
g=aceleracin de la gravedad, cm/s
2


f x ( ) :=
h 1 :=
a 0 :=
b a h + :=
f x ( )
b a x
f


Observar que la grfica se realiza en el
intervalo [a , b] de amplitud h=1.

Cambiando el valor de a para distintos
valores se realiza una exploracin que
nos permite localizar cada solucin.

Seguir las recomendaciones que el
profesor dar.

12
Localizar las tres primeras soluciones y calcular empleando Secante para la 1ra, Newton para la
2da, Aproximacin Sucesiva para la 3ra, con un error de la milsima.


b. La ecuacin para calcular la curva elstica resultante en una viga uniformemente sujeta
a una carga distribuida creciente linealmente, como se muestra en la figura, est dada
por la frmula ) 2 (
120
4 3 2 5 0
x L x L x
EIL
w
y + =
Numricamente halle el punto de mxima deflexin, y dicha deflexin. Emplear
L=450 cms, E=50000 kN/cm
2
, I=30000 cm
4
y w
0
=1.75 KN/cm


















Evaluacin de Laboratorio N 3

Las preguntas medirn las siguientes competencias

Resolucin numrica de una EDO de orde 1 y superior con C.I.
Resolucin numrica de una EDO con C.F.
























W
0





L

13


INGENIERA MATEMTICA II
Laboratorio N 6
Aproximacin

OBJETIVO GENERAL:

Aproximar el comportamiento de un fenmeno mediante un polinomio

COMPETENCIAS:

Calcula la tabla de diferencia dividida.
Construye el polinomio de interpolacin.
Planifica y construye el ajuste lineal mediante un polinomio.
Planifica y construye el ajuste lineal no lineal.

Empleando adecuadamente los archivos del aula virtual

1. Dada la siguiente informacin:

x 1 2 2.5 3 4 5
f(x) 1 5 7 8 2 1

a. Calcule la tabla de diferencia dividida con todos los datos.
b. Escriba el polinomio de interpolacin de grado 1, 2, 3 que permita aproximar f(3.4), estime el
error correspondiente.

2. Para cada caso que sigue, emplear las herramientas del aula virtual, segn la necesidad.

a. El esfuerzo cortante, en kips por pie cuadrado (kpc) de nueve muestras tomadas a distintas
profundidades en un estrato de arcilla son:

Prof.(m)
1.9 3.1 4.2 5.1 5.8 6.9 8.1 9.3 10.0
Esf.(kpc)
0.3 0.6 0.4 0.9 0.7 1.1 1.5 1.3 1.6

Estime el esfuerzo cortante a una profundidad de 4.5 m.

b. Se desea determinar la relacin entre el flujo de un ro y la precipitacin en la zona donde se
encuentra el ro, para ello se dispone de la siguiente informacin

Preci. cm)
88.9 101.6 104.1 139.7 132.1 94.0 116.8 121.9 99.1
Flujo m
3
/s
114.7 172.
0
152.9 269.0 206.4 161.4 175.8 121.9 99.1

Use la mejor curva de ajuste y predecir el flujo anual de agua si la precipitacin es de 120 cm


Evaluacin de Laboratorio N 4

Las preguntas medirn las siguientes competencias

Resolucin numrica de una ecuacin no lineal mediante diferentes mtodos






14



INGENIERA MATEMTICA II
Laboratorio N 7
Cuadratura

OBJETIVO GENERAL:

Determinar el rea de regiones del plano

COMPETENCIAS:

Calcula la integral definida en forma aproximada cuando se conoce la funcin empleando
diferentes mtodos con la ayuda del asistente, estimando el error.
Calcula la integral empleando informacin discreta y estima el error.


Para desarrollar los ejercicios haga las siguientes plantillas.































Para los datos discretos adaptar la plantilla con la recomendacin del profesor.





Trapecio Cerrado
f x ( )
1
1 x
3
+
:=
a 0 :=

b 2 :=

n 10 :=

k 0 n .. :=

h
b a
n
:=

x
k
a k h + :=
TC
h
2
1
n
k
f x
k 1
( )
f x
k
( )
+
( )

=
:=

Valor TC 1.089495609733241 =
Trapecio Abierto
f x ( )
1
1 x
3
+
:=
a 0 :=

b 2 :=

n 10 :=

k 0 n .. :=
h
b a
n
:=

x
k
a k h + :=


TA
h
2
1
n
k
f x
k 1
h
3
+
|

\
|
|
.
f x
k
h
3

\
|
|
.
+
|

\
|
|
.

=
:=

Valor
TA 1.090172258335061 =
Simpson Cerrado
f x ( )
1
1 x
3
+
:=
a 0 :=

b 2 :=

n 4 :=

k 0 n .. :=
h
b a
n
:=

x
k
a k h + :=
SC
h
6
1
n
k
f x
k 1
( )
4f x
k 1
h
2
+
|

\
|
|
.
+ f x
k
( )
+
|

\
|
|
.

=
:=


Valor SC 1.090097506764173 =

Simpson Abierto
f x ( )
1
1 x
3
+
:=
a 0 :=

b 2 :=

n 4 :=

k 0 n .. :=
h
b a
n
:=

x
k
a k h + :=

SA
h
3
1
n
k
2f x
k 1
h
4
+
|

\
|
|
.
f x
k 1
h
2
+
|

\
|
|
.
2f x
k
h
4

\
|
|
.
+
|

\
|
|
.

=
:=


Valor SA 1.089920959442364 =


15




2.-





















3.- Determinar el rea de la regin comprendida entre la carretera y el cause del ro
mostrado en la figura, la malla es de un cm por lado y la escala es 1/10 000.


























Para el diseo estructural de un velero de
carreras es necesario determinar la fuerza
total del viento ejercida sobre el mstil,
empleando la fuerza distribuida f es
necesario hallar esta fuerza total mediante la
integral:
}
=
30
0
) ( dz z f F donde
15 /
5
200 ) (
z
e
z
z
z f

|
.
|

\
|
+
=
Adems la lnea de accin F se calcula
evaluando la integral
}
}
=
30
0
30
0
) (
zf(z)dz
dz z f
d .
Determine F y d con un error del 1%





T
30 pies







3 pies







RIO





Carretera
16



4. La figura muestra la seccin transversal de un canal, los valores numricos estn en
metros, determinar el rea de esta seccin transversal, y estime el error.



















Evaluacin de Laboratorio N 5

Las preguntas medirn las siguientes competencias

Construir un polinomio interpolante y estimar el error.
Construir una curva de ajuste adecuado.





























1.8 3.9 5.8 6.5 6 5.7 4.9 3.4 1.7








0 10 m 20 m
17





INGENIERA MATEMTICA II
Laboratorio N 8
Cubicacin



1. Dada la siguiente informacin: Calcular
} }
50
0
50
0
) , ( dxdy y x f (cubicar)

Z=f(x,y)
x \ y
0 5 10 15 20 25 30 35 40 45 50
0 1.5 2.5 4.2 6.5 9.2 14.5 16.0 10.5 7.5 4.6 2.5
5 3.0 5.1 6.5 8.1 9.7 11.5 12.5 10.2 8.5 6.5 3.7
10 4.5 6.9 7.3 9.5 10.5 8.7 9.5 7.5 6.2 5.6 4.2
15 7.2 8.3 8.5 8.9 11.1 9.2 10.5 10.9 8.6 4.5 7.5
20 8.0 9.6 10.5 13.2 12.3 13.7 14.0 11.0 9.5 3.5 6.5
25 9.5 10.5 12.5 14.2 15.0 16.5 15.0 12.1 10.1 5.5 6.2
30 7.6 9.0 10.5 11.5 12.0 12.5 13.5 11.5 8.6 6.8 5.5
35 5.5 7.6 9.2 10.3 11.1 9.7 8.9 11.0 7.5 6.2 3.5
40 4.2 5.5 7.5 8.5 9.6 8.6 7.0 8.6 6.5 5.5 4.0
45 3.5 4.5 6.3 8.0 9.0 8.0 6.7 6.5 6.0 5.0 3.6
50 3 4.0 5.6 7.6 8.5 7.2 6.1 5.2 4.5 3.0 2.5





Evaluacin de Laboratorio N 6

Las preguntas medirn las siguientes competencias

Mapeo y Flujo.
Resolucin numrica de una EDO.
Construir una curva interpolante y de ajuste.
Cuadratura y cubicacin