
% Setzen der Problemparameter
V0 = 330;  % gefordertes Volumen der Dose [cm^3] 

% Setze die Zielfunktion
f = @(x) Dose_Zielfunktion(x,V0);

% Setzen von Optionen
options = optimset('LargeScale','off','Display','iter-detailed','GradObj','on','DerivativeCheck','on');

% Setzen einer Startschätzung
% Bedeutung der Variablen
% x(1) = r = Radius [cm]
x0 = 2.5;

% Aufruf des Optimierungsverfahrens
[xstar,fstar,exitflag,output,grad,hessian] = fminunc(f,x0,options);
hstar = V0 / (pi * xstar^2);


