Dear sage user,
Could you help me to find z = z(t)?, if:
((2.867e28) / (sqrt (r^3 - (2.1e13)^3 )) ) == (( (2.121e8)*(sqrt(2*(z^2) -1))) / ( z - sqrt (z^2 - 1)))
with,
r == ((1.97774e13)*(sqrt(sqrt(t)))) + ((3.55214e12)*(sqrt(t)))
1 < t < 6000
z(t= 1) = 100
r (t= 1) = 2.332954e13
thanks!!
I couldn't solve it!!!
A differential equation would involve the derivative z'(t) which you don't have in your equation(s).
If you substitute the value of r into the first equation you would have a single non-linear equation relating z and t. You may be able to solve this equation using the solve command. You don't have any unknown parameters in the expression for r or the equation involving z so you might not be able to solve the equation and have the values at t=1 that you specify be true.
You might try something like:
r,z,t = var('r,z,t')
equation = ((2.867e28) / (sqrt (r^3 - (2.1e13)^3 )) ) == (( (2.121e8)*(sqrt(2*(z^2) -1))) / ( z - sqrt (z^2 - 1)))
equation = equation.subs(r == ((1.97774e13)*(sqrt(sqrt(t)))) + ((3.55214e12)*(sqrt(t))))
solve(equation, t)
When I tried this I got:
Is (2*z^2-1)*(z-sqrt(z^2-1)) positive, negative, or zero?
meaning that Maxima (the solver underneath Sage) can't solve the equation without more information.
Asked: 2011-08-08 12:48:57 +0100
Seen: 669 times
Last updated: Aug 08 '11
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
This isn't a system of linear equations, nor a system of ordinary differential equations, unless you've made a typo. Try rephrasing your question, or looking at the Sage documentation about solving linear or differential equations.
could you tell me,what kind of equation do you think is it?. I found this equation from physical principles of shock wave to describe the temporal evolution of relativistic blastwaves. z is the lorentz factor ( http://en.wikipedia.org/wiki/Lorentz_factor ) and r is the shock radius. Thank you for your comment!