function [s] = sse(y,yp)
ye = y - yp;
s = sum(ye .* ye);