MATLAB: Index exceeds the number of array elements (10), How to solve this error?

Question

Error in Main_Program (line 101)

Tg(xx)=(((a*Tw(xx)^2-b*Tw(xx)+c)*Tw(xx)+(Ar*Ta(xx)*hw)+Ar*Ts(xx)*(0.048*Ta(xx)-9))/((a*Tw(xx)^2-b*Tw(xx)+c)…

global N;
global ee;
dt=1;
tf=24;
nm=12;
t(1)=0;
np=(tf-t(1))/dt;
dt1=(t(1)+1)/dt;
Tw0(1)=0;
l=0.05;
a=0.02612;
b=15.76;
c=2392;
d=0.048;
e=3.8213;
w=0.0;
s=5.67*10^-8;
Ta1=[15 16.5 18 22.5 24 26 27 27 26.5 25.5 24];
Ta2=[19 23 25.5 28.5 29.5 30 32 32 31.5 31 30];
Ta3=[25.5 29.5 32 35 37 37 37.5 38 37 36 34];
Ta4=[29.5 33 35 37 38.5 38.5 38.5 39 38.5 38.5 38];
Ta5=[30 32.5 34.5 37 38.5 39.5 41 41.5 41 40 39.5];
Ta6=[33 35.5 37 39.5 41 41 42 42.5 42 41.5 40.5];
Ta7=[29.5 32.5 34 35.5 36.5 37 38 38 33 33 32.5];
Ta8=[23 25.5 28.5 30 31 32 33 33 33.5 33 32.5];
Ta9=[29 29 30 31 34 35.5 36.5 37 38 38 37];
Ta10=[33 34 35 38 39.5 40 40 39 38 37];
Ta11=[23.5 25 28 29.5 31 32 34 35 34 32.5 31];
Ta12=[20 23.5 25 26 29 31 31.5 30 30 30.5 30];
RR=[15;46;74;105;135;166;181;212;243;273;304;349];
sw=[20 23 24 25 28 28 25 25 25 25 23 20];
Z=1;
Z1=10;
V1=4;
Z0=0.03;
V=(log2(Z/Z0)/log2(Z1/Z0))*V1;
hw=2.8+3*V;
Aw=1;
Ag=1.46;
Ar=Ag/Aw;
for L=1:nm;
    N=RR(L);
GTOTAL=g_t();
if L==1
    disp(['Variation of temperature productivity and solar intensity'])
    disp(['of January'])
Ta=Ta1';
elseif L==2
    disp(['Variation of temperature productivity and solar intensity'])
    disp(['of February'])
Ta=Ta2';
elseif L==3
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of March'])
Ta=Ta3';
elseif L==4
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of April'])   
Ta=Ta4';
elseif L==5
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of May'])
Ta=Ta5';
elseif L==6
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of June']) 
Ta=Ta6';
elseif L==7
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of July']) 
Ta=Ta7';
elseif L==8
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of August'])   
Ta=Ta8';
elseif L==9
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of September']) 
Ta=Ta9';
elseif L==10
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of October'])
Ta=Ta10';
elseif L==11
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of November']) 
Ta=Ta11';
elseif L==12
     disp(['Variation of temperature productivity and solar intensity'])
     disp(['of December']) 
Ta=Ta12';    
end
ee=sw(L);
Tw(1)=ee;
Ts=0.0552*Ta.^(1.5);
for xx=1:11
    kk(xx)=xx+7;
    q(xx)=xx;
    %Ts=[10 10 9.9 9.5 8.6 8.8 9.5 10.4 11.6 12.4 13.5 15.6 15.7 15.5 14.2
    %14 13.1 12.2 11.5 11.2 11 10.3 10 10];
    Tg(xx)=(((a*Tw(xx)^2-b*Tw(xx)+c)*Tw(xx)+(Ar*Ta(xx)*hw)+Ar*Ts(xx)*(0.048*Ta(xx)-9))/((a*Tw(xx)^2-b*Tw(xx)+c)...
        +(Ar*hw)+Ar*(0.048*Ta(xx)-9))-e);
    Ti(xx)=Tw(xx)/2.0+Tg(xx)/2.0;
    K(xx)=0.0244+0.7673*10^-4*Ti(xx);
    Cw(xx)=999.2+0.1343*Ti(xx)+0.01*10^-4*Ti(xx)^2-6.758*10^-8*Ti(xx)^3;
    Ub(xx)=K(xx)/1;
    Eg=0.98;
    x=647.27-(Tw(xx)+273.15);
    a1=3.2437814;
    b1=5.86826*10^-3;
    c1=1.1702379*10^-8;
    d1=2.1878462*10^-3;
    xg(xx)=647.27-(Tg(xx)+273.15);
    Pw(xx)=165960.72*10^-(x*(a1+b1*x+c1*x^3)/((Tw(xx)+273.15)*(1+d1*x)));
    Pg(xx)=165960.72*10^(-(xg(xx)*(a1+b1*x+c1*xg(xx)^3))/((Tg(xx)+273.15)*(1+d1*x)));
    Pw1(xx)=(101300/760)*Pw(xx);
    Pg1(xx)=(101300/760)*Pg(xx);
    y=Tw();
    hfg(xx)=3044205.5-1679.1109*(Tw(xx)+273)-1.1425*(Tw(xx)+273)^2;
    hrwg(xx)=0.9*s*(Tw(xx)^2+Tg(xx)^2)*(Tw(xx)+Tg(xx));
    hcwg(xx)=0.884*((Tw(xx)-Tg(xx))+((Pw1(xx)-Pg1(xx))/(268900-Pw1(xx)))*Tw(xx))^(1/3);
    hewg(xx)=(9.15*10^-7*hcwg(xx)*(Pw1(xx)-Pg1(xx))*hfg(xx))/(Tw(xx)-Tg(xx));
    hrgs(xx)=Eg*s*(Tg(xx)^2+Ts(xx)^2)*(Tg(xx)+Ts(xx));
    hcgs(xx)=hw*(Tg(xx)-Ta(xx))/(Tg(xx)-Ts(xx));
    Ui(xx)=hrwg(xx)+hcwg(xx)+hewg(xx);
    Uo(xx)=hw*(Tg(xx)-Ta(xx))/(Tg(xx)-Ts(xx))+hrgs(xx);
    Ut1(xx)=(1/(Ui(xx)+1/(Ar*Uo(xx))));
    Ut(xx)=1/Ut1(xx);
    Eue(xx)=(hewg(xx)/Ut(xx))*Ui(xx)*(Tw(xx)-Tg(xx));
    Mw(xx)=(Eue(xx))*10^8/hfg(xx);
    p(xx)=Eue(xx);
    Eff(xx)=Mw(xx)*hfg(xx)*10^-6/GTOTAL(xx);
    Tw(xx+1)=(10+(Ta(xx)+Tw(xx))/2)+(Tw(xx)+(Eff(xx)*GTOTAL(xx))+(Ut(xx)*Ts(xx)))+(Ub(xx)*Ta(xx))/((Cw(xx)/Aw)+Ut(xx)+Ub(xx))*w;
    z11(xx)=Tw(xx);
    t(xx+1)=t(xx)+dt
    dx1b(xx)=Aw/Cw(xx)*(Eff(xx)*GTOTAL(xx)-Ut(xx)*(Tw(xx)-Ts(xx))-Ub(xx)*(Tw(xx)-Ta(xx)));
    Tw(xx+1)=Tw(xx)+dx1b(xx)*dt1;
    Tw1(xx)=(Tw(xx)*Tw0(1)+y(xx))*dt1;
    dx1e(xx)=Aw/Cw(xx)*(Eff(xx)*GTOTAL(xx)-Ut(xx)*(Tw(xx+1)-Ts(xx))-Ub(xx)*(Tw(xx+1)+Ta(xx)));
    dTw(xx)=(dx1b(xx)+dx1e(xx))/2;
    Tw(xx+1)=Tw(xx)+dTw(xx)*dt;
    Tw(xx+1)=Tw1(xx);
end
    format short;
    Tg1=Tg'
    Tw2=Tw1';
    dd=kk';
    Mw1=Mw'*10^-5;
    disp([' Time(hr) WatTemp GlasTemp ambienttemp SkyTemp Productivity x10^2 ']);
    disp([dd Tw2 Tg1 Ta Ts Mw1]);
    
    if L==1
        figure(11),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(12),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(13),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(14),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(15),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        
    elseif L==2
        figure(21),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(22),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(23),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(24),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(25),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        
    elseif L==3
        figure(31),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(32),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(33),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(34),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(35),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        
    elseif L==4
        figure(41),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(42),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(43),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(44),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(45),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        
    elseif L==5
        figure(51),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(52),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(53),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(54),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(55),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
    
    elseif L==6
        figure(61),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(62),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(63),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(64),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(65),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
      
    elseif L==7
        figure(71),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(72),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(73),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(74),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(75),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        
    elseif L==8
        figure(81),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(82),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(83),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(84),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(85),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        
    elseif L==9
        figure(91),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(92),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(93),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(94),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(95),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        
    elseif L==10
        figure(101),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(102),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(103),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(104),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        figure(105),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
        
    elseif L==11
        figure(111),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
        figure(112),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
        figure(113),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor'
            

            
            
            

    

        

	
    

        
Best Answer
    

    

								

								
Five lines above it, you initialise the for loop as: 
for xx=1:11
Perhaps: 
    for xx=1:10
will eliminate the error
						    
						
		
    
        

    

Related Solutions
MATLAB: Hello, I have three complex differential equations (dTb/dt , dTw/dt , dTg/dt) and I solved it using runge kutta method , matlab run successfully but all result on commend window is “NoN”,I’m sure from the equations.

								

								
Hey Muhammad, Is it possible for you to attach your three complex equations as they were given? The reason why I'm asking is that most people misrepresent mathematical expressions in their matlab code.
						    
						
MATLAB: I want to find the determinant for 5.5 matrix

								

								
>> syms s1 s2 s3 s4 s5 s6 s7 s8 s9 s10
>> M = [0 -s1 -s2 -s3 -s4; s1 0 -s5 -s6 -s4; s2 s5 0 -s8 -s9; s3 s6 s8 0 -s10; s4 s7 s9 s10 0];
>> det(M)
ans =
s1*s4^2*s8^2 - s1*s4*s7*s8^2 - s2*s4^2*s6*s8 + s3*s4^2*s5*s8 - s3^2*s4*s5*s9 - s2^2*s4*s6*s10 + s3^2*s5*s7*s9 + s2^2*s6*s7*s10 + s2*s3*s4*s5*s10 + s2*s3*s4*s6*s9 + s1*s2*s4*s8*s10 - s1*s3*s4*s8*s9 - s2*s3*s5*s7*s10 - s2*s3*s6*s7*s9 + s2*s4*s6*s7*s8 - s3*s4*s5*s7*s8 - s1*s2*s7*s8*s10 + s1*s3*s7*s8*s9
>> simplify(ans)
ans =
(s4 - s7)*(s3*s5 - s2*s6 + s1*s8)*(s2*s10 - s3*s9 + s4*s8)
						    
						
        
        

            
            
            
Related Question
Solving system of equations with parameters
How to join points within a plot with a line
How to use scatter with transparent dots
Fit non-linear scatter
How can a draw a line between two points
        
    

Answer 1

Five lines above it, you initialise the for loop as:

for xx=1:11

Perhaps:

    for xx=1:10

will eliminate the error