function f = objfun1(x)

beta = []; model = []; w = []; goal = [];
load('objfunw.mat')

f = 0;
for i = 1:length(w)
    if w(i) ~= 0
        load(['objfun' num2str(i) '.mat'])
        D = x2fx(x,model);
        f = f + w(i)*(beta'*D'-goal(i))^2;
    end
end

end