<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: times new roman,new york,times,serif; font-size: 12pt; color: #000000'><pre style="font-family:Monospaced;font-style:normal;font-size:14.0;"><span style="color:rgb(0,0,0);">I think the following code fullfil your wish...<br><br><br>n</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">0</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">20</span> <span style="color:rgb(0,0,0);">;</span><span style="color:rgb(1,168,1);font-style:italic;">// series in the expansion</span>
<span style="color:rgb(0,0,0);">B</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">1.5</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(1,168,1);font-style:italic;">//beta factor</span>
<span style="color:rgb(0,0,0);">Tm</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">1300</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(188,143,143);">274.15</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(1,168,1);font-style:italic;">//base-lithosphere temperature [C]</span>
<span style="color:rgb(0,0,0);">tau</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">62.8</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(1,168,1);font-style:italic;">//lithosphere cooling thermal decay constant [Ma]</span>
<span style="color:rgb(0,0,0);">a</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">125000</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(1,168,1);font-style:italic;">//equilibrium lithospheric thickness [m]</span>
<span style="color:rgb(0,0,0);">alpha</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">3.28e-5</span> <span style="color:rgb(0,0,0);">;</span><span style="color:rgb(1,168,1);font-style:italic;">//thermal expansion coefficient K^-1</span>
<span style="color:rgb(1,168,1);font-style:italic;">// Tm needs to be converted to Kelvin - multiply by 274.15</span>
<span style="color:rgb(0,0,0);">rho</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">3300</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">t</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">0</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">5</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">150</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(1,168,1);font-style:italic;">//lithosphere age [Ma]</span>
<span style="color:rgb(0,0,0);">G</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">6.67e-11</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(1,168,1);font-style:italic;">// LTGA with increasing age of oceanic lithosphere</span>
<span style="color:rgb(0,0,0);">c</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">8</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">G</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">alpha</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">rho</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">a</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">Tm</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(218,112,214);">%pi</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">m</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">0</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">20</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">'</span><span style="color:rgb(0,0,0);">;</span><span style="color:rgb(0,0,0);">m1</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">m</span><span style="color:rgb(92,92,92);">+</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">g</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(0,0,0);">c</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(50,185,185);">sum</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">B</span><span style="color:rgb(92,92,92);">./</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">m1</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(188,143,143);">3</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(218,112,214);">%pi</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">.*</span><span style="color:rgb(50,185,185);">sin</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">m1</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(218,112,214);">%pi</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(0,0,0);">B</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(50,185,185);">ones</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">t</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">.*</span><span style="color:rgb(50,185,185);">exp</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">m1</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">t</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(0,0,0);">tau</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">clf</span><span style="color:rgb(0,0,0);">;</span><span style="color:rgb(174,92,176);text-decoration:underline;">plot</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">t</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">g</span><span style="color:rgb(74,85,219);">)</span></pre><br>Serge Steer<br><hr id="zwchr"><blockquote style="border-left:2px solid rgb(16, 16, 255);margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><b>De: </b>"Lester Anderson" <arctica1963@gmail.com><br><b>À: </b>"International users mailing list for Scilab." <users@lists.scilab.org><br><b>Envoyé: </b>Mardi 11 Décembre 2012 08:16:01<br><b>Objet: </b>[Scilab-users] Summation query in Scilab<br><br>Hello,<br><br>I have an equation structured like in the attached, however I am unsure of the correct syntax to compute in Scilab.<br><br>One can calculate the first term of the expansion but obvioulsy it needs to go to a large number (approximating infinity eg 20).<br>
I am sure this is something very straightforward!<br><br>Any ideas welcome.<br><br>Thanks for the help<br><br>Lester<br><br>att: image of equation, test code<br>
<br>_______________________________________________<br>users mailing list<br>users@lists.scilab.org<br>http://lists.scilab.org/mailman/listinfo/users<br></blockquote><br></div></body></html>