wolffd@0: function [startf,stopf,maxsofar]=maxsum(in)
wolffd@0: 
wolffd@0: maxsofar=0;
wolffd@0: maxendinghere=0;
wolffd@0: startf=1;
wolffd@0: stopf=1;
wolffd@0: aktstartf=1;
wolffd@0: 
wolffd@0: for i=1:length(in)
wolffd@0:     if (maxendinghere+in(i) < 0)
wolffd@0:         maxendinghere=0;
wolffd@0:         aktstartf=i+1;
wolffd@0:     else
wolffd@0:         maxendinghere=maxendinghere+in(i);
wolffd@0:     end
wolffd@0:     
wolffd@0:     if(maxendinghere > maxsofar)
wolffd@0:         maxsofar=maxendinghere;
wolffd@0:         stopf=i;
wolffd@0:         startf=aktstartf;
wolffd@0:     end
wolffd@0: end