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

