% zmienną jest r
% H = 0.33/pi*r^2 , 0.33 to objętosc
% materiał na ścianki i denka kosztuje tyle samo

x = 0.01:0.01:0.5;
k = (sqrt(5)-1)/2; % wyklad s 5
a = 0.01;
b = 0.5;
xL = b-k*(b-a);
xR = a+k*(b-a);
while (b-a)>1e-6
    if ((2*pi*xL^2 + 0.66/xL) < (2*pi*xR^2 + 0.66/xR))
        b = xR;
        xR = xL;
        xL = b-k*(b-a);
    else
        a = xL;
        xL = xR;
        xR = a+k*(b-a);
    end
end
r1 = (a+b)/2
fun = (@(x) (2*pi*x^2 + 0.66/x));

figure(1);
hold on
plot(x, arrayfun(fun, x),'b'); %rysujemy wykres funkcji f dla arg x
plot(r1, arrayfun(fun, r1) , 'r*')

% materiał na denka jest dwukrotnie droższy niż na ścianki
k = (sqrt(5)-1)/2;
a = 0.01;
b = 0.5;
xL = b-k*(b-a);
xR = a+k*(b-a);
while (b-a)>1e-6
    if ((4*pi*xL^2 + 0.66/xL) < (4*pi*xR^2 + 0.66/xR))
        b = xR;
        xR = xL;
        xL = b-k*(b-a);
    else
        a = xL;
        xL = xR;
        xR = a+k*(b-a);
    end
end
r2 = (a+b)/2
fun2 = (@(x) (2*pi*x^2 + 0.66/x));

figure(2);
hold on
plot(x, arrayfun(fun2, x),'b');
plot(r2, arrayfun(fun2, r2) , 'r*')