%file name = cable_opt.m
%cable nonlinear problem
%main program
%Total Potential Energy
%Prof. Reyolando Brasil – May 2021
clear
clc
%Design variables:  x(1) = u   x(2) = v
%
Prob_data(1)=200;%vertical load V, KN
Prob_data(2)=100;%horizontal load H, KN
Prob_data(3)=1;%cable cross section, cm²
Prob_data(4)=21000;%Young's Modulus, KN/cm²
Prob_data(5)=100;%La length, cm
Prob_data(6)=100;%Lb length, cm
Prob_data(7)=50;%N0 KN
options=optimset ('LargeScale','off','TolCon',1e-8,'TolX',1e-8);
x0=[10,30];
%optimization function call
[x,f,ExitFlag,Output]=fminunc('cable_obj',x0,options,Prob_data);
disp('     u       v')
disp(x)
disp('Energia Potencial Total Mínima em KJ')
disp(f)
%