function MSE  = expp(n,numruns,numsteps,binsize,featuregenerator,alpha,sigma)    x = diag(zeros(n));    MSE = diag(zeros(numsteps/binsize));    for run = 1:n        w = diag(zeros(n));        for t = 1:numsteps            for i = 1:n                x(i) = feval(featuregenerator,i);            end            s = w'*x;            y = feval(sigma,s);            z = x(1);            w = w + alpha * (z-y)*x;            error = (z-y)^2;            MSE(ceil(t/binsize)) = MSE(ceil(t/binsize)) + error;        end    end    MSE = MSE./(binsize*numruns);