Watch this Author's files
Really great! Thanks a lot!
just one small problem with log scaled axes. To reproduce:
figure; plot(0:5, 0:5);
set(gca, 'xLim', [0 5]);
set(gca, 'xScale', 'log');
Solution: always set lower xLim > 0, if you use log scaled axes.
(I know: quite obvious! :-) ...but this case could simply come true)
I guess this exception could be trapped in function get_projection (around line 1913): if xi(1) == Inf ...