>> r=0.05:0.1:3;
S=r.\1+2*pi*r.^2;
plot(r,S)
xlabel('r')
ylabel('S')

[x,f]=fminbnd('1/x+2*pi*x^2',0,1)

x =

    0.4301


f =

    3.4873

>>