% hw2_SOLN
% Method (1)  Use of script and external function

%  where  x(1) = w
%         x(2) = H1
%         x(3) = H2
%         x(4) = a
%         x(5) = b
clc
clear
global L1 L2 h
L1 = 12;    
L2 = 10;
h = 3;
% starting values
xg(1)=min([L1,L2])/2;
xg(2)=L1/2;
xg(3)=L2/2;
xg(4)=xg(1)/2;
xg(5)=xg(1)/2;
soln = fsolve('placetd_hw2_lab2solver',xg);
disp('w = ')
disp(soln(1))
% now preparing plot
% these are just used to establish reasonable starting values
steps = 11;
ladder = linspace(10,15,steps);
for k = 1 : steps
    L1 = ladder(k);
    xg(1)=min([L1,L2])/2;
    xg(2)=L1/2;
    xg(3)=L2/2;
    xg(4)=xg(1)/2;
    xg(5)=xg(1)/2;
    soln = fsolve('placetd_hw2_lab2solver',xg);
    disp('w = ')
    disp(soln(1))
    % save values for plotting
    xplot(k)=L1;
    yplot(k)=soln(1);
end
plot(xplot,yplot)
xlabel('Length L1 (ft)')
ylabel('Alley width w (ft)')
title('placetd\_hw2\_SOLN')