![UNPRG
JOSÉ LUIS FARRO QUESQUÉN
PRACTICA DE MATLAB
1. Resolver:
2 6 5 4
6 7 8
6 3 7 9
x y z
x y z
x y z
Solución
>> A=[2 -6 5;6 -7 1;6 -3 7];
>> b=[4;8;9];
>> X=inv(A)*b
X =
1.2008
-0.0820
0.2213
2. Derivar y hallar sus gráficas
4 3
4 6 3 6y x x x
Solución
>> x=0:0.5:10;
>> y=4*x.^4 +6*x.^3 - 3*x +6;
>> dy=16*x.^3 + 18*x.^2 -3;
>> ddy=48*x.^2+18;
>> dddy=96*x;
>> plot(x,y,'r',x,dy,'k',x,ddy,'g',x,dddy,'b')
>> grid on
>> title('La función polinómica y sus derivadas')
>> legend('y','dy','ddy','dddy')
>> xlabel('Eje x')
>> ylabel('Eje y')](https://image.slidesharecdn.com/practicadematlab-180311141920/85/Practica-de-matlab-1-320.jpg)
![UNPRG
JOSÉ LUIS FARRO QUESQUÉN
3. Resolver
´ 2
'
' 4 4 5
(0) 1, (0) 2, (0) 3
x x y z
y x z
z x y z
x y z
Solución
inic='x(0)=1,y(0)=2,z(0)=3'
[x y z]=dsolve('Dx=x+2*y-z','Dy=x+z','Dz=4*x-4*y+5*z',inic,'t')
% Graficando la solución
t=linspace(0,0.5,25)
xx=eval(vectorize(x))
yy=eval(vectorize(y))
zz=eval(vectorize(z))
plot(t,xx,t,yy,t,zz)
grid
legend('x(t)','y(t)','z(t)')
title('Sistema de ecuaciones diferenciales')](https://image.slidesharecdn.com/practicadematlab-180311141920/85/Practica-de-matlab-2-320.jpg)
Practica de matlab
- 1. UNPRG JOSÉ LUIS FARRO QUESQUÉN PRACTICA DE MATLAB 1. Resolver: 2 6 5 4 6 7 8 6 3 7 9 x y z x y z x y z Solución >> A=[2 -6 5;6 -7 1;6 -3 7]; >> b=[4;8;9]; >> X=inv(A)*b X = 1.2008 -0.0820 0.2213 2. Derivar y hallar sus gráficas 4 3 4 6 3 6y x x x Solución >> x=0:0.5:10; >> y=4*x.^4 +6*x.^3 - 3*x +6; >> dy=16*x.^3 + 18*x.^2 -3; >> ddy=48*x.^2+18; >> dddy=96*x; >> plot(x,y,'r',x,dy,'k',x,ddy,'g',x,dddy,'b') >> grid on >> title('La función polinómica y sus derivadas') >> legend('y','dy','ddy','dddy') >> xlabel('Eje x') >> ylabel('Eje y')
- 2. UNPRG JOSÉ LUIS FARRO QUESQUÉN 3. Resolver ´ 2 ' ' 4 4 5 (0) 1, (0) 2, (0) 3 x x y z y x z z x y z x y z Solución inic='x(0)=1,y(0)=2,z(0)=3' [x y z]=dsolve('Dx=x+2*y-z','Dy=x+z','Dz=4*x-4*y+5*z',inic,'t') % Graficando la solución t=linspace(0,0.5,25) xx=eval(vectorize(x)) yy=eval(vectorize(y)) zz=eval(vectorize(z)) plot(t,xx,t,yy,t,zz) grid legend('x(t)','y(t)','z(t)') title('Sistema de ecuaciones diferenciales')
- 3. UNPRG JOSÉ LUIS FARRO QUESQUÉN