# How to solve this function?

1 view (last 30 days)
Ramesh on 5 Feb 2024
Edited: Matt J on 5 Feb 2024
function rate=ptMech(C,T,par)
% Obtain the parameter values
A=par.A;
E=par.E;
b=par.b;
R=8.314;
% Compute k
k=A*exp(-E/(R*T));
% Compute the result
rate=k*C^b;
end
Matt J on 5 Feb 2024
Edited: Matt J on 5 Feb 2024
What do you mean "solve" it? A function is not an equaiton.

William Rose on 5 Feb 2024
I am not sure what you mean by "how do you solve this function".
You have posted a function.
par=struct('A',1,'E',2,'b',3); % define par
T=273; C=1; % define T, C
rate=ptMech(C,T,par) % call function ptMech()
rate = 0.9991
function rate=ptMech(C,T,par)
% Obtain the parameter values
A=par.A;
E=par.E;
b=par.b;
R=8.314;
% Compute k
k=A*exp(-E/(R*T));
% Compute the result
rate=k*C^b;
end
If you call the function with appropriate parameters, as shown above, it will return an answer.
Do you want to fix some of the parameters, and then find the value of other parameters that will make ptMech equal zero? You can see that C=0 and A=0 (either one, or both) will yield ptMech=0.