<html><body><div style="color:#000; background-color:#fff; font-family:lucida console, sans-serif;font-size:10pt"><div id="yiv4159861690"><div><div style="color:#000;background-color:#fff;font-family:lucida console, sans-serif;font-size:10pt;"><div id="yiv4159861690yui_3_7_2_49_1377529034720_56"><span id="yiv4159861690yui_3_7_2_49_1377529034720_65">Hi Everyone,</span></div><div style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;" id="yiv4159861690yui_3_7_2_49_1377529034720_67"><br><span id="yiv4159861690yui_3_7_2_49_1377529034720_65"></span></div><div id="yiv4159861690yui_3_7_2_49_1377529034720_72" style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span id="yiv4159861690yui_3_7_2_49_1377529034720_65">I seek an assistance to convert the following program to MATLAB. It was originally written in
FORTRAN and now in Scilab.</span></div><div id="yiv4159861690yui_3_7_2_49_1377529034720_78" style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;
font-style:normal;"><br><span id="yiv4159861690yui_3_7_2_49_1377529034720_65"></span></div><div id="yiv4159861690yui_3_7_2_49_1377529034720_81" style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span id="yiv4159861690yui_3_7_2_49_1377529034720_65">Thanks as I anticipate your kind response/asssistance.<br><br></span><pre style="font-family:Monospaced;font-style:normal;font-size:12.0;"><span style="color:rgb(100,174,100);font-style:italic;">//</span>
<span style="color:rgb(100,174,100);font-style:italic;">//Right hand side of ODE</span>
<span style="color:rgb(100,174,100);font-style:italic;">//</span>
<span style="color:rgb(176,24,19);">function</span> <span style="color:rgb(131,67,16);font-weight:bold;">y</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(0,0,0);text-decoration:underline;">D</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">t</span><span style="color:rgb(0,0,0);">, </span><span style="color:rgb(131,67,16);font-weight:bold;">u</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(0,0,0);">V</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(131,67,16);font-weight:bold;">u</span><span style="color:rgb(74,85,219);">(</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);">H</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(131,67,16);font-weight:bold;">u</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(0,0,0);">VDOT</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">0.0107</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">H</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(188,143,143);">0.00205</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">V</span><span style="color:rgb(92,92,92);">^</span><span style="color:rgb(188,143,143);">2</span>
<span style="color:rgb(0,0,0);">HDOT</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">0.311</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(188,143,143);">0.0624</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">V</span>
<span style="color:rgb(131,67,16);font-weight:bold;">y</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(74,85,219);">[</span><span style="color:rgb(0,0,0);">VDOT</span><span style="color:rgb(0,0,0);">;</span><span style="color:rgb(0,0,0);">HDOT</span><span style="color:rgb(74,85,219);">]</span>
<span style="color:rgb(176,24,19);">endfunction</span>
<span style="color:rgb(0,0,0);">t0</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">0</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(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(255,170,0);">:</span><span style="color:rgb(188,143,143);">1200</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(100,174,100);font-style:italic;">//</span>
<span style="color:rgb(100,174,100);font-style:italic;">//First initial condition</span>
<span style="color:rgb(100,174,100);font-style:italic;">//</span>
<span style="color:rgb(0,0,0);">V1</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">3.40</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(0,0,0);">H1</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">2.05</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">y0</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(74,85,219);">[</span><span style="color:rgb(0,0,0);">V1</span><span style="color:rgb(0,0,0);">;</span><span style="color:rgb(0,0,0);">H1</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">mprintf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(188,143,143);">INITIAL CONDITIONS: V= %f, H=%f</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">V1</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">H1</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(0,0,0);">u</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(50,185,185);">ode</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">y0</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">t0</span><span style="color:rgb(0,0,0);">,</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);text-decoration:underline;">D</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">mprintf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(188,143,143);"> TIME V H</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(50,185,185);">disp</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">[</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);">,</span><span style="color:rgb(0,0,0);">u</span><span style="color:rgb(92,92,92);">'</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">scf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1001</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">plot2d</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);">u</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(255,170,0);">:</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(50,185,185);">xtitle</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">Gravity-Flow Tank Simulation V1,H1</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">Process time, s</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">h, ft</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">scf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1001</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">plot2d</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);">u</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(255,170,0);">:</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(100,174,100);font-style:italic;">//</span>
<span style="color:rgb(100,174,100);font-style:italic;">//Second initial condition</span>
<span style="color:rgb(100,174,100);font-style:italic;">//</span>
<span style="color:rgb(0,0,0);">V2</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">2.50</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(0,0,0);">H2</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">1.2</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">y0</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(74,85,219);">[</span><span style="color:rgb(0,0,0);">V2</span><span style="color:rgb(0,0,0);">;</span><span style="color:rgb(0,0,0);">H2</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">mprintf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(188,143,143);">INITIAL CONDITIONS: V= %f, H=%f</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">V2</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">H2</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(0,0,0);">u</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(50,185,185);">ode</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">y0</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">t0</span><span style="color:rgb(0,0,0);">,</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);text-decoration:underline;">D</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">mprintf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(188,143,143);"> TIME V H</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(50,185,185);">disp</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">[</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);">,</span><span style="color:rgb(0,0,0);">u</span><span style="color:rgb(92,92,92);">'</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">scf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1002</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">plot2d</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);">u</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(255,170,0);">:</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(50,185,185);">xtitle</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">Gravity-Flow Tank Simulation V2,H2</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">Process time, s</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">h, ft</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">scf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1002</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">plot2d</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);">u</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(255,170,0);">:</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(100,174,100);font-style:italic;">//end</span></pre><br></div><div id="yiv4159861690yui_3_7_2_49_1377529034720_85" style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><br><span id="yiv4159861690yui_3_7_2_49_1377529034720_65"></span></div><div id="yiv4159861690yui_3_7_2_49_1377529034720_88" style="color:rgb(0, 0, 0);font-size:13.3333px;font-family:lucida console, sans-serif;background-color:transparent;font-style:normal;"><span id="yiv4159861690yui_3_7_2_49_1377529034720_65"></span></div><div id="yiv4159861690yui_3_7_2_49_1377529034720_58"> </div><div id="yiv4159861690yui_3_7_2_49_1377529034720_60">ABDULJELEEL A. OSUNKUNLE<br>SCIENTIFIC
OFFICER<br>PROJECTS DEVELOPMENT INSTITUTE, P.M.B. 01609, EMENE INDUSTRIAL LAYOUT, ENUGU, NIGERIA.<br>+2347058408440.</div></div></div></div></div></body></html>