Changes between Version 77 and Version 78 of PolynomialExpansion


Ignore:
Timestamp:
01/27/16 17:27:58 (10 years ago)
Author:
sili
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • PolynomialExpansion

    v77 v78  
    5757Step 8: r,,2,, = r,,1,, - alpha,,1,,*Ap,,1,,  ('''Need to expand for cancellation''')
    5858
    59   r,,2,,[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,,)) - (a,,0,,p[0]+a,,1,,p[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,,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]))(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))^2^ = 0
     59  r,,2,,[0] = r,,1,,[0] - alpha,,1,,*(A[0][0]*p1[0]+A[0][1]*p1[1])
     60               =
     61   m[b,,0,,m - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,)][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])] - (a,,0,,p,,1,,[0]+a,,1,,p,,1,,[1]){[b,,0,,m - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,)]^2^ + [b,,1,,m - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,)]^2^}\\
     62   ------------------------------------------------------------------------------------------------------------\\
     63   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])]
    6064
    6165  r,,2,,[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,,)) - (a,,1,,p[0]+a,,2,,p[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,,1,,[0](a,,0,,p[0]+a,,1,,p,,1,,[1]) + p,,1,,[1](a,,1,,p,,1,,[0]+a,,2,,p,,1,,[1]))(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))^2^ = 0