Guest User

Untitled

a guest
Jan 4th, 2018
69
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.46 KB | None | 0 0
  1. Solve [M*S + (T/0.755) - (T/0.755)*cosh[M/(2T/0.755)] == 0, T]
  2.  
  3. eqn = M S+(T/0.755)-(T/0.755)*Cosh[M/(2T/0.755)]==0;
  4.  
  5. seqn = Assuming[M>0, eqn /. T->t M //Simplify];
  6. seqn //TeXForm
  7.  
  8. ode = D[seqn /. t->t[S], S];
  9. ode //TeXForm
  10.  
  11. t1 = t /. First @ Solve[seqn && 0<t<10 /. S->1, t]
  12.  
  13. sol = NDSolveValue[{ode, t[1] == t1}, t, {S, 0, 10}];
  14.  
  15. Plot[sol[t], {t, 0, 10}]
  16.  
  17. Tsol[s_, m_] := m sol[s]
  18.  
  19. Plot3D[eqn[[1]] /. T->Tsol[S, M], {S, 0, 10}, {M, 0, 10}]
Add Comment
Please, Sign In to add comment