Mathematical Modelling for simulation . Error in ODE45.

kemo()

Array indices must be positive integers or logical values.

Error in solution>kemo/kk (line 51)
dxdt_3=((-teta*em*L)/(teta+I))+j*(T/(k+T))*L-q*L*T+(r1*N+r2*C)*T((u*(L^2)*C*I)/(kappa+I))+((pie*L*I)/(gi+I))-kl*(1-exp(delta*M))*L;

Error in odearguments (line 92)
f0 = ode(t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ode45 (line 107)
odearguments(odeIsFuncHandle,odeTreatAsMFile, solver_name, ode, tspan, y0, options, varargin);

Error in solution>kemo (line 22)
[t,x]=ode45(@kk,t,[initial_t;initial_n;initial_l;initial_c;initial_i;initial_m]);

function kemo 
t=0:0.1:200; 

a=(4.31)*10^(-1); b=(1.02)*10^(-9); c=(2.9077)*10^(-13);delta=1.8328; 
d=4.5*10^(-8);ef=0.111;f=(4)*10^(-1);p=(2.794)*10^(13);pn=(6.68)*10^(-2);gn=(2.5036)*10^5;
em =9*10^(-3);teta=(2.5036)*10^(-3);q=(3.422)*10^(10);r1=(2.9077)*10^(-11);r2=(5.8467)*10^(13);pie=2.971;gi=(2.5036)*10^3;u=(4.417)*10^(-14);
kappa=(2.5036)*10^3;j=(1.245)*10^(2);k=(2.019)*10^3;miu=11.7427;phi=(2.38405)*10^(-2);
kt=0.9;kn=0.0675;kl=0.0486;kc=0.034;gamma=0.5199;

vl=(1.77)*10^10;vm=2.3869;vi=2.7859*10^6;
alphabeta=(2.25)*10^(9);beta=(6.3)*10^(-3); 

initial_t=10^7; 
initial_n=2.5*10^8;  
initial_l=5*10^4;  
initial_c=2.25*10^9;  
initial_i=1073;  
initial_m=0; 

[t,x]=ode45(@kk,t,[initial_t;initial_n;initial_l;initial_c;initial_i;initial_m]);  


figure (1)  
semilogy(t,x(:,6),'b','LineWidth',2);  
title('Grafik M terhadap t','Fontsize',15);  
xlabel('waktu(hari)','Fontsize',15); 
ylabel('M(t)(mg/liter)','Fontsize',15);  
axis([0 30 10^0 10^12]); 
grid on 

figure(2)  
semilogy(t,x(:,1),'b',t,x(:,2),'g',t,x(:,3),'r',t,x(:,4),'c',t,x(:,5),'y','LineWidth',1); 
xlabel('waktu(hari)','Fontsize',15); 
% ylabel('','Fontsize',15); 
legend('Sel tumor T(t)','Sel NK N(t)','Sel T-CD8+L(t)','Sellimfosit C(t)','IL-2 I(t)')  
axis([0 200 10^0 10^12]); 
grid on 

    function dxdt=kk(t,x)  
        T=x(1); 
        N=x(2); 
        L=x(3); 
        C=x(4); 
        I=x(5); 
        M=x(6); 
        dxdt_1=a*T*(1-b*T)-c*N*T-d*L*T-kt*(1-exp(-delta*M))*T; 
        dxdt_2=f*((ef)*C-N)-p*N*T+((pn*N*I)/(gn+I))-kn*(1-exp(delta*M))*N;
        
        dxdt_3=((-teta*em*L)/(teta+I))+j*(T/(k+T))*L-q*L*T+(r1*N+r2*C)*T((u*(L^2)*C*I)/(kappa+I))+((pie*L*I)/(gi+I))-kl*(1-exp(delta*M))*L;
        dxdt_4=beta*(alphabeta-C)-kc*(1-exp(-delta*M))*C; 
        dxdt_5=-miu*I+phi*C;  
        dxdt_6=-gamma*M+vm; 
        
        dxdt=[dxdt_1;dxdt_2;dxdt_3;dxdt_4;dxdt_5;dxdt_6]; 
    end 
end

Can help me to fix that ?

Error in ode45 (line 115)

odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

Error in kemo (line 20)

[t,x]=ode45(@kk,t,[initial_t;initial_n;initial_l;initial_c;initial_i;initial_m]);

[John D'Errico]

When you get an error, show the COMPLETE error message, thus everything in red. Merely saying there is an error, even if you show the line it happens in, is not any real help.

If you want help, then make it easy to get help.

[Torsten]

dxdt_3=((-teta*em*L)/(teta+I))+j*(T/(k+T))*L-q*L*T+(r1*N+r2*C)*T((u*(L^2)*C*I)/(kappa+I))+((pie*L*I)/(gi+I))-kl*(1-exp(delta*M))*L;

here

There is an arithmetic operator missing after the "T".