CAPITULO IV

INTRODUCCIÓN A MATLAB

4.1 OBJETIVOS

a) Conocer y aprender el entorno de trabajo de MATLAB.

b) Aprender a manejar las variables y asignarles valores.
c) Manejar los operadores matemáticos y su prioridad, así como la aplicación de
algunas funciones de biblioteca.
d) Aprender a trabajar con varios tipos de formatos.
e) Realizar programas sencillos utilizando archivos m

4.2 BASES TEÓRICAS

4.2.1 Introducción
Matlab (MATrix LABoratory. Es una herramienta importante para cualquier tarea que
requiera cálculos matriciales, ya sea que involucren ecuaciones, sistemas de ecuaciones,
análisis de datos, etc. y/o la visualización grafica de los mismos.
Matlab trabaja con un solo tipo de objetos: matrices, o sea que todas las variables
representan matrices, así que los escalares son matrices de 1x1 y los vectores son
matrices de una fila (1x n), o de una columna (nx1).

4.2.2 Elementos del entorno gráfico de Matlab.

Hay dos modos de trabajo con MATLAB:

Trabajo interactivo (Workspace), donde el usuario realiza una consulta (escribe una
operación) y Matlab la ejecuta cuando el usuario presiona <Enter>.

Trabajo programado (M-files), donde el usuario genera uno o varios archivos con
conjuntos de instrucciones MATLAB, que se pueden ejecutar repetidas veces (con distintos
datos) desde el modo interactivo. De esta forma, el usuario puede incrementar las
funciones disponibles en MATLAB, añadiendo las suyas propias.

MATLAB se puede arrancar como cualquier otra aplicación de Windows, haciendo click
dos veces en el icono correspondiente en el escritorio o por medio del menú Inicio. Al
arrancar MATLAB se abre una ventana. Ésta es la vista que se obtiene eligiendo la opción
Home Layout/Default, aunque esta configuración puede ser cambiada fácilmente.

1
Editor: Ruta: Donde MATLAB tiene Menú Layout: Donde cambias la Help: (F1) Manual de Ayuda para utilizar
Donde almacenados tus algoritmos. vista del escritorio en MATLAB todos los recursos que tiene MATLAB
escribes tus
algoritmos.

Workspace: Espacio
Current Directory: Command Window: Espacio donde donde permite
Donde MATLAB puedes escribir comandos como: observar las
tiene almacenados cálculos (suma, resta, etc.), desarrollar variables definidas
tus algoritmos. gráficas (figuras) llamar funciones de en el comando o
ayuda. mediante un
algoritmo.

Command History: Espacio donde

puedes observar los comandos hechos
anteriormente.

2
4.2.3 Operaciones aritméticas básicas.

Operación Símbolo Ejemplo

Suma + 6+2=8

Diferencia - 6-2=4

Multiplicación * 6*2=12

División / 6/2=3

Elevado a la potencia ^ 6^2=36

La prioridad de las operaciones:

1. Exponenciación
2. Multiplicación y división
3. Suma y Resta

Paréntesis: Para paréntesis anidados, el más interno es el primero que se

ejecuta.

4.2.4 Variables
Los datos se pueden almacenar en algún lugar de la memoria, asignando
valores a nombres de variables

1. Los nombres de las variables deben de comenzar con una letra y pueden tener
hasta 31 caracteres, que pueden ser números, letras, guión (-), subrrayado(_)
2. La variable a es distinta de la variable A, MATLAB diferencia entre mayúsculas
y minúsculas
3. Los nombres de las variables deben tener significado, primero se utilizan letras
como a, x, y, pero cuando los programas son más complejos se llamarán
posicion, velocidad, etc.
4. Los nombres de las variables no pueden coincidir con las palabras reservadas
por MATLAB: sin, sqrt exp, etc
5. No están permitidos espacios entre caracteres

3
4.2.4.1 Variables Predefinidas

Representa el valor calculado de una expresión y que se puede reutilizar

ans
en el siguiente cálculo. >> 2*5+1.5 ans=11.5000

pi Representa el número pi

eps Es la más pequeña diferencia entre dos números. Precisión de los

números en coma-flotante

i, j √−1 , o el número complejo 0+1.0000i

Representa el infinito, por ejemplo cuando se divide un número entre cero,

inf
se muestra ∞. >>a=1/0 Warning: Divide by zero a =inf
No es un número. Resultado de una expresión indefinida, por ejemplo al
NaN
dividir cero entre cero. >>b=0/0 Warning: Divide by zero b=NaN

4.2.5. Comandos de ayuda

Comando Acción
>> prompt Indica que MatLab está listo para aceptar órdenes.
Para poder visualizar cuales variables han sido ya
Who
introducidas.
Nombre de la variable Para conocer el contenido de una variable.
clc Limpia la pantalla de comandos.
clear Para borrar borran todas las variables existentes.
clear nombrevariable Para borrar una variable de la memoria.
Para obtener ayuda. Nos despliega todo un listado de
help
temas.
Si se teclea help seguido de un comando, función o
help
archivo de matlab, aparecerá una descripción del
nombre_comando
comando.

Se desplazaran los comandos anteriores. También Se

pueden utilizar las flechas para localizar un comando y
, , y
modificarlo y al presionar la tecla Intro se ejecutará el
comando modificado.
Si se inicia una línea con el símbolo % se interpretará
%
como un comentario. Ejemplo: % método simple.
Se utiliza cuando se escribe una instrucción en matlab y
; (punto y coma)
no se desea ver los resultados desplegados.
Se utiliza para detener el proceso de cualquier tipo de
Control + C
cálculo, gráfico o impresión.

4
4.2.6 Archivos m
Matlab puede ejecutar una sucesión de instrucciones almacenadas en archivos.
Estos archivos se denomina archivos m debido a que su extensión por default es
m.

Los archivos m pueden ser de comandos o de funciones.

4.2.6.1 Archivos m de comandos

Estos archivos consisten en una sucesión de instrucciones de Matlab y son
guardados como un simple archivo de texto (solo que su extensión es m y no txt).
Pueden ser usados para introducir gran cantidad de datos o cuando se desea
cambiar el valor de una o más variables y reevaluar una serie de órdenes.

Para iniciar un archivo m se hace click en new del menú Home y se selecciona
New Script.
Aparecerá una nueva ventana para la edición de texto donde se teclea la lista de
comandos de Matlab. Para guardarlo como archivo m en la unidad de
almacenamiento se hace click en Save/Save as del menú Editor y se escribe el
nombre del archivo, la extensión m se genera automáticamente.

Escribiendo simplemente el nombre del archivo en la ventana command window

se ejecutaran, línea por línea las órdenes contenidas en el archivo.

2 A=26
3 B=18
4 C=54
5 D=A+C*B
Botón para
abrir un nuevo Botón de ejecutar
4
4 Editor.
4
54
5
5
5

5
4.2.6.2 Archivos m de función
Estos archivos permiten que el usuario pueda definir sus propias funciones. Son
similares a los de comandos, son archivos de texto y tienen también una extensión
m.

La primera línea activa tiene el siguiente formato obligatorio:

function[v1,v2,…,vn]=nombre_archivo(a1,a2,…,ak)

donde v1,v2,…vn son las variables de salida y a1,a2,…ak son los argumentos o
variables de entrada, que serán asignadas en el mismo orden al convocar la
función.

4.2.7. Algunos formatos de variables

Comando Función Ejemplo

Es el formato por defecto;
format short 2/3 = 0.6667
punto fijo y 4 decimales
format short e Notación científica con 4 2/3 = 6.6667e-001
decimales.
format long punto fijo y 14 decimales 2/3 = 0.66666666666667
format long e notación científica con 15 2/3 =6.666666666666666e-001
decimales
format rat Aproximaciones al número 2/3 = 2/3
racional más cercano.
format bank agrega 2 decimales 290/7=41.43

Representación de los números en MATLAB

Número Mantisa-exponente MATLAB

2412.6 2.4126×103 2.4126e3

0.00002 2×10-5 2e-5

20000 2×104 2e4

6
4.2.8. Funciones Elementales
MATLAB dispone de numerosas funciones predefinidas, algunas de ellas son las
siguientes:

Comando Función

abs(x) Valor absoluto de x

sqrt(x) Raíz cuadrada de x

nthroot(x,n) Raíz n de x. Si n=3 sería la raíz cúbica de x

sign(x) Devuelve -1 si x es menor que 0, 0 si x es igual a 0, y 1 en el

otro caso

rem(x,y) Devuelve el resto de la división x/y Ejemplo: rem(26,4)

devuelve 2

exp(x) Calcula ex

log(x) Calcula lnx, logaritmo natural

log10(x) Calcula log10x, logaritmo en base 10 de x

sin(x), Calcula el seno de x, x en radianes (o grados)

sind(x)

cos(x) Calcula el coseno de x, x en radianes (o grados)

cosd(x)

tan(x) Calcula la tangente de x, x en radianes (o grados)

tand(x)

asin(x) Calcula el arco seno, donde x es un número entre -1 a 1. La

asind(x) función devuelve el ángulo en radianes entre -π/2 y +π/2 (o
grados)

acos(x) Calcula el arco coseno, donde x es un número comprendido

acosd(x) entre -1 a 1. La función devuelve el ángulo en radianes entre 0
y π (o grados)

atan(x) Calcula el arco tangente. La función devuelve el ángulo en

7
atand(x) radianes entre -π/2 y +π/2 (o grados)

atan2(y,x) Calcula el arco tangente de y/x. Devuelve el ángulo en radianes

entre -π y +π, dependiendo del signo de x e y

sinh(x) 𝑒 𝑥 −𝑒 −𝑥
Calcula 2

cosh(x) 𝑒 𝑥 +𝑒 −𝑥
Calcula 2

tanh(x) 𝑒 𝑥 −𝑒 −𝑥
Calcula 𝑒 𝑥 +𝑒 −𝑥

asinh(x) Calcula 𝑙𝑛(𝑥 + √𝑥 2 + 1)

acosh(x) Calcula 𝑙𝑛(𝑥 + √𝑥 2 − 1)

atanh(x) 1+𝑥
Calcula 𝑙𝑛√1−𝑥 𝑝𝑎𝑟𝑎 |𝑥| ≤ 1

Funciones de redondeo:
La función round aproxima un número decimal al entero más próximo.
>> round(2.52) >> round(2.49) >> round(2.5)
ans = 3 ans = 2 ans = 3

La función ceil redondea hacia arriba, la función floor redondea hacia abajo
>> floor(2.51) >> ceil(2.4)
ans = 2 ans = 3

La función fix elimina la parte decimal

>> fix(2.49) >> fix(-2.4)
ans = 2 ans = -2

8
4.3 PRE-LABORATORIO

Operaciones básicas
1. Elabore un archivo m ((Nombre del archivo: Pre4_1_ccd )
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Pasos: hacer click en New del menú File y seleccione M-file. A continuación
aparecerá una nueva ventana para la edición del texto donde teclearás las
siguientes instrucciones:

Resultados:
A = 31
B =1.3554
C =14.6931
D =0.0000e+00 + 1.1778e+03i
E =8.7304e+05 + 8.9210e+05i
F =-1.3333
G =4.6416e-08

De acuerdo a como vayas escribiendo cada ecuación ejecuta el programa y verifica

que los resultados sean los correctos.

9
2. Elabore un archivo m ((Nombre del archivo: Pre4_2_ccd)
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Resultados:

A= 12
B= 0.0000 + 2.0000i
C= 9
D= 0.5000
E= 0.5000
F= 1.0472
G= 60.0000
H= 2
I= 3.141592653589793
J= 3.1416
K= 50.1429
L= 50.142857142857146
M= 5.0143e+01
N= 5.014285714285715e+01
O= 50.14
P= 2/3
Q= 7/6
R= 1/3
S= 5/6
>>

Nota: con format short, Volvemos al formato numérico de origen del Matlab

De acuerdo a como vayas escribiendo cada ecuación ejecuta el programa y verifica

que los resultados sean los correctos.

10
3. Elabore un archivo m (Nombre del archivo: Pre4_3_ccd )

Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

11
4. Elabore un archivo m (Nombre del archivo: Pre4_4_ccd )
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Si un automóvil está viajando con una rapidez de 28 m/s, ¿el conductor está
excediendo el límite de velocidad de 55.0 millas/h?

Solución:
Para resolver este problema, los metros deben convertirse a millas y los segundos
a horas, aplicando factores de conversión.
1 Milla= 1609 m.
1 min = 60 seg.
1 hr = 60 min.

Conversión de m/s a millas/s

𝑚 𝑚 1.00 𝑚𝑖𝑙𝑙𝑎𝑠
28.0 = (28.0 ) ( ) = 1.74 𝑋 10−2 𝑚𝑖𝑙𝑙𝑎𝑠/𝑠
𝑠 𝑠 1609 𝑚

Conversión de segundos a horas

𝑚𝑖𝑙𝑙𝑎𝑠 𝑚𝑖𝑙𝑙𝑎𝑠 𝑠 𝑚𝑖𝑛 𝑚𝑖𝑙𝑙𝑎𝑠
1.74 𝑋 10−2 = (1.74 𝑋 10−2 ) (60.0 ) (60.0 ) = 62.6
𝑠 𝑠 𝑚𝑖𝑛 ℎ ℎ
Código en Matlab:

12
4.4 DESARROLLO DE LA PRÁCTICA

1. Elabore un archivo m (Nombre del archivo: Dprac4_1_ccd )

Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Calcular la 𝜌 (densidad) de un disco cuyo diámetro es de 2.5 𝑥1010 picometros,

con un espesor o altura de 8 𝑥106 𝑛𝑚, el disco tiene un peso de 34.5 𝑔𝑟.

Datos:
Diámetro=2.5x1010 picometros 1 picometro=1x10-12 metros
Altura=8x106 nanómetro 1 nanómetro=1x10-9 metros
Masa=34.5 gramos 1 metro=100 centímetros

Formulas:
𝑚𝑎𝑠𝑎 𝑔𝑟
𝜌𝑑𝑖𝑠𝑐𝑜 = (𝑐𝑚3 ) 𝑣 = 𝜋𝑟 2 ℎ
𝑣𝑜𝑙𝑢𝑚𝑒𝑛
Convertir el diámetro a centímetros:
𝑚 𝑐𝑚
𝐷 = (2.5𝑥1010 𝑝𝑚) (1𝑥10−12 ) (100 ) = 2.5𝑐𝑚
𝑝𝑚 𝑚
Convertir la altura a centímetros:
𝑚 𝑐𝑚
ℎ = (8𝑥106 𝑛𝑚) (1𝑥10−9 ) (100 ) = 0.8𝑐𝑚
𝑛𝑚 𝑚
Calcular el volumen:
2
2.5
𝑣 = 𝜋 ( 𝑐𝑚) (0.8𝑐𝑚) = 3.927𝑐𝑚3
2
Calcular la densidad:
34.5 𝑔𝑟
𝜌= = 8.7854 𝑔𝑟/𝑐𝑚3
3.927𝑐𝑚3

13
 2. Elabore un archivo m (Nombre del archivo: Dprac4_2_ccd )
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

𝑦𝑎𝑟𝑑𝑎𝑠
Convertir 1000𝑐𝑚/𝑠𝑒𝑔3 a ℎ 𝑆𝑒𝑚𝑎𝑛𝑎 𝑀𝑖𝑛
. Haciendo la conversión por pasos ejemplo:
cm a yardas, seg a hrs, seg a semana, seg a min. Indicando con un texto cada
conversión por pasos. Ir tomando el valor anterior para cada operación.

Datos:
c=1000 cm/s3 1 hr=60 min
1 yarda=91.44 cm 1 día=24 hr
1min=60 seg 1sem=7 días

𝑐𝑚 1 𝑦𝑎𝑟𝑑𝑎 𝑦𝑎𝑟𝑑𝑎𝑠
𝑐 = (1000 ) ( ) = 10.9361
𝑠 3 91.44 𝑐𝑚 𝑠3

𝑦𝑎𝑟𝑑𝑎𝑠 𝑠 𝑚𝑖𝑛 𝑦𝑎𝑟𝑑𝑎𝑠

𝑐 = (10.9361 ) (60 ) (60 ) = 3.9370𝑒 + 04
𝑠3 𝑚𝑖𝑛 ℎ𝑟 𝑠 2 ℎ𝑟

𝑦𝑎𝑟𝑑𝑎𝑠 𝑠 𝑚𝑖𝑛 ℎ𝑟 𝑑𝑖𝑎

𝑐 = (3.9370𝑒 + 04 ) (60 ) (60 ) (24 ) (7 )
𝑠 2 ℎ𝑟 𝑚𝑖𝑛 ℎ𝑟 𝑑𝑖𝑎 𝑠𝑒𝑚

𝑦𝑎𝑟𝑑𝑎𝑠
𝑐 = 2.3811𝑒 + 10
𝑠 ℎ𝑟 𝑠𝑒𝑚

𝑦𝑎𝑟𝑑𝑎𝑠 𝑠 𝑦𝑎𝑟𝑑𝑎𝑠
𝑐 = (2.3811𝑒 + 10 ) (60 ) = 1.4287𝑒 + 12
𝑠 ℎ𝑟 𝑠𝑒𝑚 𝑚𝑖𝑛 ℎ𝑟 𝑠𝑒𝑚 𝑚𝑖𝑛

14
 3. Elabore un archivo m (Nombre del archivo: Dprac4_3_ccd )
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

El Sol irradia 385𝑥1024 𝐽/𝑠 de energía, los cuales se generan mediante reacciones
nucleares que convierten la materia en energía. Realizar la conversión de J/s a J,
con la ecuación de Einstein determine cuánta materia se debe convertir en
energía para producir esa cantidad de radiación en un día. 𝐸 = 𝑚𝑐 2 , Rapidez de la
luz c = 3𝑥108 𝑚/𝑠.

Datos:
E=385x1024 j/s 1 hr=60 min
c=3x108 m/s 1 dia =24 hr
1 min=60 seg

formula: E=mc2

385𝑥1024 𝐽 60𝑠 60𝑚𝑖𝑛 24ℎ𝑟

𝐸=( )( )( )( ) (1𝑑𝑖𝑎) = 3.3264 x 1031 𝐽
𝑠 𝑚𝑖𝑛 ℎ𝑟 𝑑𝑖𝑎

Despejando m de la ecuación de energía:

𝐸
𝑚=
𝐶2
Y sustituyendo valores:
3.3264 x 1031 𝐽 3.6960 x 1014 𝐽
𝑚= 2 =
3 x 108 𝑚 𝑚2 /𝑠 2
( )
𝑠

15
 4. Elabore un archivo m (Nombre del archivo: Dprac4_4_ccd )
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Los túneles de viento juegan un importante papel en la comprensión del

comportamiento de las aeronaves de alto rendimiento. Con la finalidad de
interpretar los datos del túnel de viento, los ingenieros requieren entender cómo
se comportan los gases. La ecuación básica que describe las propiedades de los
gases es la ley del gas ideal, Dónde:

Formula: PV = nRT
Datos: P = 100 kpa
V =1000 𝑚3
n = número de kmoles de gas en la muestra,
R = constante de gas ideal, 8.314 kpa m3/kmol K,
T = 300 K

se sabe que el número de kmoles de gas es igual a la masa del gas dividida por la
masa molar (también conocida como peso molecular)

Formula: n = m/MW
Datos: m = masa en kg
MW =29 kg/kmol. (Peso molecular)

Encuentre la masa m del aire en el túnel de viento.

𝑃𝑉 (100𝑘𝑝𝑎)(1000𝑚3 )
𝑃𝑉 = 𝑛𝑅𝑇 𝑑𝑒𝑠𝑝𝑒𝑗𝑎𝑛𝑑𝑜 𝑛 = 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑦𝑒𝑛𝑑𝑜 𝑛 =
𝑅𝑇 8.314𝑘𝑝𝑎 𝑚3
( ) (300𝐾)
𝑘𝑚𝑜𝑙 𝐾
𝒏 = 𝟒𝟎. 𝟎𝟗𝟑𝟎 𝒌𝒑𝒂
𝑚 𝑘𝑔
𝑛= 𝑑𝑒𝑠𝑝𝑒𝑗𝑎𝑛𝑑𝑜 𝑚 = (𝑛)(𝑀𝑊) 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑦𝑒𝑛𝑑𝑜 𝑚 = (40.0930𝑘𝑝𝑎) (29 )
𝑀𝑊 𝑘𝑝𝑎
𝒎 = 𝟏. 𝟏𝟔𝟐𝟕𝒆 + 𝟎𝟑 𝒌𝒈

16
4.5. POST-LABORATORIO
1. Elabore un archivo m (Nombre del archivo: Post4_1_ccd )
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Calcular el volumen de un cilindro. La fórmula para el cálculo es 𝒗 = 𝝅𝒓𝟐 ℎ. Hacer

el cálculo para h =1, h=5, h=12. Donde r=3.

2. Elabore un archivo m (Nombre del archivo: Post4_2_ccd )

Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Las poblaciones de bacterias tienden a expandirse exponencialmente. Esto es

donde
𝑷 = 𝑷𝟎 𝒆𝒓𝒕

P = población actual,
P0 = población original,
r = tarifa de crecimiento continua, expresado como fracción
t = tiempo.
Si originalmente se tienen 100 Bacterias que se reproducen a una tasa de
crecimiento constante de 90% (r = 0.9) por año, encuentre cuántas Bacterias
tendrá al final de 10 años.

17
 3. Elabore un archivo m (Nombre del archivo: Post4_3_ccd )
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

La predicción del clima es un proceso complicado, incluso con los mejores datos.
Una ecuación que usan los meteorólogos es la ecuación Clausius-Clapeyron, esta
ecuación se emplea para determinar la relación entre presión de vapor de agua de
saturación y la temperatura atmosférica. La presión de vapor de agua de
saturación se puede usar para calcular la humedad relativa, un componente
importante de la predicción del clima, cuando se conoce la verdadera presión
parcial del agua en el aire.

La ecuación Clausius-Clapeyron es
𝑝° ∆𝐻𝑣 1 1
ln( )=( )∗( − )
6.11 𝑅𝑎𝑖𝑟 273 𝑇
Dónde:
P0 = presión de vapor de saturación para el agua, en mbar, a temperatura T,
∆𝐻𝑣 = calor latente de vaporización para el agua, 2.453 𝑥106 J/kg,
Rair = constante de gas para el aire húmedo, 461 J/kg
T = temperatura en grados kelvin (K).

Es raro que las temperaturas en la superficie de la Tierra sean menores que -60 °F
o mayores que 120 °F.
Use la ecuación Clausius-Clapeyron para encontrar la presión de vapor de
saturación para temperaturas en este rango (T=-60°F y T= 120°F).

Nota la temperatura debe estar en K.

5(𝐹 − 32)
𝐾= + 273.15
9

18
 4. Elabore un archivo m (Nombre del archivo: Post4_4_ccd )
Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Suponga que queremos conocer el porcentaje de cada nutriente en una porción de

un alimento cuyo peso P se da en gr y cuya composición es: Proteína (Prot),
Carbohidrato (CH) y Grasa (Gr).

Usando los valores:

Prot=30
CH=45
Gr=15
P=115

Y las ecuaciones:
porcPro = (Prot/P)*100;
porcCar = (CH/P)*100
porcGra = (Gr/P)*100

Calcular y mostrar los resultados obtenidos.

5. Elabore un archivo m (Nombre del archivo: Post4_5_ccd )

Nota: ccd, es la clave que tienes asignada en el centro de cómputo docente.

Utilice los datos del problema anterior para calcular el porcentaje calórico
correspondiente a cada nutriente (porCalP, porCalCH y porCalGr) y el porcentaje
calórico total (Totcal).

Usando las equivalencias a calorías de cada nutriente como sigue:

calPro =4
calCH =4
calGr =9

y las ecuaciones:
CalCH= CH*calCH
CalGr = Gr*calGr
CalPro = Prot*calPro
Totcal =CalCH+CalGr+CalPro
porCalP = (CalPro/Totcal)*100
porCalCH = (CalCH/Totcal)*100
porCalGr = (CalGr/Totcal)*100.

19