% Labuppgift 3 c, minimera plåtåtgång

x1=0.01:.1:3;
x2=0.01:.1:3;
%x=[x1; x2];
[X1,X2] = meshgrid(x1,x2);
V=X1.^2.*pi.*X2-1;
A=2*pi.*(X1.^2+X1.*X2);

% nivåkurvor för A samt nivåkurva för V=1
figure(1)
contour(X1,X2,V,[0,0],'r')
hold on
lv=linspace(0.1,40,20);
contour(X1,X2,A,lv)

% funktionsvärdena för A över bivillkoret
figure(2)
a=2*pi.*x1.^2+2./x1;
plot(x1,a)

% startapproximation
x0 = [.6; 1.1]; 

options = optimset('LargeScale','off');
[x,fval] = fmincon(@(x)fun_3_c_funk(x),x0,[],[],[],[],[],[],@(x)fun_3_c_bivillk(x),options)