Guest User

Untitled

a guest
Mar 18th, 2018
79
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.98 KB | None | 0 0
  1. Remove["Global`*"];
  2. cr = 1.8*10^7/0.5*omega^2*Exp[-omega];
  3. r55 = 2/Pi*Integrate[cr*Cos[omega*t], {omega, 0, 10}]
  4.  
  5. (*r55 = 1/(Pi (1.` + t^2)^3)
  6. 2 (7.200000000000001`*^7 -
  7. 2.1600000000000006`*^8 t^2 + (-199396.49151683348` -
  8. 317073.10946119425` t^2 - 130751.79771595637` t^4) Cos[
  9. 10 t] + (238622.0308316204` t + 388986.59820497024` t^3 +
  10. 163439.7471449455` t^5) Sin[10 t]) *)
  11.  
  12.  
  13. Remove["Global`*"];
  14. omega = 1.2;
  15.  
  16. i1 = 25640;
  17. ia = 4*10^6;
  18. c2 = 1.4513*10^5;
  19. t5 = 16000*Cos[omega*t];
  20. tpto = 2.8*10^6*Theta'[t];
  21.  
  22. r55 = 1/(Pi (1.` + t^2)^3)
  23. 2 (7.200000000000001`*^7 -
  24. 2.1600000000000006`*^8 t^2 + (-199396.49151683348` -
  25. 317073.10946119425` t^2 - 130751.79771595637` t^4) Cos[
  26. 10 t] + (238622.0308316204` t + 388986.59820497024` t^3 +
  27. 163439.7471449455` t^5) Sin[10 t]) /. {t -> t - Tau};
  28.  
  29. s = NDSolve[{(i1 + ia)*Theta''[t] +
  30. Integrate[r55*Theta'[Tau], {Tau, -Infinity, t}] +
  31. c2*Theta[t] == t5 + tpto, Theta'[0] == 1, Theta[0] ==
  32. 0}, Theta, {t, 0, 10}]
Add Comment
Please, Sign In to add comment