Raster and PETH Simple Version

Raster and PETH for neuroscientists and others :)

event1=rand(100,1)*1000;%for neuroscientists, this is the timestamp of the node;
event1=sort(event1);
event2=[event1+10+rand(100,1)*2;event1+10+rand(100,1)*4; event1+10+rand(100,1)*4; event1+100+rand(100,1)*10]; % for neuroscientists, this is your spikes :)

[rst,nodewndw]= RSTPTHSMPL(event1,event2,50,50);
[h1,h2]=PLTRSTPTHSMPL(rst,'Event2',2,50,50,'around Event1');

% More Info, refer to .m file