home *** CD-ROM | disk | FTP | other *** search
/ Celestin Apprentice 4 / Apprentice-Release4.iso / Languages / RLaB 1.18c / examples / ode_examples.r < prev    next >
Encoding:
Text File  |  1994-10-02  |  701 b   |  38 lines  |  [TEXT/ttxt]
  1. //
  2. //  Integrate the Van der Pol equation
  3. //
  4.  
  5. vdpol = function ( t , x ) 
  6. {
  7.   local(xdot);
  8.  
  9.   xdot[1;1] = x[1] * (1 - x[2]^2) - x[2];
  10.   xdot[2;1] = x[1];
  11.   return xdot;
  12. };
  13.  
  14. vdpol_out = function ( t , x )
  15. {
  16.   return [t; x; 2*x];
  17. };
  18.  
  19. // system( "/usr/ucb/ps -aux | grep rlab" );
  20. t0 = 0;
  21. tf = 20
  22. x0 = [0; 0.25];
  23. tic();
  24. // out = ode23( vdpol, t0, tf, x0, vdpol_out);
  25. // rfile ode78
  26. // out = ode78( vdpol, t0, tf, x0, vdpol_out );
  27. // rfile ode4
  28. // out = ode4( vdpol, t0, tf, x0, vdpol_out );
  29. out = ode (vdpol, t0, tf, x0, , , , vdpol_out);
  30. printf(" Elapsed integration time: %i\n", toc() );
  31.  
  32. plgrid ();
  33. ptitle ( "RLaB ODE Example (vdpol)" );
  34. xlabel ( "Time (seconds)" );
  35. ylabel ( "States" );
  36.  
  37. plot( out );
  38.