Info

This question is closed. Reopen it to edit or answer.

Hi, what does it mean 'Not Enough Input Arguements'? Even though i have provided the values, but i am still unable to find the reason for error.

1 view (last 30 days)
function [dy r s] = getrate(t,y,param) % function for rate equations
getname; s = zeros(param.n_species,param.n_rxn); r = zeros(param.n_rxn,1); %disp(s);
%yc is a cell array for concentrations yc = num2cell(y); [D2, M, D2M, R, D2R, AKT_P, AKT, GSK3, GSK3_P]= deal(yc{:}); %disp(y);
% ---> D2 r(1,1) = param.k(1,1); s(param.nD2,1) = 1;
% D2 ---> r(2,1) = param.k(2,1)*D2; s(param.nnD2,2) = -1;
% ----> M r(3,1) = param.k(3,1); s(param.nM,3) = 1;
% M ----> r(4,1) = param.k(4,1)*M; s(param.nM,4) = -1;
% D2 + M ---> D2M r(5,1) = param.k(5,1)*D2*M; s(param.nD2M,5) = 1; s(param.nD2,5) = -1; s(param.nM,5) = -1;
% D2M ----> r(6,1) = param.k(6,1)*D2M; s(param.nD2M,6) = -1;
% D2 + R ----> D2R r(7,1) = param.k(7,1)*D2*R; s(param.nD2R,7) = 1; s(param.nD2,7) = -1; s(param.nR,7) = -1;
% D2R ----> r(8,1) = param.k(8,1)*D2R; s(param.nD2R,8) = -1;
% ----> AKT_P r(9,1) = param.k(9,1); s(param.nAKT_P,9) = 1;
% AKT_P ----> r(10,1) = param.k(10,1)*AKT_P; s(param.nAKT_P,10) = -1;
% AKT_P + [D2M] ----> AKT r(11,1) = (param.k(12,1)*D2M)/(D2M + param.k(11,1)); s(param.nAKT,11) = 1; s(param.nAKT_P,11) = -1; s(param.nD2M,11) = -1;
% the reaction that will happen when the drug is absent; k(13,1) = 0; % AKT -----> AKT_P r(12,1) = param.k(14,1)*AKT; s(param.nAKT,12) = -1; s(param.nAKT_P,12) = 1;
% AKT + [D2R] ----> AKT_P; the reaction that will occur in the presence of Risperidone , K(13,1) = 1 r(13,1) = (param.k(13,1)*D2R)/(D2R + param.k(14,1)); s(param.nAKT_P,13) = 1; s(param.nAKT,13) = -1; s(param.nD2R,13) = -1;
% ----> GSK3 r(14,1) = param.k(15,1); s(param.nGSK3,14) = 1;
% GSK3 + [AKT_P] ----> GSK3_P r(15,1) = (param.k(17,1)*AKT_P)/(AKT_P + param.k(16,1)); s(param.nGSK3_P,15) = 1; s(param.nGSK3,15) = -1; s(param.nAKT_P,15) = -1;
% GSK3 ----> GSK3_P r(16,1) = param.k(18,1)*GSK3_P; s(param.nGSK3_P,16) = 1; s(param.nGSK3,16) = -1;
% GSK3_P -----> r(17,1) = param.k(17,1)*GSK3_P; s(param.nGSK3_P,17) = -1;
dy = s*r; %keyboard; return;

Answers (0)

This question is closed.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!