Guest User

Untitled

a guest
Mar 18th, 2018
93
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.95 KB | None | 0 0
  1. ClearAll[Sum1];
  2. Attributes[Sum1] = {HoldFirst};
  3. Sum1[f_, {ks_List, ksummax_Integer?Positive}] :=
  4. With[{perms =
  5. Flatten[Permutations /@
  6. IntegerPartitions[#, {Length@ks}, Range[0, #]] & /@
  7. Range[2, ksummax], 2]}, Total[f /. (Thread[ks -> #] & /@ perms)]];
  8.  
  9. [CapitalPhi][dmax_, ksummax_][n_, s_,
  10. t_List] := (1/6*
  11. Sum[t[[i + 1]] t[[j + 1]] t[[n - i - j + 1]], {i, 0, n}, {j, 0,
  12. n - i}] +
  13. Sum1[Sum[
  14. No[d][Sequence @@ Table[k[i], {i, 2, n}]]*
  15. Product[t[[i + 1]]^k[i]/k[i]!, {i, 2, n}]*E^(d t[[2]])*s^d, {d,
  16. 0, dmax}], {k[#] & /@ Range[2, n], ksummax}])
  17.  
  18. In[218]:= [CapitalPhi][2, 4][3, s, Array[t, 4, 0]] // Simplify
  19.  
  20. Out[218]= 1/24 (4 t[1]^3 + 24 t[0] t[1] t[2] + 12 t[0]^2 t[3] +
  21. 12 t[3]^2 No[0][0, 2] + 4 t[3]^3 No[0][0, 3] +
  22. t[3]^4 No[0][0, 4] + 24 t[2] t[3] No[0][1, 1] +
  23. 12 t[2] t[3]^2 No[0][1, 2] + 4 t[2] t[3]^3 No[0][1, 3] +
  24. 12 t[2]^2 No[0][2, 0] + 12 t[2]^2 t[3] No[0][2, 1] +
  25. 6 t[2]^2 t[3]^2 No[0][2, 2] + 4 t[2]^3 No[0][3, 0] +
  26. 4 t[2]^3 t[3] No[0][3, 1] + t[2]^4 No[0][4, 0] +
  27. 12 E^t[1] s t[3]^2 No[1][0, 2] + 4 E^t[1] s t[3]^3 No[1][0, 3] +
  28. E^t[1] s t[3]^4 No[1][0, 4] + 24 E^t[1] s t[2] t[3] No[1][1, 1] +
  29. 12 E^t[1] s t[2] t[3]^2 No[1][1, 2] +
  30. 4 E^t[1] s t[2] t[3]^3 No[1][1, 3] +
  31. 12 E^t[1] s t[2]^2 No[1][2, 0] +
  32. 12 E^t[1] s t[2]^2 t[3] No[1][2, 1] +
  33. 6 E^t[1] s t[2]^2 t[3]^2 No[1][2, 2] +
  34. 4 E^t[1] s t[2]^3 No[1][3, 0] +
  35. 4 E^t[1] s t[2]^3 t[3] No[1][3, 1] + E^t[1] s t[2]^4 No[1][4, 0] +
  36. 12 E^(2 t[1]) s^2 t[3]^2 No[2][0, 2] +
  37. 4 E^(2 t[1]) s^2 t[3]^3 No[2][0, 3] +
  38. E^(2 t[1]) s^2 t[3]^4 No[2][0, 4] +
  39. 24 E^(2 t[1]) s^2 t[2] t[3] No[2][1, 1] +
  40. 12 E^(2 t[1]) s^2 t[2] t[3]^2 No[2][1, 2] +
  41. 4 E^(2 t[1]) s^2 t[2] t[3]^3 No[2][1, 3] +
  42. 12 E^(2 t[1]) s^2 t[2]^2 No[2][2, 0] +
  43. 12 E^(2 t[1]) s^2 t[2]^2 t[3] No[2][2, 1] +
  44. 6 E^(2 t[1]) s^2 t[2]^2 t[3]^2 No[2][2, 2] +
  45. 4 E^(2 t[1]) s^2 t[2]^3 No[2][3, 0] +
  46. 4 E^(2 t[1]) s^2 t[2]^3 t[3] No[2][3, 1] +
  47. E^(2 t[1]) s^2 t[2]^4 No[2][4, 0])
  48.  
  49. LHS = 2 Composition[D[#, t[1]] &, D[#, t[2]] &,
  50. D[#, t[3]] &] @[CapitalPhi][2, 4][3, s, Array[t, 4, 0]] //
  51. Simplify // Expand;
  52. RHS = Composition[D[#, t[1]] &, D[#, t[1]] &, D[#, t[1]] &,
  53. D[#, t[2] ] &, D[#, t[2]] &,
  54. D[#, t[2] ] &]@[CapitalPhi][3, 6][3, s, Array[t, 4, 0]] -
  55. Composition[D[#, t[1] ] &, D[#, t[1] ] &,
  56. D[#, t[2] ] &]@[CapitalPhi][2, 4][3, s, Array[t, 4, 0]]*
  57. Composition[D[#, t[1] ] &, D[#, t[2] ] &,
  58. D[#, t[2] ] &]@[CapitalPhi][2, 4][3, s, Array[t, 4, 0]] //
  59. Simplify // Expand;
  60.  
  61. coeffLHS = Flatten@CoefficientList[LHS, {Sequence @@ Array[t, 4, 0], E^t[1], s}];
  62. coeffRHS = Flatten@CoefficientList[RHS, {Sequence @@ Array[t, 4, 0], E^t[1], s}];
  63.  
  64. lim = Min[Length@coeffLHS, Length@coeffRHS];
  65. varsLHS = Variables[coeffLHS[[1 ;; lim]]];
  66. varsRHS = Variables[coeffRHS[[1 ;; lim]]];
  67. Solve[LogicalExpand[coeffLHS[[1 ;; lim]] == coeffRHS[[1 ;; lim]]],
  68. Union[varsLHS, varsRHS]]
Add Comment
Please, Sign In to add comment