programacion

UNIVERSIDAD DE CUENCA FACULTAD DE CIENCIAS QUÍMICAS ESCUELA DE INGENIERÍA AMBIENTAL
Realizado por: Jenny Patiño
Programación y métodos numéricos
1.1) Ejecute los siguientes enunciados
a= [1 2 3; 4 5 6]'
b= [9;7;5;3;1]
c= b(2:4)
d=b(4:-1:1)
e= sort (b)
f= [3,b']
>>
a = 1 4
2 5
3 6


b =
9
7
5
3
1
c =
7
5
3
d =
3
5
7
9

e =
1
3
5
7
9


f =
3 9 7 5 3 1


1.2) Cree un archive m de funsión, fun_es(x), que calculi la siguiente función:
Y=0.5*EXP(X/3)-X.^2*SIN(X)
function y=fun_es(x)
y=(0.5.*(exp(x./3)))-((x.^2).*sin(x))
y=fun_es(3)
y =0.0891
y=fun_es([1 2 3])
y = -0.1437 -2.6633 0.0891
1.3)Repita la tarea del problema 1.2 para la función:
Y=sen ( x ) log(x+1)-x^2,x>0
function y=fun_lg(x)
y=(sin(x).*log(1+x))-(x.^2);
>>
y=fun_lg(3)
y = -8.8044
y=fun_lg([1 2 3])
y = -0.4167 -3.0010 -8.8044

1.5) Se supone que el lectorya creo fun_es y fun_lg de los problemas 1.2 ; 1.3. Ahora cree una función f_es(x) q:
Que pregunte el nombre de la función que quiere evaluar
Permita al usuario teclear el nombre de la función evalué la función con feval y devuelva los valores funcionales y pruebe su f_es calculando fun_es(3) y fun_lg(3)

function y=f_es(x)
y=input('Ingrese el nombre de la función que desea evaluar= ', ‘s’);if a==(fun_es)
y=feval(fun_es,3)
elseif a==(fun_lg);
y=feval(fun_lg,3)
else
return
end
end

1.9) se tiene un vector: A=[1 2 3 4 5 6 7 8 9 0], escriba un guion que imprima el contenido del vector empleando el comando fprint en un ciclo tal que el aspecto de la salida sea:
El vector a es:
[1, 2, 3, 4, 5, 6, 7, 8, 9, 0]
Termino la impresión.
En la salida, coloque una coma y dosespacios en blanco entre dos números. Los elementos de a de

A=[1 2 3 4 5 6 7 8 9 0];
fprintf('El vector A es \n')
fprintf(' ')
fprintf('%3.0f,',A)
fprintf(']')
fprintf('\n')
fprintf('Terminó la impresión \n ')
El vector A es
[ 1, 2, 3, 4, 5, 6, 7, 8, 9, 0,]
Terminó la impresión
1.12)Cree un archivo m de función llamado fun_xb que evalué la siguienteserie:

Los valores de x y n se pasan a la funsion mediante argumentos. Prueba la función comprando con cálculos manuales para x=1 y n=4. La serie en cuestión es una expansión mclaurin truncada de log(1+x) y converge para . Sabiendo esto, pruebe su función con valores de x selectos como x=-0.5y 0.5 con n=1, 2, 3, 5, 10, 20, 50 y compare el resultado con log(1+x).
x=input('Ingrese el valor dex= ');
for a=1:n
sum=0;
Q(a)=fun_xa(x);
sum=sum+Q(a);
end
y=sum;
disp(y);
function y=fun_xb(x)
a=input('Ingrese el valor de n= ');
for a=1:2:n
for b=2:2:n
y=((x.^(n))./a)-((x.^(b))./b)
end
end
>>Ingrese el valor de x= 1
Ingrese el valor de n= 4
>> y= 0.6931
(7.2) Encuentre en forma aproximada todas las soluciones positivas de las siguientes ecuacionesutilizando el método gráfico:
a) Tan(x)-x+1=0, 0 %Tan(x)-x+1=0
>>x=-4:0.1:4;
>>y=tan(x)-x+1;
>>plot(x,y)
>>grid
>>xlabel('x');ylabel('y=tan(x)-x+1')

Raíz más o menos en x=-1.2





b)Sen(x)-0.3=0, x>0
>> x=-4:0.1:4;
>>y=sin(x)-0.3*exp(x);
>>plot(x,y)
>>grid
>> xlabel('x'); ylabel('y=sin(x)-0.3*exp(x)')










raíz más o menos en x=-3.05, x=0.45, x=1c) 0.1-5-x+4+=0
>> %0.1*x^3-5*x^2-x+4+exp(-x)=0
>> x=-4:0.1:4;
>> y=0.1.*x.^3-5.*x.^2-x+4+exp(-x);
>> plot(x,y)
>> grid
>> xlabel('x');ylabel('y=0.1.*x.^3-5.*x.^2-x+4+exp(-x)')

Raíz más o menos en x=-1.4, x=0.9








d)ln(x)-0.2-1=0
>> %log10(x)-0.2.*x.^2+1=0
>> x=-4:0.1:4;
>> y=log10(x)-0.2.*x.^2+1;
>> plot(x,y)
Warning: Imaginary parts of complex X and/or Y...