Guest User

Untitled

a guest
Jan 16th, 2018
104
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.42 KB | None | 0 0
  1. deq = {pr[x] + a y'[x] + b pr'[x] - c pr''[x] == 0,
  2. e pt[x] + e2 pr'[x] + l pt'[x] + b2 pr[x]/x - dd pt''[x] == 0,
  3. v0 pr'[x] + v0 pr[x]/x - dd y'[x]/x - dd y''[x] == 0, y[L] == y0,
  4. pt[L] == pt0, pr[L] == pr0, y'[L] == v0 pr[L]/dd,
  5. pr'[L] == (pr[L] + v0 y[L])/(2 dd),
  6. pt'[L] == v0 (b2 pr[L] + e pt[L])/(2 dd)} /.
  7. Thread[Rule[{a, b, c, e, e2, l, b2, dd, L, pt0, pr0, v0},
  8. ConstantArray[1, 12]]];
  9.  
  10. DSolve[deq, {pr, pt, y}, x]
Add Comment
Please, Sign In to add comment