Don't like ads? PRO users don't see any ads ;-)
Guest

TAREA EDO ode45

By: a guest on May 2nd, 2012  |  syntax: MatLab  |  size: 0.26 KB  |  hits: 95  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. I=[0,60]; % Intervalo
  2. N=1000; % Divisiones
  3. X01=[3;1]; % Condiciones iniciales 1
  4. X02=[0.5;0]; % Condiciones iniciales 1
  5.  
  6. edo = @(t1,X1) F(t1,X1,mu1);
  7. [t1,X1]=ode45(edo,I,X01);
  8. x1=X1(:,1);
  9. y1=X1(:,2);
  10.  
  11. hold on
  12. plot(x1,y1)
  13. title('sol con ode45');
  14. hold off