Hausbau;

% Zu zeichnenden Bereich fuer d1 und d2 festlegen
d1_vector = linspace(2,20,25);
d2_vector = d1_vector;

% Funktionswerte z = K(d1,d2) berechnen
for i1 = 1:length(d1_vector)
  d1 = d1_vector(i1);
  for i2 = 1:length(d2_vector)
    d2 = d2_vector(i2);
    z(i2,i1) = f([d1;d2]);
  end
end

% Grafik erstellen
figure(1); clf
surf(d1_vector,d2_vector,z,'FaceAlpha',0.3);
xlabel('Dicke d_1');
ylabel('Dicke d_2');
zlabel('Zielfunktionswert K(d_1,d_2)');
title('Zielfunktion K(d_1,d_2) fuer das Hausbau-Beispiel');

% Loesung zusaetzlich einzeichnen
hold on
plot3(xstar(1),xstar(2),fstar,'ro','MarkerFaceColor','r','MarkerSize',10);


