clc; clear all; close all; t=0:0.1:20; opt = optimset('Algorithm','interior-point') x = fmincon(@optPendControl,[-0.2 -0.5],[],[],[],[],[],[],@const1,opt) %x = ga(@optPendControl,2) %x = particleswarm(@optPendControl,4)