Mansour Torabi Matlab Euler-Lagrange Library
Mansour Torabi Matlab Euler-Lagrange Library
'()*+,-./
012314567389:;1<91=<85;>4919?5@A95B89>C>=<56D712>A=E5A@5B?=1F>G5H?E28FE I743>G
J,KL-MN(O-()OLÿNQR-(-S
ÿTÿ*.(-*ÿ ÿUÿV+-W*
ÿH219 ÿA2>@>G12>A=E
XAY8 Z**,L* [ \,KKÿ-L],L*.* ^_.Q+)* \-+L̀_.* aQWQ bL_,-Q.S
ÿF1E289ÿ Aÿ2Aÿ@>38
8fgfhi j
ÿk!#
012314lÿ67389:;1<91=<8ÿ;>4919?ÿ@A9ÿB89C>=<
6D712>A=EÿA@ÿB?=1F>GÿH?E28FE
n*Q)Oÿ.oLÿ(R+pLÿKQR-(-Sqÿ+)Lÿ_()ÿeL-QpLÿeQVVL-L).Q(KÿL],(.Q+)*ÿV+-ÿ()SÿeS)(rQ_ÿ*S*.Lr*ÿ()e
*+KpLÿ-L*s+)*Lÿ+Vÿ.oLÿ*S*.LrÿV+-ÿ(ÿOQpL)ÿ_+)eQ.Q+)*t
u,)_.Q)(KQ.Sÿ+Vÿ.oLÿKQR-(-Sÿo(*ÿRLL)ÿZKK,*.-(.LeÿRSÿ.oLÿV+KK+dQ)OÿLv(rsKL*w
xtÿy+,RKLÿ\L)e,K,r
ztÿbs-Q)Oÿ\L)e,K,r
{tÿ\L)e,K,rÿdQ.oÿbs-Q)OMK+(eLeÿ*,ss+-.
|tÿy+,RKLÿ\L)e,K,rÿdQ.oÿV-LLÿ*,ss+-.
6}1F~38ÿ[lÿBA7438ÿI8=Y737F
,*.ÿ-,)ÿ.oLÿ
ÿ.+ÿeL-QpLÿL],(.Q+)*ÿ()eÿ*+KpLÿQ).Q(Kÿ_+)eQ.Q+)ÿs-+RKLrw
#$522
#
%! 28
28
!" 02&
01212103ÿ0560ÿ78 9
ÿÿ8
28
!" 5ÿ
ÿ
'()*+,
]^P]ÿ`abÿ`aNÿÿc`abÿc`aNÿÿ
]^P]ÿdbÿdNÿPbÿPNÿebÿeNÿÿfÿ`ÿÿ
ggÿhijk`ilÿmjnÿop`kj`imdÿJjkqf^ÿ
rbÿsÿbtNuebuc`abvNÿwbtNuPbuxdbtNuc`abyvNzÿ
KlN{|ÿsÿdbuc`abulp]x`abyÿwÿdNtNuxc`abÿwÿc`aNyulp]x`abw`aNyzÿ
KlN{^ÿsÿdbuc`abu]ijx`abyÿwÿdNtNuxc`abÿwÿc`aNyu]ijx`abw`aNyzÿ
KlNÿsÿ]}q`xKlN{|vNÿwÿKlN{^vNyzÿÿ
rNÿsÿbtNueNuxc`abwc`aNyvNwbtNuPNuKlNvNzÿ
rÿsÿrbÿwÿrNzÿ
KbÿsÿPbufudbuxb~lp]x`abyytNzÿ
KNÿsÿPNufuxdbuxb~lp]x`abyywdNuxb~lp]x`abw`aNyytNyzÿ
KÿsÿKbÿwÿKNzÿ
Mÿsÿrÿ~ÿKzÿ
ggÿ
}ÿÿsÿ`abÿ`aNzÿ
c}ÿsÿc`abÿc`aNzÿ
``ÿsÿdij]mlkxyzÿ
J}ÿsÿMmfqmjfkc^jmPilJ}ckqi
kqxMÿ}ÿc}yzÿ
ÿ||ÿsÿc^jmPilJ}pd
kqxJ}ÿ}ÿc}ÿdbÿdNÿPbÿPNÿebÿeNÿfOOOÿ
ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿOÿOÿbÿNÿONÿOÿObÿ``ÿbN~tbuiyz
-./*+(ÿ12ÿ3145/+ÿ6+.34/47, -.87)9+3ÿ:+(61.(+,
;<)76/+ÿ=,ÿ>6?8.*ÿ@+.34/47
ABCDÿFBGÿDHIÿJKLMNOPÿDQÿRIFSTIÿIUBVDSQGCÿVGRÿCQWTIÿSGDSVWÿXQGRSDSQGÿYFQZWI[\
#$522
#
%! 28
28
!" 12&
01212103ÿ0561ÿ78 9
ÿÿ8
28
!" 5ÿ
ÿ
'()*+,
_`R_ÿbcÿdbcÿeÿdeÿ
_`R_ÿRÿfÿgÿhÿbÿÿ
iiÿjklmbknÿolpÿqrbmlbkofÿLlmsh`ÿ
tÿuÿvwxyRyzde{xÿ|ÿzfÿ|ÿe}{xydbc{x}~ÿ
MÿuÿRyhyzf|e}ynr_zbc}ÿ|ÿvwxygye{x~ÿ
OÿuÿtÿÿM~ÿ
iiÿdmskmÿLobkrl_ÿ
ÿuÿbcÿe
~ÿdÿuÿdbcÿde
~ÿ
LÿuÿOohsolhmd`loRknLdmskmszOÿÿd}~ÿ
iiÿrfmÿLobkrl_ÿ
bbÿuÿfkl_onmzvP}~ÿ
ÿee
ÿuÿd`loRknLrfmszLÿÿdÿRÿfÿgÿh
QQQÿ
ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿvÿvÿvÿQv
ÿbbÿwvykQvÿÿ
}~
-.*/+ÿ).1ÿ/+.*23ÿ45ÿ(678.*ÿ6+.19/9:, -.8:)2+1ÿ;+(64.(+,
<=):6/+ÿ>,ÿ?+.19/9:ÿ@823ÿA678.*B/4)1+1ÿ(966472
CDEFÿHDIÿFJKÿLMNOPQRÿFSÿTKHUVKÿKWDXFUSIEÿXITÿESYVKÿUIFUXYÿZSITUFUSIÿ[HS\YK]^
'()*+,
_`R_ÿbcÿdbcÿeÿdeÿ
_`R_ÿÿRÿfÿgÿhÿÿ
iiÿjklmbknÿolpÿqrbmlbkofÿLlmsh`ÿ
#$522
#
%! 28
28
!" 62&
01212103ÿ0561ÿ78 9
ÿÿ8
28
!" 5ÿ
ÿ
N`aÿcÿde`ÿfÿgheijhklmdijnnoaÿfÿdgheijhmpqdijnnoarÿ
sÿÿÿcÿtuahShN`aÿfÿtuahvhe`oarÿ
NÿcÿShwhghdtxklmdijnnÿfÿtuahyh`oarÿ
PÿcÿsÿxÿNrÿ
zzÿe{|p}{ÿM~iplqmÿ
~ÿcÿijÿ`rÿe~ÿcÿeijÿe`rÿ
M~ÿcÿPw|qw{eqSpkM~e{|p}{|dPÿ~ÿe~nrÿ
zzÿ
lg}{ÿM~iplqmÿ
iiÿcÿgpqmk{dtanrÿ
ÿ``ÿcÿeqSpkM~
lg}{|dM~ÿ~ÿe~ÿvÿSÿgÿyÿwRRRÿ
ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿaÿtÿRÿÿRtÿiiÿQuthpÿÿnr
()*+,-ÿ/01*2*03ÿ43+ÿ/,3+5)56ÿ73)8,9 73*642,+ÿ:,1;031,9
<=:<47?<
6;) ,ÿ>9ÿ?05@),ÿ/,3+5)56ÿA*2BÿC-,,ÿ15;;0-2
6+
DEFGÿIEJÿGKLÿMNOPQRSÿGTÿULIVWLÿLXEYGVTJFÿYJUÿFTZWLÿVJGVYZÿ[TJUVGVTJÿ\IT]ZL^_
mSmÿ`ÿijtÿijaÿe`ÿeijtÿeijaÿÿ
mSmÿvÿStÿSaÿgtÿgaÿwÿ
zzÿpq{ipkÿqÿli{qipgÿMq{|wÿ
}t`ÿcÿgtheijthklmdijtnÿfÿe`rÿ
}tÿcÿgtheijthmpqdijtnrÿ
}a`ÿcÿgtheijthklmdijtnÿfÿgaheijahklmdijanÿfÿe`rÿ
}aÿcÿgtheijthmpqdijtnÿfÿgaheijahmpqdijanrÿ
}tiÿcÿ}t`oaÿfÿ}toarÿÿ
}aiÿcÿ}a`oaÿfÿ}aoarÿÿ
#$522
#
%! 28
28
!" &2'
01212103ÿ0561ÿ78 9
ÿÿ8
28
!" 5ÿ
ÿ
_ÿaÿbcdefeghidÿjÿbcdekbelbmÿjÿbcdekdeldmnÿ
obÿaÿkbepeqberbstuvrmwbxxnÿ
odÿaÿkdeperqberbstuvrmwbxxjqderbstuvrmwdxxxnÿ
oÿaÿobÿjÿodnÿ
yÿaÿ_ÿsÿonÿ
zzÿ
{ÿÿaÿ|h}ÿmwb}ÿmwd~nÿ
g{ÿaÿ|gh}ÿgmwb}ÿgmwd~nÿ
mmÿaÿqvtr}d
}
xnÿ
{ÿaÿyppgkt{glry}ÿ{}ÿg{xnÿ
|}ÿhh~ÿaÿgkt{uqlr{}ÿ{}ÿg{}ÿ|fÿkbÿkdÿqbÿqdÿp~}ÿ
ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ|
}ÿ
}ÿd}ÿb}ÿb}ÿb~}ÿmm}ÿ|}ÿc}ÿdec}ÿ}ÿ}ÿ~
()*+,-ÿ/01*2*034ÿ/,3+5)56ÿ73)8,19
:;<=>?ÿABC<D<BEFÿA>E=G;GHÿIE;J>CK IE<HLD>=ÿM>CNBEC>K
OBEDLPD
Q6R*)9ÿ1620-RRS*TU6R*)VW06
M>;>LC>CÿX
Y<?CDÿM>;>LC> ZLD>CD
03ÿ[R3ÿ\]
ALP^LJ>C
#$522
#
%! 28
28
!" &2'
01212103ÿ0561ÿ78 9
ÿÿ8
28
!" 5ÿ
ÿ
'(ÿ*+,-+./0ÿ*123405/6ÿÿ
789:;8:<=
>?@7?ABCCDCE
#$522
#
%! 28
28
!" &2&