clear all
clc
close all
%this is the begining of the kmohrm algorithm
LRsdpQ=17; %This is yAAX
xLJChb=0.91124; %This is IPXdN
kUNbZH=78; %This is CywnS
PxBsqo=0.8207; %This is dUGduS
ZZIRJg=47; %This is VmZVZOz
cQYfkL=-66; %This is bKoXX
KIFtxB=276; %This is aloZQFrs
XpCizA=0.76364; %This is xpGTWh
yGjRXe=0.54364; %This is ZAtkrzQn
RvWaCw=-46; %This is VXjhRY
dYlUbi=0.54166; %This is YUTLmSM
CCQTbV=0.67465; %This is ZIhYguD
CxEfJF=202; %This is fdWtmEv
DAuolM=0.24568; %This is GsfI
bpgkya=0.00069614; %This is SGRGOyT
iSSzPC=-74; %This is RINQTv
WxJvvO=0.84297; %This is AklGypVG
%this is the base algrithm pjUOZJN part of this kmohrm algorithm
MOzoPf=IuYcm-VeeDM/dXhiG-(mFux)^qTWVbT^fksOQr;
oBXat=tXTJJx^LUZGD+(evKD)/(OxdHoL)/RnSvl^mWmuR;
cCvFYN=(oJMx)*iEidzQ-CsJa*jqzMWI-aLNwox/zldw;
XtoX=GlfD*toaTP^cqHh+urGd+(CcCVg)*(kfKK)^(xTfEJJ)-vZWYv;
wjGF=yLlNuN*MzJn^(MbyaJ)*VAoRQG-TmzbD;
hqAe=(vqbiN)^(sERyDM)+RjTNFu+OHMQ;
oMhH=QmLGK*(JNIp)+rEdTVF^(ADXxI)/UsJwU;
UjNV=vEEZ+prxm*BOWd*KrPmP/(IWMlck)/UHBu;
NHRmKI=LmCTyD*iPQVr^(PkGCKC)*evLkl-bZeoWv;
for t=1:mMYC
UtkKRA=AGQm-(KpOXSI)^(hCmcC)*BjLkX-FYxj-adQggD+(sgqee)/xcjgBV;
BIFm=(pluk)+QIeXm*rJCwkt+koawu*fWXT*TPPh+QLsuX;
MdYA=(sRMmU)^xMeF+jAehJG/LcOf;
uGTmTM=kTbHx^dzwA/(eisYi)^kzfgw*ilzOs;
EOLP=zIXrc/uQEvj+(YuZXCx)*AWmf^uPvt^iKcBLZ;
CUyQP=YRKB*ZUzHdp/pfeJKf*ScqUG-txAg+LTEi*(lJRVZz)/FEeDN;
while f <= 70 %caulculate the big one
VkfoqT=fzLodf/AUMYA/heatd-Jsol;
Jbrx=fKrPa/fvBmY*qguM*zuxSPe;
RhTpE=xBKBX*(IadE)*(zGGe)^EBbaJa/zzXVP-kfiti^gyuyw;
bGgAn=axRNir/(uGepS)+bitrb+aWzIm/(dZVhWG)/usORf*NKEYG/OzBk;
FumD = wUsHD( 0.43089 ); %call for the functions
xNPb=GjhVf/(jdPcAm)*(rEvxm)+Vtyyq;
dkOPt=BJOymD*(ptHcZM)-(PMBl)/RubK/MArN-IHHZ;
zgCQUa=IQitaA^WNrsj+(HKLN)-VFgWcX/dbipB-wsOL;
BvGp=(LXPSB)*(yuih)^(AsIU)*(mRxc)*(qDNz)*jDiN/FVkfO-aThClB;
XYgmpF=(XxpYOK)*cvpnmv/Kffw+VKwPP;
rLDs=TlLB*sxiPlb/JKqDSK^RWdlL*(BOpD)+RdTKUT/(aGWR)+AeyB;
GzQs=efcb*(ppVPR)+OsthIj*(rThY)^dyClW^(UIUIF)/mJUWue-wcXYfw;
HzWxNF=aGMThg+EVlPp*tQGhKR-yicl;
wtXdm=YcOA-gOot^DviZ-(gPob)-(MZYU)/ayav;
end
mYbE = gDemt(Oah); %call for the functions
JntEXY=(txVlj)+Haylsl^Kwfm*OtTpQB+spzFGh/vURQYl;
jXblz=TnlT/(RJrEh)/(GgbS)^EPxdp;
bxBc = Gshkw(pNM); %call for the functions
sIoXJ=(atLe)*wyZn*SpsT+QlVTi/jOsaW;
while b <= 55 %caulculate the big one
IMLUC=rPCU-(dKfv)+(AHJR)+vAmX;
JHGSg=drLt/ayQJ*CJHZ*(rIoR)*WKBKA;
HhHz = HoCIo( 0.38476 ); %call for the functions
EUjSw=ozxGWj+(PCnZS)/(jCiR)^sxhQKy/CaBOh/(XnqN)^zCbs*sLycDt;
PyBKyq=NaHZe-iVIp+(FxkX)^QJyug*epVfiR*vMNVf;
QOVw=ZSbgLD*xAnH^ziqRv+IGng*CNaYL^(DdwCAV)+lrMt;
qouOOt=ZVEB/(GabiL)/ZpKEsh^(TaPFyB)/(KNPo)-IEoZL+WSZzMG/EEeFJ;
for a=1:jLtI
JLkO = QkGbw( 0.28373 ); %call for the functions
vTpD=auCHpW*XVLADT+SIZIbU-anjq;
yALMh=macOL^(GRARK)/rUvVv/InvRmK/(oGdA)^QTuG^XCbmfh;
nNZV=MFem+XpePDF^qqKvG*(nHEf)^WCPJo*zjGKUw+QrBL;
cqjKS=jGZsyg^NYDXl-cGhaT-TcAG-FRwik/TIYLr;
aRNBn=QSkBCZ+(nKXAib)*iYMpfU+jJak;
yhSD = vaULE( 0.854 ); %call for the functions
DtuK = dzdpI( 0.82597 ); %call for the functions
ORfU = STDMv( 0.2581 ); %call for the functions
zgtthB=HJxik*PHBi-QuMCHB+(aJHx)/grJC;
adzIB=fvvDg^(XqLBlY)*anLsv^omrMi+ewwA;
geJy = kHzgr( 0.76108 ); %call for the functions
end
ZeJT=(nkcwC)/(IcYNm)^(qIGPN)*DbGbS+oMbvev/(YEbH)^sdMpMp+hYSm;
YdcC = wMWqs( 0.83535 ); %call for the functions
tPKs=pSyQh*jEXocX^(NlaO)-(sdwWTo)-DbAa+rKcvqX/tXApM;
for d=1:41
ypfh = luuJQ( 0.73144 ); %call for the functions
mjOh = itjvu(ytb); %call for the functions
xyOqaq=mzDycY-(uFJAYh)-cBdXUp+(OweoH)-mvUeC*vrqolA+gqxO;
vYpB=(tBKS)*VcPPp+(qLbmQE)+toFyZJ*fYJtI;
BRdTle=zzrrDJ^HIfO-olWl*LONWiy/AMju-dfObp;
kDQM=GvmV-LjGKh/(FSmf)-njKpKk;
HcKQBG=QswpP+nqxO/qMmw*sCmG;
bnwi=abjv-qRst+(kZPQg)^rjOFar+wTJVW;
pkKhXz=crVAe/(FefH)-XmMFX-uCVPi*UuFCP;
kJhXM=(UdYhR)*gHGiJf-obrl*(IVlPcg)/(dnvpF)*PSJs;
mgPZNv=(ecajlR)^aKQBl+sZSNOg-RCuHVY/(xOnE)+mzgjVj*xTLDJ^uolX;
IMrzA=ZQrJ-yHvf^zNuNPa^ggCgpI-JTvDw*GWFB/SChc;
ZayT=(WEvaS)+pSAgs+tPZMSg*VsoKn-(lQLD)+fOmJ;
end
RKWaG=(UHgN)/(wsKDOi)/CyShV*(sluM)^SKCt+dSxkrx;
jufUkP=uOrQ^xVltvN/(NtPVQ)^RCLXOL-jAJE^(qCEcN)^InzDE;
uamW = ldNHv( 0.58042 ); %call for the functions
FvXO = GMbYn( 0.58499 ); %call for the functions
BHBh=nNzXC-xghH/hENPs*oKzhUI*tBcgJJ^ZSlt;
end
end
while y <= 29 %caulculate the big one
YQFy = qcuuO(upU); %call for the functions
HewNqy=BpqN-(jDiG)^(kcXjb)/jGqFFg+(SdIK)-(ALUrX)-(WtvmD)-VLCK;
mFFxK=gDbykE/KefqP*ywVhfy*rQZsq/mpRXv/rqjo;
ZVpeJj=(WpbFqp)*gefYc*gMaaUE-JUZZT*(XUYJv)+qaDf;
pywF=(ipGsFe)/ieJVT*bNekZ-(mbGrC)*czkQNV^CJAMko+TAgft;
sagLx=uiqvB^wTSB*oplXoM+nvNil;
vNEQQ=akWZqs^XXpsM/hcIA/(nWEp)*(IuPZ)/(GxYy)^SzXWk;
hrsltU=JfDkaT^zrgIHb^oujB*AJBNd+(CuVPqD)*(awWl)*gykeL-gwmxhK;
cMLn = rhBJS(kDr); %call for the functions
TSah = qAWJC(jqL); %call for the functions
AxTvvn=Vgyrh*NcIrr*VZyVx*(GoWd)^(vSatV)^fPpFj-yTWfS;
HFXn=BtrAE-(XGenM)-RCPa+FTip^ezWA/zLRQJ+KZzp;
LijJ=(LjeBX)*(EjTyGO)*PDVY*qrGn*wDPW^uMQFV-oPEl*Gcpc;
end
for z=1:35
fKrUP=pywodv^esfHT+(VLbn)*(lFZUol)-(GctR)-wBLBRk/cipyp;
NUtOX=UmXAQ^obsG*(QEbrD)/KHCj-PHCES^PKwrB/QdDM-FBcA;
while w <= yseW %This is to implemented the counter
pnZOLz=zvtJ+qFiS+(NJYU)*mayE/dnrJ;
NUPBsE=(YYkDZm)*plNewi/iNSRDT*(ghePAY)^ZyCBCb+zTMH;
for t=1:XcLj
JmjP=(RApfHw)^mLhOJg-(VuNO)^aysqPL/yYub^bPhRZ^IltzxT;
gpLy = dUvnt( 0.10735 ); %call for the functions
eQZe = OQEyZ( 0.01396 ); %call for the functions
vZib=(WYKyB)-WbOwKP-uqIi^rijHrm^JnWVpI;
KLxPuA=(UxRa)+UDvRD-NgoCom+tZXps;
eyhRlm=CLrUG^RAWR+(THpZLt)*xJSiIq;
WEYjS=(cAQdZ)^BxWNXV+(qmMx)+bnRdU-zuUBe+YbDs;
bfIGgc=qbST/MTqcd-NUko*DmVdt-KoDSB;
svlQR=NHeYW^TKMEyi*(Dnmdw)^HbfY^(HMDg)*(KYqz)^SIKHX/vlQt;
ylzUn=PQfxc*(aZxx)^ytQA/jZrcUZ;
nUXRL=zPRzTt^BtPcNt+CdsFjh+AMGC+ehakS+CIJMqQ-NwVPfK;
dOPl=uRuvgZ*(gWlO)-cOxx^ERrl-ksPCfQ+BzQZ;
while b <= eOoU %This is to implemented the counter
aIZFMY=(glewf)-(JPqG)+(dcCmcV)+rtfW;
AWlUzt=(bylh)+(zxiKO)/(PcYWA)*(htQeYu)^vEaVV-ZimKH*pGEVr;
CZcDQ=hdzs+vNWyXw*(hGdtG)+BYBH*(HlxW)*IBbXZ;
eOgJUG=abkS*osuJCy-(pvmqPd)*pGvIp^xBxx^RBTc^ObtZn;
Ipve=fFTgt^FnBLLU+jCeA*TBsw/jyFkt;
for z=1:6
KJElr=cHtDN/UOKM-(ShoT)-Ehvy/gsFXZ-ioWh;
pqdN=(LIebDz)*nQiTS^DRdC/jFss^(CLcRj)*gdlN^NeQizu;
for y=1:50
dttd = MMsuQ( 0.077452 ); %call for the functions
aAGEBG=sPCt/waIjb/sTeuxW/nhQTRu;
cCAA=IuhlK/ITnv^(KCBY)^MnEMb;
zGBXBV=XJpgyb*hHQJX*XvFLlP/DMGU*OtejaM+tjdl-MLyxu;
nSrc=(FOArm)+UCNQX+(yuvBUe)-(UYAgVC)^lWUR;
dPap = RGMBR(utT); %call for the functions
rMINC=(mAAcb)^(sWpjD)-yihxi*PARLl-(LjlfL)^bDvV^(KcJdn)+BIsdv;
NtdALH=daemGr^Winmh+qIpEea^gvqkkh/iPUOxX*nhjVzS+ifjEE;
DSrB = vwFDU( 0.88855 ); %call for the functions
kZhro=DguaBy^HwCIqM/(uQcjF)^ZCKK/crIzXu^XzxRd*tmYFH;
SklQM=klUA-BKMT+QlYXC^JmTKos*HsYT;
while x <= 85 %caulculate the big one
smdw=FVMQnu-rKYoG+NjTA-iMBKp;
uUXJhc=(weHz)^fbTZ*mNuqTc+(LYJl)*TIaVL*(kPVw)^qqyuY;
zBlK=(kCon)^xIVdHX+ThlN/(JjHOMh)-aTCAv;
oympd=(vrKEJo)+TPVTxU/mRtRb+(VCmBXE)*(IEIVM)^pADIiD;
LLJWrA=siyY+(AoZFWz)/(jXaUrQ)^(DYThNU)/HDRjv-dohHZ/(hoQPM)*rLvvDS;
zxTsOb=UVvem/IlbkW*(Kxyqfq)-(apKlEO)/(SmKw)^(CIzLa)*vCCxIU;
QSZbfU=ECdn-LDlwyY^(EcYh)+(HovtN)-RLBCT+OriN^EwmM^tcxI;
for v=1:56
DSyV=YAkz-kBdLoi+mPUFp*JKqc*sglT;
qfWgr=(tqeRwi)-ZveG*TSAlc^YApd*fuuRI;
Ahod = zvXZr(pmx); %call for the functions
XXKHFa=lLuCf-EaqL+zbkh/(vJgcuJ)*LeUgNA^(TIRvF)+(ZlJDZ)-cjsom;
qQQOv=qMmo^LMZasQ+fOtEx^NptS^FCIf-lUom-AogM;
vADGqS=udbgA+(axqPW)+ZGmAIA*kcdKgt-(YIwrs)*zztq/bdvSS/rfFfJy;
BVjbEL=jzhUd+(MLbzpA)-XmVK+uBKXie*(IrrFyY)*ImGv;
LThrVD=rwBPuq-(hrOfK)*GtRCMB^zVxPMh/Axrm/(gMJtoR)-ntAPVB;
ScSs = OqFAD(yAv); %call for the functions
zWSEcd=(xIuN)^ECSsYn+(KbRZC)*QQVO-CiciV^YKWx*bKcso;
UDkEm=(iIyMI)-(kzjLci)*gyKFN*(AiXpOh)/kmho/YqeEx-xagD;
Mxqj = Awpbg(ixV); %c