i want to insert mulit resolution wavelet to A wavelet network pid code

resolution=(upper -lower)/(wl-1); for i=1:wl; if(i==1); translation(i)=lower; else translation(i)= translation(i-1)+resolution; end;

end; yy(k)=x(k);

ap(k)=0; for q=0:M; for i=1:wl; zm(i)=(t-translation(i)-q)/dilation(i); ap(k)=ap(k)+(zm(i)*exp(-0.5*zm(i)^2))*w(i)*yy(k)*ff(q+1);

                end;
            end;
 gama(k)=ap(k)/yy(k);
   app(k)=0;        
    for q=1:N;
          app(k)=app(k)+(ap(k)-q)*v*fb(q);
    alpha(k)=app(k);
    end;
            appp(k)=ap(k)+app(k);

            err=yy(k)-appp(k);

            for q=0:M;
                for i=1:wl;
                    zz(i)=zz(i)+(zm(i)*exp(-0.5*zm(i)^2))*ff(q+1);
                    ww(i)=ww(i)+err*yy(k)*zz(i);
                    zzz(i)=zzz(i)+w(i)*(1/dilation(i))*(zm(i)^2-1)*exp(-0.5*zm(i)^2)*ff(q+1);
                    bb(i)=bb(i)+err*yy(k)*zzz(i);
                    aa(i)=aa(i)+err*yy(k)*zzz(i)*((t-translation(i))/dilation(i));
                end;
            end;

            for q=0:M;
            for i=1:wl;

            cc(q+1)=cc(q+1)+(zm(i)*exp(-0.5*zm(i)^2))*w(i)*yy(k)*err;
           end;
            end;
           for q=1:N;
           dd(q)=dd(q)+(app(k)-q)*err*v;
            end;

    %%%%%%%%%%%%%%%%%%   appp(k) represents y^(k) %%%%%%%%%%%%%%%%%%%%%
     e(k)=y(k)-appp(k);
     error=r(k)-appp(k);
     E=E+err^2;

     if (k==1)
         P(k+1)=Pr;
         I(k+1)=In;
     else
         P(k+1)=P(k-1)+MP*e(k)*gama(k)*(er(k)-er(k-1));
         I(k+1)=I(k-1)+MI*e(k)*gama(k)*er(k);
     end;
     if (k==1)
         D(k+1)=Dr;
     elseif (k==2)
         D(k+1)=D(k-1)+MD*e(k)*gama(k)*(er(k)-2*er(k-1));

     else
       D(k+1)=D(k-1)+MD*e(k)*gama(k)*(er(k)-2*er(k-1)+er(k-2));
     end;
 %%%%%%%%%%%%%%%%%%%%%%%%%%%end to the t loop %%%%%%%%

     for i=1:wl;
            translation(i)=translation(i)+mb*bb(i);
            dilation(i)=dilation(i)+ma*aa(i);
            w(i)=w(i)+mw*ww(i);

     end;

        for q=0:M;
            ff(q+1)=ff(q+1)+mc*cc(q+1);
        end;
        for q=1:N;
            fb(q)=fb(q)+md*dd(q);
        end;

        E=0.5*E;
        end; 
        iteration=iteration+1
        ISE(iteration)=E

   end