load(numericalio)$ load(lbfgs)$

data:read_matrix("~/Dropbox/preprintsNEW/MNU15Augsburg/Ferienh_rohdaten.csv")$

vars:first(data); data:rest(data);

prognose(datensatz):=
block([Personen:datensatz[1],Wfl:datensatz[2],Meerentfern:datensatz[3],Pool:datensatz[4],Hund:datensatz[5]],
  ps*Personen+w*sqrt(Wfl)+m*Meerentfern+p*Pool+h*Hund);

abstandsumme(data):=sum((last(data[i])-prognose(data[i]))^2, i,1,length(data));

optimum:lbfgs(abstandsumme(data),  [ps,w,m,p,h], [1,1,1,1,1],   1e-2, [-1,0] );

subst(optimum,prognose([5,100,200,1,1]));