programacion informe 1

FACULTAD DE CIENCIAS Químicas
ESCUELA DE Ingeniería AMBIENTAL

Materia:

Programación y métodos numéricos

Tema:
Ejercicios (capítulo 1, 2 y 7)







CAPÍTULO 1
Página 32
(1.1) Ejecute los siguientes enunciados:
>>a = [1 2 3; 4 5 6]'
a =
1 4
2 5
3 6
>> b = [9; 7; 5; 3; 1]
b =
9
7
5
3
1
>> c = b (2:4)
c =7
5
3
>>d = b (4:-1:1)
d =
3
5
7
9
>> e=sort (b)
e =
1
3
5
7
9
>> f= [3, b’]
f =
3 9 7 5 3 1
(1.2) Cree un archivo M de función, fun_es(x), que calcule la siguiente función:
y = 0.5ex/3 – x2senx
El argumento debe aceptar tanto un escalar como un vector. Pruebe su función tecleando en MATLAB.En el editor
function y=fun_es(x)

y=(0.5*exp(x/3))-(x.^2)*sin(x)

En el comand window

>> fun_es(3)
y =
0.089060841690718
ans =
0.089060841690718
En el editor
function y=fun_es(x)
y=0.5*exp(x./3)- x.^2.*sin(x)
end

En el comand window



(1.3)Repita la tarea del problema 1.2 para la función:

Denote la función por fun_lg(x)
(a)
En el editor
function y =fun_log(x)
f=input('ingrese la funcion:','s');
x=input('ingrese valores x > 0')
y=sin(x)*log(1+x)-x.^2;
En el comand window
ingrese la funcion:fun_log(3)
ingrese valores x > 03
x =
3
ans =
-8.8044
(b)
En el editor
function y = fun_log(x)
f=input('ingrese la funcion:','s');
x=input('ingrese valores x > 0')
y=sin(x).*log(1+x)-x.^2;
en el comand window
ingrese lafuncion:fun_log([1 2 3])
ingrese valores x > 0([1 2 3])

x =
1 2 3


ans =
-0.4167 -3.0010 -8.8044
(1.4)
(a) Escriba un archivo M de función que calcule la solución de

Su forma general es r_cuad(a,b,c) donde a, b y c pueden ser vectores.(b) pruebe la función con a=3, b=1, c=1, r_cuad(a,b,c). (c) Pruebe la función con a= [3 1 2], b=[1 -4 9], c=[1 3 -5].
En el editor
functionr_cuad(a,b,c)
a=input('ingrese el valor de a');
b=input('ingrese el valor de b');
c=input('ingrese el valor de c');
x1=(-b+sqrt(b.^2-4*a*c))./(2*a)
x2=(-b-sqrt(b.^2-4*a*c))./(2*a)

En el command window

ingrese el valor de a3
ingrese el valor de b1
ingrese el valor de c1
x1 =
-0.1667 + 0.5528i
x2 =
-0.1667 - 0.5528i
Segunda parte
En el editor
function r_cuad(a,b,c)a=input('ingrese el valor de a');
b=input('ingrese el valor de b');
c=input('ingrese el valor de c');
x1=(-b+sqrt(b.^2-4.*a.*c))./(2.*a)
x2=(-b-sqrt(b.^2-4.*a.*c))./(2.*a)

En el command window

ingrese el valor de a[3 1 2]
ingrese el valor de b[1 -4 9]
ingrese el valor de c[1 3 -5]
x1 =
-0.1667 + 0.5528i 3.0000 0.5000
x2 =
-0.1667 - 0.5528i 1.0000-5.0000

(1.9) se tiene un vector:

Escriba un guion que imprima el contenido del vector empleando el comando fprintf en un ciclo tal que el aspecto de la salida sea:
El vector A es

Termino la impresión.
En el editor
A=[1 2 3 4 5 6 7 8 9 0]
fprintf('el vector A es\n');
fprintf('[')
fprintf('%1.0f, ',A)
fprintf(']\n');
fprintf('terminó la impresion')

En elcomand window

A =
1 2 3 4 5 6 7 8 9 0
el vector A es
[1, 2, 3, 4, 5, 6, 7, 8, 9, 0, ]
terminó la impresion>>



CAPITULO 2
Página 92
En las siguientes tareas de trazado de figuras, imprima leyendas de ejes y leyendas para explicar el significado de cada curva si hay másde una. También deberá imprimir el titulo de la figura en cada grafica que elabore. Utilice labelx, labely, title, text, gtext.
(2.1)Grafique las siguientes funciones en el dominio que se indica:
(a)
En el comand window
>> x= (0:0.05:4*pi)';
>> y=sin(x)/(1+cos(x));
>> plot(x,y)
>> xlabel('x'); ylabel('y')
>> title('GRAFICA 1')
>> grid
GRÁFICO

(b)
(c)
En el comand window
>> x=...