Hola Ivan.
%=====================
% dX/dt=-6*X X(0)=1
% dY/dt=6*X -8*Y Y(0)=0
% dZ/dt=8Y Z(0)=0
function S=ode23_caballero
opciones=odeset('RelTol',1e-5,'AbsTol',[1e-5 1e-5 1e-5]);
[T,V]=ode45(@ecua3,[0 1],[1 0 0 ],opciones);
S=[T V];
figure(gcf)
plot(T,V(:,1),T,V(:,2),T,V(:,3));
[x,y,z]=dsolve('Dx =-6*x', 'Dy = 6*x-8*y', 'Dz=8*y','x(0)=1', 'y(0)=0', 'z(0)=1');
legend(['x(t)=',char(x)],['y(t)=',char(y)],['z(t)=',char(z)])
grid
function dv=ecua3(t,v)
dv=zeros(3,1);
dv(1)=-6*v(1);
dv(2)=6*v(1)-8*v(2);
dv(3)=8*v(2);
%============================
EJECUCION.
>> S=ode23_caballero
S =
0 1.0000 0 0
0.0033 0.9802 0.0195 0.0003
0.0067 0.9608 0.0382 0.0010
0.0100 0.9418 0.0559 0.0023
0.0133 0.9231 0.0729 0.0040
0.0221 0.8761 0.1134 0.0105
0.0308 0.8314 0.1489 0.0197
0.0395 0.7890 0.1798 0.0312
0.0482 0.7488 0.2065 0.0447
0.0570 0.7103 0.2296 0.0601
0.0658 0.6738 0.2493 0.0769
0.0746 0.6392 0.2658 0.0951
0.0834 0.6063 0.2794 0.1143
0.0927 0.5733 0.2911 0.1356
0.1021 0.5420 0.3003 0.1577
0.1114 0.5125 0.3071 0.1804
0.1208 0.4845 0.3119 0.2036
0.1307 0.4566 0.3150 0.2284
0.1406 0.4302 0.3163 0.2535
0.1505 0.4053 0.3161 0.2786
0.1604 0.3819 0.3145 0.3036
0.1710 0.3584 0.3115 0.3301
0.1816 0.3364 0.3073 0.3563
0.1921 0.3157 0.3022 0.3820
0.2027 0.2963 0.2963 0.4074
0.2140 0.2769 0.2892 0.4339
0.2254 0.2587 0.2816 0.4598
0.2367 0.2417 0.2734 0.4849
0.2480 0.2258 0.2649 0.5093
0.2602 0.2099 0.2555 0.5347
0.2724 0.1951 0.2458 0.5591
0.2846 0.1813 0.2361 0.5826
0.2968 0.1685 0.2263 0.6052
0.3100 0.1557 0.2158 0.6285
0.3232 0.1438 0.2054 0.6508
0.3364 0.1329 0.1952 0.6720
0.3496 0.1227 0.1852 0.6921
0.3641 0.1126 0.1746 0.7128
0.3785 0.1032 0.1644 0.7324
0.3929 0.0947 0.1546 0.7508
0.4073 0.0868 0.1451 0.7681
0.4232 0.0789 0.1352 0.7859
0.4391 0.0717 0.1258 0.8025
0.4550 0.0652 0.1169 0.8179
0.4709 0.0593 0.1085 0.8322
0.4887 0.0533 0.0997 0.8470
0.5064 0.0479 0.0915 0.8606
0.5242 0.0431 0.0839 0.8730
0.5419 0.0387 0.0769 0.8844
0.5620 0.0343 0.0695 0.8962
0.5822 0.0304 0.0627 0.9068
0.6023 0.0269 0.0566 0.9165
0.6225 0.0239 0.0510 0.9251
0.6459 0.0207 0.0451 0.9341
0.6693 0.0180 0.0399 0.9421
0.6928 0.0157 0.0352 0.9491
0.7162 0.0136 0.0311 0.9553
0.7412 0.0117 0.0272 0.9611
0.7662 0.0101 0.0237 0.9662
0.7912 0.0087 0.0207 0.9706
0.8162 0.0075 0.0180 0.9745
0.8412 0.0064 0.0157 0.9779
0.8662 0.0055 0.0137 0.9808
0.8912 0.0048 0.0119 0.9834
0.9162 0.0041 0.0103 0.9856
0.9372 0.0036 0.0092 0.9872
0.9581 0.0032 0.0082 0.9887
0.9791 0.0028 0.0072 0.9899
1.0000 0.0025 0.0064 0.9911
y tambien te va a salir una grafica.
saludos
JOSE JEREMIAS CABALLERO
ASESOR DE PROYECTOS CON MATLAB
PROFESOR DE METODOS NUMERICOS CON MATLAB
PROGRAMADOR MATLAB
[email protected]