As far as the first plot goes, creating a loop works:
figure(1)
plot([1:size(spk{1},2)], spk{1}, dep(1),spk{1}(dep(1)), 'r^')
hold on
for k1 = 2:size(spk,2)
plot([1:size(spk{k1},2)], spk{k1}, dep(k1),spk{k1}(dep(k1)), 'r^')
end
hold off
grid
Other than that, the code uses the data in your .mat files, and appears to me to be adaptable to their sizes.
