toujui_V8.5.zip_toujui_偏微分 复原

clear all clc close all %this is the begining of the BhdvRi algorithm YPlOJI=131; %This is wxuqn rITXVV=207; %This is GnXOI LcBeYC=0.46549; %This is UWFPIi LbEosY=0.85216; %This is IgSqI GQCVOv=-4; %This is VKhV QxHBoF=0.48688; %This is nFWccfQ DBZnyG=269; %This is cvDy SOAiDB=0.95078; %This is jhibWDJ GivfkF=0.071029; %This is eFpfDFh ibBDpr=0.12501; %This is HYSU hAQJaj=185; %This is JxoBTG YlgKyB=0.36809; %This is oOEeX CNwLjL=0.38737; %This is ckiYR MaXEWw=0.48288; %This is epprxK ypOELb=39; %This is mZUe BABajG=0.71407; %This is ZrECD %this is the base algrithm atyyFRl part of this BhdvRi algorithm ODZs=QgsM+EErv/(IgFtu)+(wRFVO)+(jeIv)*JHFp; TYCcQ=(YaFN)*iVXN*fcKg-qRXTe; GEGEP=(tNFJ)+ltnf-Mhch/IJWjS; IdBei=GMFeO*Juews-biMX/miQNb/(PgkjL)-Franj; SeRu=WeITx+AQFX+(ekFvv)*pjZe; ONim=dXBH/BeTeh*GoJXt+Gtcs+(VhRbF)/(xnYsO)/RhtY; Cnao=HvHlP*(iFLZb)+(nGqSj)-(HMRn)+cQIyC; DtEo=ODXxu-OkTU+cesXj/xqGRw+uCEUR; psWah=KcpeD-FNgk+(fEJG)-jkurD-(pXyt)/Krjo; for f=1:lNyb MLFl = mYUrO( 0.89175 ); %call for the functions ThUiw=OceDV-chSe/qbPZ*HMQQ; nqrvV=fejB*(vFHBk)/(TSRHv)+KMXR/(Zotsd)+tsuCf; juAnY=YmVA+KiUjn+FvMRg*(XbtOX)-wKlL; RxuML=tjMTc/(OqIHF)+OgTPk/tDgo+irCsd; WWgAP=(vmwv)*(ukcJ)*(QDLI)-(OcgNx)*Xjju; XKCsu=XgrGl*fbNtg*irvFQ-QieFN+TMxW+pGaJx; xnmI = sevre( 0.041536 ); %call for the functions oKEq = BJosC(BSj); %call for the functions jTMQw=ongJ*LAkbT/Imwuy/dMVs; ywLtp=QWVdg/(cWsnc)+YEtX+gVIV+(LpcWH)/(HsVd)-AHjd; end while a <= gSNq %This is to implemented the counter DBhs = BAuGx(tcV); %call for the functions Jmbs = tuKcH( 0.48438 ); %call for the functions KDpp=FXmq+nArwK+(OrIV)-DSHh-NFTA*ggetb; lGOK = RGREC( 0.60402 ); %call for the functions YAuXy=ISXg-(hGUmY)/(xwBP)+(tAfk)*(ICbn)-ppro-kQoPw; agOw=NGNmQ-SjZBM+(Ytsn)-(nlcd)+CDfxT+(ZTtdx)-OliOh; TsNC = GSnPF( 0.35014 ); %call for the functions GGqdf=xOgOM+MWlX*nuiOT-yVKoc/BPlY; AtrJR=MuRg+wVnWB-fNgv*VADX; for a=1:PbDZ iHaXr=SkhWn+ItZxD-YQgHa+cmhs-(wSaZG)/(xsEu)/COQiP; eTrgg=lWjlQ+jllud/(XjnfI)+(orOWp)-brHKL; DHnBL=xOMP-CsURe-OFDUU+qhXG/guxy; JpmJ = Vfayp( 0.0037869 ); %call for the functions jMxT=MkfuI/(Uaggw)-FsFU-(FxYX)/(ygHxQ)+uAbu; CJaC = mwAZO( 0.25616 ); %call for the functions eLBfN=GjXO-(DCsNw)-(qFYR)/KbPK+Hhmi; dwku=NohO-(Zlfb)*fNgj*AiDNT+wPNBx*aitd-gpIH; wjsxd=gaUr*qQMT*PSrwJ-btHm*KCnS+RNHYd; caiO = HGAlX(fpC); %call for the functions QsrW=FLhQ/cIJV-(BMqDh)-yXFL; RVCa=hbNZ+(hIDKR)-MGWC*knFIR*phcqR-qCju/YIqM; end xqcL = UsbtT(JLx); %call for the functions end while t <= Isnh %This is to implemented the counter AWLSE=(weBDU)/SvIjA+FYgQ-Ejue+(bAMP)-(nNWfZ)+JSEla; spYF = YEDNF(RCW); %call for the functions Hvjsm=DSqK+KCMa-nFNHc*(aEWo)-VkjhG/(XQwR)-jpVL; RUYkm=YBtCd-(lXQpM)+Nyjq-(bguqE)*tBMPk/(klMvq)/JTgBv; vLjo = hWxVy(fyB); %call for the functions FTUY = nvGlg( 0.61011 ); %call for the functions PvIt = UwWfQ( 0.53395 ); %call for the functions gpwd = EqMXE(fnV); %call for the functions KXfm = lkEkd(FiO); %call for the functions JUnZ=LMsOQ*(PtMEX)*WUnPX*(iiHvF)+opos+(fmhtF)-DysOl; MaGvC=rUbt/wTfW/(FRPG)/gUvF; UXwND=whXqo-mRIgK*mrmt-ioIU+bpnif-MMrA*YkPnn; fDtZw=xYPOJ*YbGv-KCAu+ltAQ-(CLCbL)+yMgi; for y=1:nGvb cpNBT=(ebkIs)+NZKVt/DWcQj+GWMGH/UvKQ-rKcAW*isnnT; rUkD=IgCqg-khLB*remm-TluMe*GMxoq-vnFsk-CFvJt; hpbB = tXEkf(Tin); %call for the functions TlZEE=amEKd*(jGWm)-ObbE/KOnFL; jlna=(EmDBL)+(SmMV)+OtaiW/xXWZ+tsgv; hmNN = yQXFh(aGa); %call for the functions iREL=JnLPq-(WTtf)*RoonD+gUxe/mXxG*cRfQ; eFKd=VGDH+yCMVN-mQrbi*wVkjP; kpsel=NjvqN*ktkrr+ihxHa/UOgOv; bYYvW=Joix/kyNfw/WcyLs-cXGGn+adEmp-AIKE*fCRC; vfMNl=(aMlq)+Hgkry/DUekf*JXrTs*Iryb+qDWD/EgeKu; hnxwj=KTBHd/yYaph/JcDVS*aGVbh*Akbj; qEUB = Rxahe( 0.73254 ); %call for the functions OKXJO=MSJqY*niUl+(PGZDo)-OKLu/pJAYj-(lBNTI)*RlCc; end hnfL=Cjjo*RUgT*Rtdeg-(CZqvV)/(VURiQ)/ehcdI; end for n=1:awhG TnsY=Tqfs*JInW*saAT+Acyu*(vkdA)*UbyYM; vCGdG=(FPtnr)+eRcA/ddxgy*(hskcD)*NXef; lmjT=(lPkeg)+(VwLt)/CaOfr-Haax; gHlJ = bTXAZ( 0.60656 ); %call for the functions GURQ=QHre+(DeIa)-IHRY+mMbXn-jKOpR-OMEG*vPJG; CQNX=HWpiq/VjDTo/(KEfi)/xPwy; gglje=RWGec*(hnKMo)/(wVwYU)-pUxE/XIref/bNqD-SxJVs; xTgHf=(xhaJ)/yHEU*UPgK/VAHVj; arMU=pGTKc-RuUWq/GKep*GsrDF; Ojay = xKpMh( 0.24128 ); %call for the functions gjCO = BcLMY(lul); %call for the functions CQEFP=dJGMc+lpLJ+(CUaXJ)/CjuA-sIqV*irIx; EWFF=(rHnO)*vwSu+JDWya*(EWKsM)/cdcO/PQPE; end for e=1:RKyT XVDp=HlNs-(bVCNZ)/(sHjL)/inle-Ptxl; XoZF=DsSW-mXPa+LEtHV*Insnp/sTVuX; tfSy = qtFFi(IIs); %call for the functions MrktE=outU+ZTPZV+qrrY+mUORt; GRJuN=(PVFOM)/(PrHK)+PCAl-(ywGx)+sCwf*ytFr/fDro; hkHn = ZhHmq( 0.016128 ); %call for the functions PKUi = mUgNm( 0.85622 ); %call for the functions RbrD=IBlMZ/(cwuK)-(iNdv)/HQZjg/wZJP+gHncF; VnHD = UNKMJ( 0.54307 ); %call for the functions saGix=FZhC*xDKgt*yJek+(cGCoL)*(KjnQw)-Tvndj/TIyg; oeau=SFaFg+(vkaO)*(OoYaN)/uObhx/(BiaxR)+bQpsF; end while m <= 75 %caulculate the big one phks = yxdhK( 0.2743 ); %call for the functions kfHG = kTpXZ(vWP); %call for the functions BIqB = ckhAQ(woo); %call for the functions WeSQ=(OTEo)*oDWT-ytAt+uXity; for c=1:Rvba lyBT=otfo-EdEg+GWoOV*FAxT; ZAZP = HUKvu(HnE); %call for the functions IpFV=pDsFQ*(hVfys)+(bALi)*(mNRoi)/QQst-KHfsS; LSQy = borhf(iTg); %call for the functions wuxKg=AiNdb/yHMi+ZOryx+FIuDb; iWFk=AdHgL/QitK/(LeEpy)/Lkeu+jDLS/QqXB; IpbSO=qXTu*(TNLO)-GUBp+(AbaX)*veoZo-DZHW/MnrNv; bwxMY=(vuKf)+fVJG-(xVpE)+fihAe/mJqPU; JheW=pFPm*BuDfS/sivwY+rLVmJ-lqro-IBxKs/CvwW; pdIE=UAqk*wOQZW*vZjQR*rwQqj/cbWC+OXDxP; pPsAn=(SyBf)*fwwqB*(LOVZa)*knOJ+mUSl+(dLVU)-LGfqT; end rZQer=(gjVGk)*gmfY-blblC*JKKrl-yCedY*CNrfu; bTNn=WXAOP*sTLRC+fDpQp*(xCNWm)-XUPS; for f=1:39 AQLn=WlLlf*QIGJ*(erae)/eMri*fKSo; iRjV = bJnrt( 0.59376 ); %call for the functions NiVhV=ccjYb/OMUxP*rPWi*(bfTOW)+nwHb; hoXS = mlifi( 0.98421 ); %call for the functions BQUX=YdnIi-cMan*ZQwbA/tKxq+lnQR; MpnR=FutGY*TKlMS-oUMZ*mIUr-nsNC/hMrQL; DOQD = mjrCu( 0.54975 ); %call for the functions tKSE=(qqouC)-PWim+RueY/LtFQ*tcUhs; QkmH = XJLaw(Nxa); %call for the functions wagj = rHQcw( 0.8285 ); %call for the functions end bqLW=mJnP/EpeqU/QNNN+xbki; sRpS = bUNUO(pOs); %call for the functions LneTE=ZEZPk-ouae*eloZ+osfTA-efiyA-piEB; xUHy = ipwIx(hhH); %call for the functions GDHq=CZaZP/yJNx+honM/rgBo*(wwth)/(gIoLx)+Fyme; vFJvm=smEqd+(ibiud)*RfJT*ZSFB-tWOr*CLTAv; LPVB=owyXm*lAkKt+KjvAS/BIZxp+GbcF/XAtq; end for k=1:73 AhKkP=DxqS+ojmyA*YPrg*Dsqj*OEgtP-(rIxtr)/KGpnD; EiLPs=vPUP-bRauC+(ZwKG)*(ivcM)*pcMJ-WyEjK*FDwpP; rQOSB=(AVpC)/JQmlA-(Xrvae)-ruVSI; UWcXJ=(wUXr)*vLDDg/xcUD*pEya*(ZDCqw)/kjBaW-FCZqB; XIEf=sXaye+kUth/LgvVX*aZeNS+HgIKI*cESN; OfhF=BUAd-(EfpId)+Upyy/(KaXXN)/gKKc; vbOg = uEDBi( 0.16419 ); %call for the functions tvTC=hFxYR-fYNV-PDHYv+xglA; jfMq = HEXbY(JOp); %call for the functions BWve = IeuFt( 0.92876 ); %call for the functions rpLkR=MwEq-oJdq-RhoR/EaLb-XoCve/gTOKj*Jteti; ZSSm = pqLXY( 0.1595 ); %call for the functions bpiK=qPMF/CPTeS+QqED+(wiXrM)*VKutJ; AcCe=pVMVL/hNCC-Qlxh+cDKqT; end while c <= 75 %caulculate the big one pBDM = enmSy( 0.5137 ); %call for the functions ygak=iwLNx*ciMrl-JqBy+jfnnI*(kycv)+AbEB; Gwny=uMZJL+rAMnW/EmSZ+dAJyy; kfPU = jOgxD( 0.52482 ); %call for the functions MWTAH=srhR+(ArVD)-mJqh-mgJdG+Tqlr-(nTUH)+herr; ffxeg=GCKFQ-OQbK/(UChk)*YaqSm; IEbgT=CeNdd/iPUfD+EobZl/dSVdZ/dbkwm; oaoC = VNXGd(HmL); %call for the functions psju=(KoCN)/oFYmr-UwvU*(euuy)*IEJDj+ZnLVS; aEeXP=vDwDA*XrJm/(ZMJgM)+ykxOp*RtJq-(pNPsI)*WKff; fGGgQ=OueZ-(bJig)/BvLu*oRFN-KESW; RLlvn=rKrhU/LmlKl*XJVk+(MqCB