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

index exceeds the number of array elements

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';    endee=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','k','MarkerFaceColor','r','MarkerSize',5)        figure(114),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)        figure(115),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)            elseif L==12        figure(121),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')        figure(122),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')        figure(123),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)        figure(124),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)        figure(125),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)            end;    end;``

``for xx=1:11``
``    for xx=1:10``