<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><meta name="ProgId" content="Word.Document"><meta name="Generator" content="Microsoft Word 12"><meta name="Originator" content="Microsoft Word 12"><link rel="File-List" href="file:///C:%5CDOCUME%7E1%5CADMINI%7E1%5CLOCALS%7E1%5CTemp%5Cmsohtmlclip1%5C01%5Cclip_filelist.xml"><link rel="themeData" href="file:///C:%5CDOCUME%7E1%5CADMINI%7E1%5CLOCALS%7E1%5CTemp%5Cmsohtmlclip1%5C01%5Cclip_themedata.thmx"><link rel="colorSchemeMapping" href="file:///C:%5CDOCUME%7E1%5CADMINI%7E1%5CLOCALS%7E1%5CTemp%5Cmsohtmlclip1%5C01%5Cclip_colorschememapping.xml"><style>
<!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;
mso-font-charset:0;
mso-generic-font-family:roman;
mso-font-pitch:variable;
mso-font-signature:-1610611985 1107304683 0 0 159 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;
mso-font-charset:0;
mso-generic-font-family:swiss;
mso-font-pitch:variable;
mso-font-signature:-1610611985 1073750139 0 0 159 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{mso-style-unhide:no;
mso-style-qformat:yes;
mso-style-parent:"";
margin:0cm;
margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-ascii-font-family:Calibri;
mso-ascii-theme-font:minor-latin;
mso-fareast-font-family:Calibri;
mso-fareast-theme-font:minor-latin;
mso-hansi-font-family:Calibri;
mso-hansi-theme-font:minor-latin;
mso-bidi-font-family:"Times New Roman";
mso-bidi-theme-font:minor-bidi;
mso-fareast-language:EN-US;}
span.EmailStyle15
{mso-style-type:personal;
mso-style-noshow:yes;
mso-style-unhide:no;
mso-ansi-font-size:11.0pt;
mso-bidi-font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-ascii-font-family:Calibri;
mso-ascii-theme-font:minor-latin;
mso-fareast-font-family:Calibri;
mso-fareast-theme-font:minor-latin;
mso-hansi-font-family:Calibri;
mso-hansi-theme-font:minor-latin;
mso-bidi-font-family:"Times New Roman";
mso-bidi-theme-font:minor-bidi;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
mso-default-props:yes;
mso-ascii-font-family:Calibri;
mso-ascii-theme-font:minor-latin;
mso-fareast-font-family:Calibri;
mso-fareast-theme-font:minor-latin;
mso-hansi-font-family:Calibri;
mso-hansi-theme-font:minor-latin;
mso-bidi-font-family:"Times New Roman";
mso-bidi-theme-font:minor-bidi;
mso-fareast-language:EN-US;}
@page Section1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;
mso-header-margin:36.0pt;
mso-footer-margin:36.0pt;
mso-paper-source:0;}
div.Section1
{page:Section1;}
-->
</style>
<p class="MsoNormal">Hi all<br>
<br>
I am busy programming a penicillin growth inside a fed batch reactor model in
SCILAB.<br>
<br>
Now I have 9 ordinary differential equations where i am supposed to get 6
graphs from the 9 ordinary differential equations.<br>
<br>
Now when using only the 6 main ordinary differential equations, then i am
getting six graphs but the graphs is wrong because the other 3 ordinary
differential equations is going to have a impact on the six graphs. As soon as
i am programming the other 3 ordinary differential equations my model doesnt
want to run graph 7 till 9.<br>
<br>
it runs but only shows title name and y label and x label.<br>
<br>
the error i am getting says the following.<br>
<br>
"<span style="color: red;">plot2d(t,y(7,:))<br>
!--error 21 <br>
Invalid index.<br>
t line 58 of function f called by : <br>
y=ode(x0, 0, t, f);"<br>
</span><br>
i am not familiar with SCILAB so can someone assist me and explaining if i
should programm the last 3 ordinary differential equations differently?</p><p class="MsoNormal"><br></p><p class="MsoNormal">the 9 equations are as follows.</p><p class="MsoNormal"><br></p><p class="MsoNormal">dx(1) = (((mu)*(x(1)))-((x(1))/(x(6)))*(dx_6)); // biomass concentration X<br>
dx(2) = z*(mu*x(1)-(F*x(1))/x(6))+QQ; // hydrogen ion concentration H+<br>dx(3) = mupp*(x(1))-(K)*(x(3))-((x(3))/x(6))*(dx_6); // Penicilin concentration P<br>dx(4) = -mu/Yxs*x(1)-mupp/Yps*x(1)-mx*x(1)+Fsf/x(6)-x(4)/x(6)*dx_6; // Substrate concentration S<br>
dx(5) = -mu/Yxo*x(1)-mupp/Ypo*x(1)-mo*x(1)+kla*(cll-x(5))-x(5)/x(6)*dx_6; //dissolved oxygen concentration Cl<br>dx(6) = F+Fab+Floss; //Culture Volume V<br>dx(7) = ((rq1(x(1))*x(6))+(rq2*(x(1))*(x(6)))); // Heat of reaction Qrxn</p>
<p class="MsoNormal">dx(8) = ((((F/sf)(Tf-(x(8))))+(1/(x(6)*pcp)))*(QT)); // Temperature T</p><p class="MsoNormal">dx(9) = (a1*(dx_1)+a2*x(1)+a3); // CO2 evolution, CO2</p><p class="MsoNormal"><br></p><p class="MsoNormal">
As you can see that eqaution 7 till 9 has some effects on other ordinary differential eqautions.</p><p class="MsoNormal"><br></p><p class="MsoNormal"><br></p>