Guest User

Untitled

a guest
May 26th, 2018
130
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.17 KB | None | 0 0
  1. [Beta][M] = 0;
  2. [Psi] = JacobiDN[u, m] Exp[I [Omega] x];
  3. [Alpha][j_] := [Alpha][
  4. j] = -I Integrate[([Psi] Conjugate[[Beta][j]] -
  5. Conjugate[[Psi]] [Beta][j]), u] + 2^(j - 1) f[j]
  6. [Beta][j_] := [Beta][j] = [Psi] [Alpha][j + 1] +
  7. I/2 D[[Beta][j + 1], u]
  8.  
  9. $Assumptions =
  10. x [Element] Reals && [Omega] [Element] Reals && m >= 0 && m < 1 &&
  11. u [Element] Reals && M [Element] Integers &&
  12. f[_] [Element] Reals
  13.  
  14. In[76]:= q = M - 4 ;
  15. FullSimplify[TrigToExp[ComplexExpand[[Beta][q]]]]
  16.  
  17. Out[77]= 2^(-4 + M) E^(
  18. I x [Omega]) ((f[-3 + M] - m f[-1 + M]) JacobiDN[u, m] +
  19. I m (-f[-2 + M] + m f[M]) JacobiCN[u, m] JacobiSN[u, m])
  20.  
  21. In[79]:= FullSimplify[TrigToExp[ComplexExpand[[Alpha][q]]]]
  22.  
  23. Out[79]= 2^(-5 +
  24. M) (f[-4 + M] +
  25. m (-m f[M] - 2 (f[-2 + M] - m f[M]) JacobiCN[u, m]^2))
  26.  
  27. In[81]:= B =
  28. 2^(-4 + M) E^(
  29. I x [Omega]) ((f[-3 + M] - m f[-1 + M]) JacobiDN[u, m] +
  30. I m (-f[-2 + M] + m f[M]) JacobiCN[u, m] JacobiSN[u, m]);
  31. Integrand = FullSimplify[Conjugate[B] [Psi] - Conjugate[[Psi]] B];
  32. A = FullSimplify[- I Integrate[Integrand, u] + 2^(M - 4 - 1) f[M - 4]]
  33.  
  34. Out[83]= 2^(-5 +
  35. M) (f[-4 + M] + 2 m (-f[-2 + M] + m f[M]) JacobiCN[u, m]^2)
Add Comment
Please, Sign In to add comment