Changes between Version 54 and Version 55 of PolynomialExpansion


Ignore:
Timestamp:
01/26/16 23:31:26 (10 years ago)
Author:
sili
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • PolynomialExpansion

    v54 v55  
    3838  alpha,,1,, = ((b,,0,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,))^2^ + (b,,1,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,))^2^) / (p[0](a,,0,,p[0]+a,,1,,p[1]) + p[1](a,,1,,p[0]+a,,2,,p[1]))(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))^2^
    3939
    40 Step 8: r,,2,, = r,,1,, - alpha,,1,,*Ap,,1,,  ("has to expand for cancellation")
    41   Ap,,1,, = a,,0,,p[0]+a,,1,,p[1] + a,,1,,p[0]+a,,2,,p[1]
     40Step 8: r,,2,, = r,,1,, - alpha,,1,,*Ap,,1,,  ('''Need to expand for cancellation''')
    4241
    43   r[0] = 0
     42  r[0] = (b,,0,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - alpha,,1,,*(a,,0,,p[0]+a,,1,,p[1])
    4443
    45   r[1] = 0
     44  r[1] = (b,,1,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - alpha,,1,,*(a,,1,,p[0]+a,,2,,p[1])
    4645
    4746Step 9: x,,2,, = x,,1,, + alpha,,1,,*p,,1,,