% % dbss.m % function ss = dbss(params,data) global Re Pr Nu = data.ydata; q1 = params(1); q2 = params(2); q3 = params(3); res = q1*(Re.^(q2)).*(Pr.^(q3)) - Nu; ss = res'*res;