matrix dimension must agree and the size is different on left and right side. I need to iterate to observe cooling. Please help !

R_in =[0.02;1.0]; %total E Resistance

Keff = 1.5e-3; % W/K

rho = 10e-6; % ohm.meter

Sab = 200e-6; % Volt per Kelvin

N = [75];

% Ti = 200; % intiial temp

% Hhx = 0.004; % heat transfer Coef forced Air circulation

Ta = [50;10] ; % cooling side or being pumped in Kelvin

alpha =220e-6;% V/Keff seebecKeff coeff

I_in = [2:30];

V_in = [2:30];

W = V_in.*I_in;

T_pump_cold = 200.15; % final required temperature

T_sink = 173.15;

Tb = 298.15; % hot side initial temperature

DT = T_pump_cold - T_sink;

dx = 0.04; % in m

Q_c = 2*N*((Sab*T_pump_cold.*I_in) - 2*N*((Keff*rho)./R_in)*(DT)-0.5*(I_in.^2).*R_in)

COP = Q_c/W;

A_cr = 0.00002; % in m2

t(1) =1;

T(1)=T_pump_cold ;

Time2=(1:t:12); hours

for j = (1:t:1000)

[T(j+1)] = ((dx.*Q_c)./Keff.*Time2.*A_cr)+T(j);

end