Hello, It seems to me that rk incorrectly increments the time variable... (%i1) fpprintprec:5;rk(x+t,x,1,[t,1,3,0.5]); (%o1) 5 (%o2) [[1.0, 1.0], [0.5, 2.4453], [1.0, 4.5036], [1.5, 8.2207], [2.0, 14.672]] (%i3) I use Maxima 5.30.0 (rk is coded in Lisp) - GCL 2.6.7 (Debian testing). With Maxima 5.29 (rk is a macro), there is no problem... JM Sarlat