% ECE 321 Lecture 8
% Hanging Chain of length 10
% Calling routine
% y = -[1:5]'/2;
% cost_chain(y)
%[y,e] = fminsearch('cost_chain', y)

function [ J ] = cost_chain( Z )
 Y = zeros(11,1);
 for i=1:5
     Y(i+1) = Z(i);
     Y(11-i) = Z(i);
 end
 PE = sum(Y);
 L = 0;
for i=2:length(Y)
 L = L + sqrt(1 + (Y(i) - Y(i-1))^2);
end
 E = (12 - L);
 J = PE + 100*E^2;
 plot([0:length(Y)-1],Y,'.-');
 ylim([-5,1]);
 clc
 disp([J,L]);
 pause(0.01);
end