<rss xmlns:itunes="http://www.itunes.com/dtds/podcast-1.0.dtd" version="2.0">
  <channel>
    <title>MaplePrimes - answers and comments on Question, Matlab To Maple</title>
    <link>http://www.mapleprimes.com/questions/37168-Matlab-To-Maple</link>
    <language>en-us</language>
    <copyright>2026 Maplesoft, A Division of Waterloo Maple Inc.</copyright>
    <generator>Maplesoft Document System</generator>
    <lastBuildDate>Tue, 09 Jun 2026 16:20:06 GMT</lastBuildDate>
    <pubDate>Tue, 09 Jun 2026 16:20:06 GMT</pubDate>
    <itunes:subtitle />
    <itunes:summary />
    <description>The latest answers and comments added to the Question, Matlab To Maple</description>
    <image>
      <url>http://www.mapleprimes.com/images/mapleprimeswhite.jpg</url>
      <title>MaplePrimes - answers and comments on Question, Matlab To Maple</title>
      <link>http://www.mapleprimes.com/questions/37168-Matlab-To-Maple</link>
    </image>
    <item>
      <title>Explanation, please?</title>
      <link>http://www.mapleprimes.com/questions/37168-Matlab-To-Maple?ref=Feed:MaplePrimes:Matlab To Maple:Comments#answer65333</link>
      <itunes:summary>&lt;p&gt;Rather than just trying to understand what your MATLAB code is supposed to do and then determining how to best implement that algorithm in Maple, could you just tell us what problem you are trying to solve?&lt;/p&gt;
&lt;p&gt;Doug&lt;/p&gt;
&lt;pre&gt;
---------------------------------------------------------------------
Douglas B. Meade  &amp;lt;&amp;gt;&amp;lt;
Math, USC, Columbia, SC 29208  E-mail: mailto:meade@math.sc.edu
Phone:  (803) 777-6183         URL:    http://www.math.sc.edu&lt;/pre&gt;</itunes:summary>
      <description>&lt;p&gt;Rather than just trying to understand what your MATLAB code is supposed to do and then determining how to best implement that algorithm in Maple, could you just tell us what problem you are trying to solve?&lt;/p&gt;
&lt;p&gt;Doug&lt;/p&gt;
&lt;pre&gt;
---------------------------------------------------------------------
Douglas B. Meade  &amp;lt;&amp;gt;&amp;lt;
Math, USC, Columbia, SC 29208  E-mail: mailto:meade@math.sc.edu
Phone:  (803) 777-6183         URL:    http://www.math.sc.edu&lt;/pre&gt;</description>
      <guid>65333</guid>
      <pubDate>Wed, 17 Jun 2009 21:21:51 Z</pubDate>
      <itunes:author>Doug Meade</itunes:author>
      <author>Doug Meade</author>
    </item>
    <item>
      <title>dear Sir,
 
     Iam</title>
      <link>http://www.mapleprimes.com/questions/37168-Matlab-To-Maple?ref=Feed:MaplePrimes:Matlab To Maple:Comments#answer65334</link>
      <itunes:summary>&lt;p&gt;dear Sir,&lt;/p&gt;
&lt;p&gt;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; Iam waiting yours replay.&lt;/p&gt;
&lt;p&gt;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;Thanks.&lt;/p&gt;
&lt;p&gt;Ahmed&lt;/p&gt;</itunes:summary>
      <description>&lt;p&gt;dear Sir,&lt;/p&gt;
&lt;p&gt;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; Iam waiting yours replay.&lt;/p&gt;
&lt;p&gt;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;Thanks.&lt;/p&gt;
&lt;p&gt;Ahmed&lt;/p&gt;</description>
      <guid>65334</guid>
      <pubDate>Sun, 21 Jun 2009 18:45:29 Z</pubDate>
      <itunes:author>ahmedluss</itunes:author>
      <author>ahmedluss</author>
    </item>
    <item>
      <title>More explanation</title>
      <link>http://www.mapleprimes.com/questions/37168-Matlab-To-Maple?ref=Feed:MaplePrimes:Matlab To Maple:Comments#answer65335</link>
      <itunes:summary>&lt;p&gt;Yes, more explanation.&amp;nbsp; Ahmedluss, can you&amp;nbsp;explain a little more what your code does.&amp;nbsp; Some of us do not have Matlab.&amp;nbsp; Yes Maple can export and convert code to Matlab but it's apparent and unfortunate Matlab doesn't&amp;nbsp;convert to Maple.&amp;nbsp; My first quick guess at what your code does, is some sort of insect&amp;nbsp;growth&amp;nbsp;population scenario.&amp;nbsp;&lt;/p&gt;</itunes:summary>
      <description>&lt;p&gt;Yes, more explanation.&amp;nbsp; Ahmedluss, can you&amp;nbsp;explain a little more what your code does.&amp;nbsp; Some of us do not have Matlab.&amp;nbsp; Yes Maple can export and convert code to Matlab but it's apparent and unfortunate Matlab doesn't&amp;nbsp;convert to Maple.&amp;nbsp; My first quick guess at what your code does, is some sort of insect&amp;nbsp;growth&amp;nbsp;population scenario.&amp;nbsp;&lt;/p&gt;</description>
      <guid>65335</guid>
      <pubDate>Mon, 22 Jun 2009 04:26:31 Z</pubDate>
      <itunes:author>Christopher2222</itunes:author>
      <author>Christopher2222</author>
    </item>
    <item>
      <title>Dear,
   I want to used</title>
      <link>http://www.mapleprimes.com/questions/37168-Matlab-To-Maple?ref=Feed:MaplePrimes:Matlab To Maple:Comments#answer65336</link>
      <itunes:summary>&lt;p&gt;Dear,&lt;/p&gt;
&lt;p&gt;&amp;nbsp;&amp;nbsp; I want to used these codes for optimiziation by using&amp;nbsp; Partical Swarm Optimization (POS) algorithm,&lt;/p&gt;
&lt;p&gt;&lt;b&gt;&lt;font face="Times-Bold" size="2"&gt;
&lt;p align="left"&gt;Algorithm:&lt;/p&gt;
&lt;font face="Times-Bold" size="2"&gt;
&lt;p align="left"&gt;&amp;nbsp;&lt;/p&gt;
&lt;/font&gt;&lt;/font&gt;&lt;/b&gt;&lt;/p&gt;
&lt;p align="left"&gt;&lt;font face="Times-Roman" size="2"&gt;Original PSO.&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;D &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;dimensions in the search space.&lt;/font&gt;&lt;b&gt;&lt;font face="Times-Bold" size="2"&gt;loop&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;
&lt;p align="left"&gt;3: For each particle, evaluate the desired optimization fitness function in&lt;/p&gt;
&lt;p align="left"&gt;4: Compare particle&amp;rsquo;s fitness evaluation with its&lt;/p&gt;
&lt;p align="left"&gt;5: Identify the particle in the neighborhood with the best success so far, and assign itsindex to the variable g.&lt;/p&gt;
&lt;/font&gt;&lt;i&gt;&lt;font face="Times-Italic" size="2"&gt;D &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;variables.&lt;/font&gt;&lt;i&gt;&lt;font face="Times-Italic" size="2"&gt;pbest&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;. If current value is better than &lt;/font&gt;&lt;i&gt;&lt;font face="Times-Italic" size="2"&gt;pbest&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;, then set &lt;/font&gt;&lt;i&gt;&lt;font face="Times-Italic" size="2"&gt;pbest&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;equal to the current value, and &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;p&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;equal to the current location &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;in &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;D&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;-dimensional space.&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;
&lt;p align="left"&gt;6: Change the velocity and position of the particle according to the following equation(see notes below):&lt;/p&gt;
&lt;/font&gt;&lt;font lang="JA" face="MTEX" size="2"&gt;
&lt;p align="left"&gt;&amp;nbsp;&lt;/p&gt;
&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;
&lt;p align="left"&gt;&amp;nbsp;&lt;/p&gt;
&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;v&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;larr; &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;v&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;+ &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;U(&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;0&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;,&amp;phi;&lt;/font&gt;&lt;font face="Times-Roman" size="1"&gt;1&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;)&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;otimes;&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;( &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;p&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;minus; &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="MTMI" size="2"&gt;)&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;+ &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;U(&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;0&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;,&amp;phi;&lt;/font&gt;&lt;font face="Times-Roman" size="1"&gt;2&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;)&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;otimes;&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;( &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;p&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;g &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;minus; &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="MTMI" size="2"&gt;),&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;
&lt;p align="left"&gt;&amp;nbsp;&lt;/p&gt;
&lt;p align="left"&gt;7: If a criterion is met (usually a sufficiently good fitness or a maximum number ofiterations), exit loop.&lt;/p&gt;
&lt;p&gt;8:&lt;/p&gt;
&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;larr; &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;+ &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;v&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="MTMI" size="2"&gt;.&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;(1)&lt;/font&gt;&lt;b&gt;&lt;font face="Times-Bold" size="2"&gt;end loop&lt;/font&gt;
&lt;p&gt;&amp;nbsp;the Codes in Matlab &lt;br /&gt;
&lt;br /&gt;
&amp;nbsp;&lt;/p&gt;
&lt;p&gt;% THE&lt;br /&gt;
% MATLAB&lt;/p&gt;
&lt;p&gt;function [xmin, fxmin, iter] = PSO&lt;/p&gt;
&lt;p&gt;%%%&lt;br /&gt;
success = 0;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;PopSize = 20;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;br /&gt;
MaxIt = 5000;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;br /&gt;
iter = 0;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;fevals = 0;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
maxw = 1.2;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;minw = 0.1;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;weveryit = floor(0.75*MaxIt);&amp;nbsp;&lt;br /&gt;
c1 = 0.5;&amp;nbsp;&lt;br /&gt;
c2 = 0.5;&amp;nbsp;&lt;br /&gt;
inertdec = (maxw-minw)/weveryit;&lt;br /&gt;
w = maxw;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
f = &amp;ldquo;DeJong&amp;rdquo;;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
dim = 2;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;upbnd = 5;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;wbnd = &amp;ndash;5;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
GM = 0;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
ErrGoal = 1e&amp;ndash;3;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;%&amp;nbsp;Initializing&amp;nbsp;&lt;br /&gt;
popul = rand(dim, PopSize)*(upbnd-lwbnd) + lwbnd;&lt;br /&gt;
vel = rand(dim, PopSize);&lt;/p&gt;
&lt;p&gt;%Evaluate&amp;nbsp;&lt;br /&gt;
for i = 1:PopSize,&lt;br /&gt;
fpopul(i) = feval(f, popul(:,i));&lt;br /&gt;
fevals = fevals + 1;&lt;br /&gt;
end&lt;/p&gt;
&lt;p&gt;%&amp;nbsp;Initializing&amp;nbsp;&lt;br /&gt;
% values&lt;br /&gt;
bestpos = popul;&lt;br /&gt;
fbestpos = fpopul;&lt;/p&gt;
&lt;p&gt;% Finding &lt;br /&gt;
[fbestpart,g] = min(fpopul);&lt;br /&gt;
lastbpf = fbestpart;&lt;/p&gt;
&lt;p&gt;% SWARM EVOLUTION LOOP ? START ?&lt;br /&gt;
while (success == 0) &amp;amp; (iter &amp;lt; MaxIt),&lt;br /&gt;
iter = iter + 1;&lt;br /&gt;
% Update the value of the inertia weight w&lt;br /&gt;
if (iter&amp;lt;=weveryit)&lt;br /&gt;
w = maxw &amp;ndash; (iter&amp;ndash;1)?inertdec;&lt;br /&gt;
end&lt;/p&gt;
&lt;p&gt;% VELOCITY UPDATE&lt;br /&gt;
for i=1:PopSize,&lt;br /&gt;
A(:,i) = bestpos(:,g);&lt;br /&gt;
end&lt;br /&gt;
R1 = rand(dim, PopSize);&lt;br /&gt;
R2 = rand(dim, PopSize);&lt;br /&gt;
vel = w?vel + c1?R1.?(bestpos-popul) + c2?R2.?(A-popul);&lt;/p&gt;
&lt;p&gt;%&lt;/p&gt;
&lt;p&gt;popul = popul + vel;&lt;/p&gt;
&lt;p&gt;% Evaluate the new swarm&lt;br /&gt;
for i = 1:PopSize,&lt;br /&gt;
fpopul(i) = feval(f,popul(:, i));&lt;br /&gt;
fevals = fevals + 1;&lt;br /&gt;
end&lt;/p&gt;
&lt;p&gt;% Updating&amp;nbsp;&lt;br /&gt;
changeColumns = fpopul &amp;lt; fbestpos;&lt;br /&gt;
fbestpos = fbestpos.*(~changeColumns) + fpopul.*changeColumns;&lt;br /&gt;
bestpos(:, find(changeColumns)) = popul(:, find(changeColumns));&lt;/p&gt;
&lt;p&gt;% Updating&lt;br /&gt;
[fbestpart, g] = min(fbestpos);&lt;/p&gt;
&lt;p&gt;% Checking&amp;nbsp;&lt;br /&gt;
%if&amp;nbsp;&lt;br /&gt;
if abs(fbestpart-GM) &amp;lt;= ErrGoal&lt;br /&gt;
success = 1;&lt;br /&gt;
else&lt;br /&gt;
lastbpf = fbestpart;&lt;br /&gt;
end&lt;br /&gt;
end&lt;br /&gt;
%&amp;nbsp; LOOP ? END ?&lt;br /&gt;
% Output &lt;br /&gt;
xmin = popul(:,g);&lt;br /&gt;
fxmin = fbestpos(g);&lt;/p&gt;
&lt;p&gt;function DeJong=DeJong(x)&lt;br /&gt;
DeJong = sum(x.?2);&lt;br /&gt;
&amp;nbsp;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
How can rewrite these codes&amp;nbsp;&amp;nbsp;in Maple because I used Maple only .&lt;br /&gt;
&lt;br /&gt;
I am waiting youres replay.&lt;/p&gt;
&lt;p&gt;Thank you so much&lt;/p&gt;
&lt;p&gt;Ahmed&lt;/p&gt;
&lt;p&gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;nbsp;&lt;/p&gt;
&lt;/b&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/b&gt;&lt;/i&gt;&lt;/p&gt;
&lt;p align="left"&gt;1: Initialize a population array of particles with random positions and velocities on&lt;/p&gt;
&lt;p align="left"&gt;2:&lt;/p&gt;</itunes:summary>
      <description>&lt;p&gt;Dear,&lt;/p&gt;
&lt;p&gt;&amp;nbsp;&amp;nbsp; I want to used these codes for optimiziation by using&amp;nbsp; Partical Swarm Optimization (POS) algorithm,&lt;/p&gt;
&lt;p&gt;&lt;b&gt;&lt;font face="Times-Bold" size="2"&gt;
&lt;p align="left"&gt;Algorithm:&lt;/p&gt;
&lt;font face="Times-Bold" size="2"&gt;
&lt;p align="left"&gt;&amp;nbsp;&lt;/p&gt;
&lt;/font&gt;&lt;/font&gt;&lt;/b&gt;&lt;/p&gt;
&lt;p align="left"&gt;&lt;font face="Times-Roman" size="2"&gt;Original PSO.&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;D &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;dimensions in the search space.&lt;/font&gt;&lt;b&gt;&lt;font face="Times-Bold" size="2"&gt;loop&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;
&lt;p align="left"&gt;3: For each particle, evaluate the desired optimization fitness function in&lt;/p&gt;
&lt;p align="left"&gt;4: Compare particle&amp;rsquo;s fitness evaluation with its&lt;/p&gt;
&lt;p align="left"&gt;5: Identify the particle in the neighborhood with the best success so far, and assign itsindex to the variable g.&lt;/p&gt;
&lt;/font&gt;&lt;i&gt;&lt;font face="Times-Italic" size="2"&gt;D &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;variables.&lt;/font&gt;&lt;i&gt;&lt;font face="Times-Italic" size="2"&gt;pbest&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;. If current value is better than &lt;/font&gt;&lt;i&gt;&lt;font face="Times-Italic" size="2"&gt;pbest&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;, then set &lt;/font&gt;&lt;i&gt;&lt;font face="Times-Italic" size="2"&gt;pbest&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;equal to the current value, and &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;p&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;equal to the current location &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;in &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;D&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;-dimensional space.&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;
&lt;p align="left"&gt;6: Change the velocity and position of the particle according to the following equation(see notes below):&lt;/p&gt;
&lt;/font&gt;&lt;font lang="JA" face="MTEX" size="2"&gt;
&lt;p align="left"&gt;&amp;nbsp;&lt;/p&gt;
&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;
&lt;p align="left"&gt;&amp;nbsp;&lt;/p&gt;
&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;v&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;larr; &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;v&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;+ &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;U(&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;0&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;,&amp;phi;&lt;/font&gt;&lt;font face="Times-Roman" size="1"&gt;1&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;)&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;otimes;&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;( &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;p&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;minus; &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="MTMI" size="2"&gt;)&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;+ &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;U(&lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;0&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;,&amp;phi;&lt;/font&gt;&lt;font face="Times-Roman" size="1"&gt;2&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;)&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;otimes;&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;( &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;p&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;g &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;minus; &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="MTMI" size="2"&gt;),&lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;
&lt;p align="left"&gt;&amp;nbsp;&lt;/p&gt;
&lt;p align="left"&gt;7: If a criterion is met (usually a sufficiently good fitness or a maximum number ofiterations), exit loop.&lt;/p&gt;
&lt;p&gt;8:&lt;/p&gt;
&lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;&amp;larr; &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;x&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font lang="JA" face="MTSYN" size="2"&gt;+ &lt;/font&gt;&lt;i&gt;&lt;font face="MTMI" size="2"&gt;v&lt;/font&gt;&lt;font face="MTMI" size="1"&gt;i &lt;/font&gt;&lt;font face="MTMI" size="2"&gt;.&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;/font&gt;&lt;font face="Times-Roman" size="2"&gt;(1)&lt;/font&gt;&lt;b&gt;&lt;font face="Times-Bold" size="2"&gt;end loop&lt;/font&gt;
&lt;p&gt;&amp;nbsp;the Codes in Matlab &lt;br /&gt;
&lt;br /&gt;
&amp;nbsp;&lt;/p&gt;
&lt;p&gt;% THE&lt;br /&gt;
% MATLAB&lt;/p&gt;
&lt;p&gt;function [xmin, fxmin, iter] = PSO&lt;/p&gt;
&lt;p&gt;%%%&lt;br /&gt;
success = 0;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;PopSize = 20;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;br /&gt;
MaxIt = 5000;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;br /&gt;
iter = 0;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;fevals = 0;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
maxw = 1.2;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;minw = 0.1;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;weveryit = floor(0.75*MaxIt);&amp;nbsp;&lt;br /&gt;
c1 = 0.5;&amp;nbsp;&lt;br /&gt;
c2 = 0.5;&amp;nbsp;&lt;br /&gt;
inertdec = (maxw-minw)/weveryit;&lt;br /&gt;
w = maxw;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
f = &amp;ldquo;DeJong&amp;rdquo;;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
dim = 2;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;upbnd = 5;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;wbnd = &amp;ndash;5;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
GM = 0;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;br /&gt;
ErrGoal = 1e&amp;ndash;3;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/p&gt;
&lt;p&gt;%&amp;nbsp;Initializing&amp;nbsp;&lt;br /&gt;
popul = rand(dim, PopSize)*(upbnd-lwbnd) + lwbnd;&lt;br /&gt;
vel = rand(dim, PopSize);&lt;/p&gt;
&lt;p&gt;%Evaluate&amp;nbsp;&lt;br /&gt;
for i = 1:PopSize,&lt;br /&gt;
fpopul(i) = feval(f, popul(:,i));&lt;br /&gt;
fevals = fevals + 1;&lt;br /&gt;
end&lt;/p&gt;
&lt;p&gt;%&amp;nbsp;Initializing&amp;nbsp;&lt;br /&gt;
% values&lt;br /&gt;
bestpos = popul;&lt;br /&gt;
fbestpos = fpopul;&lt;/p&gt;
&lt;p&gt;% Finding &lt;br /&gt;
[fbestpart,g] = min(fpopul);&lt;br /&gt;
lastbpf = fbestpart;&lt;/p&gt;
&lt;p&gt;% SWARM EVOLUTION LOOP ? START ?&lt;br /&gt;
while (success == 0) &amp;amp; (iter &amp;lt; MaxIt),&lt;br /&gt;
iter = iter + 1;&lt;br /&gt;
% Update the value of the inertia weight w&lt;br /&gt;
if (iter&amp;lt;=weveryit)&lt;br /&gt;
w = maxw &amp;ndash; (iter&amp;ndash;1)?inertdec;&lt;br /&gt;
end&lt;/p&gt;
&lt;p&gt;% VELOCITY UPDATE&lt;br /&gt;
for i=1:PopSize,&lt;br /&gt;
A(:,i) = bestpos(:,g);&lt;br /&gt;
end&lt;br /&gt;
R1 = rand(dim, PopSize);&lt;br /&gt;
R2 = rand(dim, PopSize);&lt;br /&gt;
vel = w?vel + c1?R1.?(bestpos-popul) + c2?R2.?(A-popul);&lt;/p&gt;
&lt;p&gt;%&lt;/p&gt;
&lt;p&gt;popul = popul + vel;&lt;/p&gt;
&lt;p&gt;% Evaluate the new swarm&lt;br /&gt;
for i = 1:PopSize,&lt;br /&gt;
fpopul(i) = feval(f,popul(:, i));&lt;br /&gt;
fevals = fevals + 1;&lt;br /&gt;
end&lt;/p&gt;
&lt;p&gt;% Updating&amp;nbsp;&lt;br /&gt;
changeColumns = fpopul &amp;lt; fbestpos;&lt;br /&gt;
fbestpos = fbestpos.*(~changeColumns) + fpopul.*changeColumns;&lt;br /&gt;
bestpos(:, find(changeColumns)) = popul(:, find(changeColumns));&lt;/p&gt;
&lt;p&gt;% Updating&lt;br /&gt;
[fbestpart, g] = min(fbestpos);&lt;/p&gt;
&lt;p&gt;% Checking&amp;nbsp;&lt;br /&gt;
%if&amp;nbsp;&lt;br /&gt;
if abs(fbestpart-GM) &amp;lt;= ErrGoal&lt;br /&gt;
success = 1;&lt;br /&gt;
else&lt;br /&gt;
lastbpf = fbestpart;&lt;br /&gt;
end&lt;br /&gt;
end&lt;br /&gt;
%&amp;nbsp; LOOP ? END ?&lt;br /&gt;
% Output &lt;br /&gt;
xmin = popul(:,g);&lt;br /&gt;
fxmin = fbestpos(g);&lt;/p&gt;
&lt;p&gt;function DeJong=DeJong(x)&lt;br /&gt;
DeJong = sum(x.?2);&lt;br /&gt;
&amp;nbsp;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
How can rewrite these codes&amp;nbsp;&amp;nbsp;in Maple because I used Maple only .&lt;br /&gt;
&lt;br /&gt;
I am waiting youres replay.&lt;/p&gt;
&lt;p&gt;Thank you so much&lt;/p&gt;
&lt;p&gt;Ahmed&lt;/p&gt;
&lt;p&gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;nbsp;&lt;/p&gt;
&lt;/b&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/i&gt;&lt;/b&gt;&lt;/i&gt;&lt;/p&gt;
&lt;p align="left"&gt;1: Initialize a population array of particles with random positions and velocities on&lt;/p&gt;
&lt;p align="left"&gt;2:&lt;/p&gt;</description>
      <guid>65336</guid>
      <pubDate>Mon, 22 Jun 2009 07:36:15 Z</pubDate>
      <itunes:author>ahmedluss</itunes:author>
      <author>ahmedluss</author>
    </item>
  </channel>
</rss>