void calculateTma(int limit)
{
   int i,j,k;
   double FullLength = 2.0*HalfLength+1.0;
     
   for (i=limit; i>=0; i--)
   {
      double sum  = (HalfLength+1)*iMA(NULL,0,1,0,MODE_SMA,Price,i);
      double sumw = (HalfLength+1);
      for(j=1, k=HalfLength; j<=HalfLength; j++, k--)
      {
         sum  += k*iMA(NULL,0,1,0,MODE_SMA,Price,i+j);
         sumw += k;

         if (j<=i)
         {
            sum  += k*iMA(NULL,0,1,0,MODE_SMA,Price,i-j);
            sumw += k;
         }
      }
      tmBuffer[i] = sum/sumw;
  
      double diff = iMA(NULL,0,1,0,MODE_SMA,Price,i)-tmBuffer[i];
      if (i> (Bars-HalfLength-1)) continue;
      if (i==(Bars-HalfLength-1))
      {
         upBuffer[i] = tmBuffer[i];
         dnBuffer[i] = tmBuffer[i];
         if (diff>=0)
            {
               wuBuffer[i] = MathPow(diff,2);
               wdBuffer[i] = 0;
            }
         else
            {               
               wdBuffer[i] = MathPow(diff,2);
               wuBuffer[i] = 0;
            }                  
         continue;
      }
      if(diff>=0)
      {
          wuBuffer[i] = (wuBuffer[i+1]*(FullLength-1)+MathPow(diff,2))/FullLength;
          wdBuffer[i] =  wdBuffer[i+1]*(FullLength-1)/FullLength;
      }
      else
      {
           wdBuffer[i] = (wdBuffer[i+1]*(FullLength-1)+MathPow(diff,2))/FullLength;
           wuBuffer[i] =  wuBuffer[i+1]*(FullLength-1)/FullLength;
      }
      upBuffer[i] = tmBuffer[i] + BandsDeviations*MathSqrt(wuBuffer[i]);
      dnBuffer[i] = tmBuffer[i] - BandsDeviations*MathSqrt(wdBuffer[i]);
   }
}

iMA(NULL,0,1,0,MODE_SMA,Price,i+j);