I always get the error I have written above. I decided to try a working code which you can find in the book Essentials of Chemical Reaction Engineering:International Edition(H.Scott Fogler). This code is also resembling the code that I wrote. And the error was the same when I copied it to my MatLab(7.12.0, R2011a)
Code:
function dYfuncvecdt = ODEfun(t,Yfuncvec);
Rs= Yfuncvec(1);
RTi= Yfuncvec(2);
Cs= Yfuncvec(3);
LTi= Yfuncvec(4);
% Explicit equations
kf = 7.2e7;
kr = .34;
keR = .03;
keC = .165;
krec = .058;
kdeg = .0022;
Vs = .013;
fR = .5;
fL = .5;
L = 10^-8;
Navgo = 6.02*10^23;
kfP = 0;
% Differential equations
dRsdt = -kf*Rs*L+kr*Cs-keC*Rs+(krec*(1-fR)*RTi)+Vs;
dRTidt = keR*Rs+keC*Cs-(krec*(1-fR)*RTi)-kdeg*fR*RTi;
dCsdt = kf*L*Rs-kr*Cs-keC*Cs;
dLTidt = keC*Cs+kfP*L*Navgo-(krec*(1-fL)*LTi)-kdeg*fR*RTi;
dYfuncvecdt=[dRsdt;dRTidt;dCsdt;dLTidt];
Error:
??? function dYfuncvecdt = ODEfun(t,Yfuncvec);
|
Error: Function definitions are not permitted in this context.
That was the first thing i've tried and I got this error. After that i tried to copy it directly. I didn't understand what is going on.
??? Input argument "Yfuncvec" is undefined.
Error in ==> ODEfun at 2
Rs= Yfuncvec(1);
Error in ==> run at 57
evalin('caller', [s ';']);
The errors say that you did not define the input variable Yfunvec when you call ODEfun
You have to call the function in the matlab command line as Kenth said
ODEfun(0.5,rand(4,1))
Is this code works you have to replace the inputs by your parameters
ODEfun(yourT,[yourRs,yourRti,yourCs,yourLTi])
No, it does not work with this definition.
??? function dYfuncvecdt = ODEfun(0.5,rand(4,1));
|
Error: Function definitions are not permitted in this context.
You must save the code into a m file.
After that in the matlab console try the code
ODEfun(0.5,rand(4,1))
The working folder must be the same where you save the m file
If this code works you can change this random input for your analysis parameters
ODEfun(yourT,[yourRs,yourRti,yourCs,yourLTi])
I appreciate every answer you wrote. I'm a new matlab user and i make too many mistakes. I hope i will learn quickly.
Hi Cansu,
attached are the m files which for the sample set of equations you put up previously.
To get the output of the model simulations you need to run the ODEsolution.m file.
Please beware that the equations presented seem to be different from the ones in the paper you are working with. therefore, you should expect different results.
Here is the ODEsolution.m file. The two files must be put in the same directory
You cant run function which has an input argument. Yfuncvec is input argument in your function. You have to call this function from outside and give this input argument.
- Badges
- Science topic