How to solve this function?

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]

What do you mean "solve" it? A function is not an equaiton.

Answer 1

[William Rose]

@Ramesh,

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.