PRIMER EXAMEN PARCIAL (1/2016)
Programación y análisis numérico (PRQ 404)
En los problemas siguientes plantee el algoritmo, elabore el diagrama de flujo (o pseudocódigo) y
realice el programa en MATLAB:
Ejercicio 1:
El valor de e se calcula con la siguiente serie de Taylor, para x=1:
𝑥 𝑥2 𝑥3
𝑒𝑥 = 1 + + + + ⋯, −∞<𝑥 <∞
1! 2! 3!
a) Calcule el valor de e, determine el error e imprima estos valores.
b) Para cualquier valor de x, calcule el valor de 𝑒 𝑥
Ejercicio 2
Generar el n-ésimo valor de la serie Fibonacci. Tome en cuenta que la serie de Fibonacci es: 0, 1, 1,
2, 3,5, 8, … el 1er término es 0, el quinto es 3, ¿el n-ésimo será…?
Ejercicio 3
Realice un solo gráfico que muestre 4 cuadrantes de gráfico que contenga lo siguiente:
1) 𝑦1 = sin(𝑥), 2) 𝑦2 = cos(𝑥), 3) 𝑦3 = ln(𝑥), 4) 𝑦4 = 𝑒 −𝑥 . Realice los gráficos en el rango [0,2π].
Por indicaciones del docente en este ejercicio solo se realizó el codificado en MATLAB
%Realizar 4 graficas al mismo tiempo%
%Para todo el cálculo se generan valores de x%
x=-10:0.01:10;
%Para el primer gráfico%
y=sin(x);
subplot(2,2,1)
plot(x,y)
grid on
title 'Funcion sin'
hold on
%Para el segundo gráfico%
y=cos(x);
subplot(2,2,2)
plot(x,y,'r')
grid on
title 'Funcion cos'
hold on
%Para el tercer gráfico%
y=exp(-x);
subplot(2,2,3)
plot (x,y,'r')
grid on
title 'Función e^-x'
hold on
%Para el cuarto gráfico%
y=log(x);
subplot(2,2,4)
plot(x,y,'b--')
grid on
title 'Función ln(x)'
hold on
%Fin%
Ejercicio 4
Se ha sumergido el uso del metano, en un cilindro como un combustible de emergencia para el
sistema de calefacción de una planta que normalmente usa gas natural (compuesto por gran parte
por metano). Se ha de mantener en reserva una cantidad suficiente de cilindros de gas para
𝑘𝑐𝑎𝑙 𝑘𝑐𝑎𝑙
suministrar 25200 durante 24 horas, si el metano produce 97200 al quemarse y se
ℎ 𝑘𝑚𝑜𝑙
entrega en un cilindro de 56.62 litros a 204 atm y 21 °C. ¿Cuántos cilindros deben ser mantenidos
en reservas? Obtener predicciones basadas en cada una de las siguientes ecuaciones de estado,
ara el inciso b) utilice el método de Newton-Raphson o medianas.
a) Ecuación de Gas Ideal
b) Ecuación de Van der Waals
PRIMER EXAMEN PARCIAL (2/2016)
Programación y análisis numérico (PRQ 404)
Ejercicio 1
Generar un vector de n componentes e imprimir el valor mayor y menor.
%definiendo el tamaño del vector
n=input('tamaño del vector ');
%ciclo de generacion de valores
for i=1:1:n
a(i)=input('introdzca elementos del vector ');
end
%ciclo de repeticion del proceso
for i=1:1:n
%ciclo de ordenamiento
for j=1:1:n-1
if a(j+1)>a(j)
c=a(j);
a(j)=a(j+1);
a(j+1)=c;
end
end
end
%impresion de valores buscados
fprintf ('los valores buscados son %f (mayor) y %f (menor)\n',a(1),a(n))
Ejercicio 2
Se tiene un reactor en el que la constante de equilibrio en la reacción A → B se representa por:
1
(1 − 𝑥) ∗ (3 + 𝑥)2
𝑘𝑝 = 1 1
𝑥 ∗ (𝑥 + 1)2 ∗ 𝑝2
Donde x es la conversión del compuesto A a B, encuentre dicha concentración.
clc
clear all
%introduciendo el valor de aproximacion entre 0 y 1 tal que sea
convergente
xo=input('introduzca el valor inicial de x ');
%definiendo los valores constantes
kp=3.6;
p=5;
T=300;
%definiendo la funcion
f=@(x,p,kp) ((((1-x)*((3+x)^(1/2)))/((x*((x+1)^(1/2))*((p)^(1/2))))-kp));
%definiendo los limites y el paso
h=10^-6;
%error arbitrario para permitir el funcionamiento del programa
error=10^-3;
%ciclo de aproximacion al valor real
while error>10^-6
g=f(xo+h,p,kp);
df=((f(xo+h,p,kp)-f(xo,p,kp))/h);
x1=xo-(g/df);
error=abs(x1-xo);
fprintf ('x es: %f; f(x)es: %f; df es: %f ;el nuevo x es: %f ;el
error es: %f\n',xo,g,df,x1,error)
xo=x1;
end
%impresion del valor buscado
Ejercicio 3
Realice la integración numérica por el método del trapecio de:
𝑥 3 ∗ sin(𝑥) ∗ 𝑒 𝑥
𝑓(𝑥) = 3
𝑥 + sin(𝑥) + 𝑒 𝑥
clc
clear all
% definiendo la funcion
f=@(x) (((x^3)*(sin(x))*exp(x))/((x^3)+sin(x)+exp(x)));
%definiendo el numero de intervalos
n=input ('introduzca el numero de intervalos deseados ');
%definiendo los limites
a=0;
b=1;
h=(b-a)/n;
%primer ciclo de generacion de valores
for i=1:1:n;
x(i)=a+(i*h);
y(i)=f(x(i));
end
s=0;
%ciclo de generacion de sumatorias
for i=2:1:n-1;
s=s+y(i);
end
%integracion
I=(h/2)*(y(1)+(2*s)+y(n));
%impresion del resultado
fprintf ('el valor de la integral es %f\n',I)
plot(x,y)
