Marduk Bolaños

Implement the multioutative property for operators

Implement the multioutative property for operators

plot2d does not overwrite xlabel

plot2d with a single function now uses xlabel

Implement the necessary limits to evaluate this definite integral

Now I can't reproduce it either. It is very strange.

By the way, the default accuracy, defined in the beginning of rungekuttaeval u := accuracycontrol(u,20,6); is larger than machine precision and very likely unnecessarily large. My experiments with the harmonic oscillator show that with low accuracy (e.g. 6) after 500 oscillations and using 14,000 steps the energy is conserved. Which might also suggest that this method is sympletic. If someone knows the algorithm behind the code it would be good to explicitly mention it in the documentation.

It turns out that rungekuttares was not written carefully. Here is the corrected version, that checks the size of the list before applying car and cdr. symbolic procedure rungekuttares(l,st); % eliminate intermediate points. if st=iterations!* then l else << for i:=1:iterations!* collect <<for j:=1:m do if (length l > 1) then l:= cdr l; car l>> >> where m=st/iterations!*;