Changes between Version 90 and Version 91 of PolynomialExpansion


Ignore:
Timestamp:
01/28/16 11:06:07 (10 years ago)
Author:
sili
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • PolynomialExpansion

    v90 v91  
    7171
    7272Step 9: x,,2,, = x,,1,, + alpha,,1,,*p,,1,,
     73  If not expand:
     74
    7375  x,,2,,[0] = x,,1,,[0] + alpha,,1,,*p,,1,,[0]
    7476
     
    8587   m^2^[p,,1,,[0](a,,0,,p,,1,,[0]+a,,1,,p,,1,,[1]) + p,,1,,[1](a,,1,,p,,1,,[0]+a,,2,,p,,1,,[1])]
    8688
     89  if expand and simplified:
     90
     91  x,,2,,[0] = (a,,2,,b,,0,, - a,,1,,b,,1,,) / (a,,0,,a,,2,, - a,,1,,^2^)
     92
     93  x,,2,,[1] = (a,,0,,b,,1,, - a,,1,,b,,0,,) / (a,,0,,a,,2,, - a,,1,,^2^)
     94
    8795assertion:   ('''Need to expand for cancellation''')
    8896  bncg[0] = A[0][0]*x,,2,,[0] + A[0][1]*x,,2,,[1] = a,,0,,(a,,2,,b,,0,, - a,,1,,b,,1,,) / (a,,0,,a,,2,, - a,,1,,^2^) + a,,1,,(-a,,1,,b,,0,, + a,,0,,b,,1,,) / (a,,0,,a,,2,, - a,,1,,^2^) = b,,0,,(a,,0,,a,,2,, - a,,1,,^2^) / (a,,0,,a,,2,, - a,,1,,^2^) = b,,0,,