function J = sumsq_allom(p,ldata,wdata)
% Function computing sum of square errors for von Bertalanffy model
model = p(1)*ldata.^p(2);
error = model - wdata;
J = error*error';
end

