%file name = tre_freq_opt.m
%redundant truss optimization for frequency constraints
%main program
%Prof. Reyolando Brasil - May 2021
clc
clear
%problem data
Prob_data(1)=1000;%density, kg/m³
Prob_data(2)=100000;%Young's Modulus, N/m²
%options
options=optimset ('LargeScale','off','TolCon',1e-8,'TolX',1e-8);
%Lower and upper bounds of design variables
Lb=[0.01 0.01];Ub=[1 1];
%initial design
x0=[0.01 0.01];
%optimization function call
[x,FunVal,ExitFlag,Output]=...
    fmincon('tre_freq_obj',x0,[],[],[],[],Lb,Ub,'tre_freq_con',options,Prob_data)
%