Recently (27.3.) we had a similar question in the thread "plotting results of rk". Try: plot2d( [discrete, a] );