Introducing LaTeX Math Typesetting

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
766 replies · 754K views
Let the function [itex]f[/itex] be continuos on the closed interval [itex][a, b][/itex], and assume that [itex]f(x) \geq 0[/itex] for all [itex]x[/itex] in [itex][a, b][/itex]. If [itex]S[/itex] is the solid of revolution obtained by revolving about the [itex]x[/itex] axis the region bounded by the curve [itex]y = f(x)[/itex], the [itex]x[/itex] axis, and the lines [itex]x = a[/itex] and [itex]x = b[/itex], and if [itex]V[/itex] is the number of cubic units in the volume of [itex]S[/itex], then

[tex]V = \pi \int^b_{a}[f(x)]^2 dx[/tex]
 
Last edited:
Physics news on Phys.org
Just a test

[tex]a^_x[/tex]
 
I think we all know Taylor expansion:

[tex]\boxed{f\left( x \right) = \sum\limits_{n = 0}^\infty {\frac{{f^{\left( n \right)} \left( a \right)}}{{n!}}\left( {x - a} \right)^n } ,\left| {x - a} \right| < R}[/tex]
 
[tex]F_{1}+F{2}[/tex]

just testing
 
Last edited:
level: 1
[tex] \sum_{j=0}^m a + cj = \frac{(m + 1)(2a + cm)}{2}[/tex]

level: 2
[tex] \sum_{j=0}^m \frac{(j + 1)(2a + cj)}{2} = \frac{(m + 1)(m + 2)(3a + cm)}{6}\\[/tex]

level: 3
[tex] \sum_{j=0}^m \frac{(j + 1)(j + 2)(3a + cj)}{6} = \frac{(m + 1)(m + 2)(m + 3)(4a + cm)}{24}\\[/tex]

level: 4
[tex] \sum_{j=0}^m \frac{(j + 1)(j + 2)(j + 3)(4a + cj)}{24} = \frac{(m + 1)(m + 2)(m + 3)(m + 4)(5a + cm)}{120}\\[/tex]

level: n
[tex] = \frac{(m + n)!(a(n + 1) + cm)}{m!(n + 1)!}[/tex]

As product
[tex] \frac{(m + n)!(a(n + 1) + cm)}{m!(n + 1)!} = (a(n + 1) + cm) \prod_{j=1}^n \frac{m + j}{j + 1}[/tex]
 
Did I mention Mathematica has an Eigenvector[Matrix] command? However, I'm not good at calculating eigenvectors so I really should do a few by hand:

The eigenvector equation is simple:

[tex] \mathbf{M}v=\lambda v[/tex]

So for:

[tex] \lambda_1=-1[/tex]

[tex] \left(<br /> \begin{array}{ccc} 1 & 0 & 3 \\<br /> 0 & -1 & 0 \\<br /> -3 & 0 & 1<br /> \end{array}<br /> \right)<br /> \left(<br /> \begin{array}{c} x \\<br /> y \\<br /> z<br /> \end{array}<br /> \right)=-1<br /> \left(<br /> \begin{array}{c} x \\<br /> y \\<br /> z<br /> \end{array}<br /> \right)[/tex]

So:

[tex] x+3z=-x\\<br /> -y=-y\\<br /> -3x+z=-z[/tex]

The middle one is easy:

[tex] 0y=0[/tex]

That means y can be anything so let y=1.
The other two:

[tex] 2x+3z=0<br /> -3x+2z=0[/tex]

The simple thing here, since we're looking for ANY eigenvector, is to just pick the zero solution and thus:

[tex] v_1=<br /> \left(<br /> \begin{array}{c} 0 \\<br /> 1 \\<br /> 0<br /> \end{array}<br /> \right)[/tex]

For:
[tex] \lambda_2=(1+3i)[/tex]

[tex] \left(<br /> \begin{array}{ccc} 1 & 0 & 3 \\<br /> 0 & -1 & 0 \\<br /> -3 & 0 & 1<br /> \end{array}<br /> \right)<br /> \left(<br /> \begin{array}{c} x \\<br /> y \\<br /> z<br /> \end{array}<br /> \right)=(1+3i)<br /> \left(<br /> \begin{array}{c} x \\<br /> y \\<br /> z<br /> \end{array}<br /> \right)[/tex]

So that's:

[tex] x+3z=(1+3i)x<br /> -y=(1+3i)y<br /> -3x+z=(1+3i)z[/tex]

For the middle one, the only way ay=y is if y=0. The other two yield:
[tex] 3z=3ix[/tex]

or:

z=ix

so let x=1 and z=i.

Thus:

[tex] v_2=<br /> \left(<br /> \begin{array}{c} 0 \\<br /> 1 \\<br /> 0<br /> \end{array}<br /> \right)[/tex]

Same dif for the other eigenvalue which yields:
[tex] v_3=<br /> \left(<br /> \begin{array}{c} 1 \\<br /> 0 \\<br /> -i<br /> \end{array}<br /> \right)[/tex]

Mathematica returns equivalent eigenvectors.
Thus we are led to the solution in matrix form:

[tex] \mathbf{x}=c_1e^{-t}<br /> \left(<br /> \begin{array}{c} 0 \\<br /> 1 \\<br /> 0<br /> \end{array}<br /> \right)+<br /> c_2e^{(1+3i)t}<br /> \left(<br /> \begin{array}{c} 1 \\<br /> 0 \\<br /> i<br /> \end{array}<br /> \right)+<br /> c_3e^{(1-3i)t}<br /> \left(<br /> \begin{array}{c} 1 \\<br /> 0 \\<br /> -i<br /> \end{array}<br /> \right)[/tex]

You ever work a problem and the answer is just as difficult as the question?
 
Last edited:
After reviewing, I wish to clear up two points in my efforts to solve this equation:

1. There is no need to directly calculate the eigenvector of [itex]\lambda_3[/itex]:

The complex conjugate of an eigenvector for [itex]\lambda_2[/itex] is an eigenvector for [itex]\lambda_3[/tex].<br /> <br /> So above I calculated the eigenvector for [itex](1+3i)[/itex] to be:<br /> <br /> [tex]v_2=\left(\begin{array}{c} 1 \\0 \\i\end{array}\right)[/tex]<br /> <br /> Therefore to calculate the eigenvector for [itex](1-3i)[/itex], conjugate the eigenvector for [itex](1+3i)[/tex]:<br /> <br /> If:<br /> <br /> [tex]v_2=\left(\begin{array}{c} 1+0i \\0+0i \\0+i\end{array}\right)[/tex]<br /> <br /> Then:<br /> <br /> [tex]\overline{v_2}=v_3=\left(\begin{array}{c} 1-0i \\0-i \\0-i\end{array}\right)=\left(\begin{array}{c} 1 \\0 \\-i\end{array}\right)[/tex]<br /> <br /> See how that works?<br /> <br /> Ok that's one.<br /> <br /> 2. The complex eigenvalues both yield the same solution! That is, the solution from (1+3i) is the same solution as that from (1-3i). Go figure. I did. So that's the reason we only need calculate the Real and Complex contribution from ONE member of each complex pair. And also, it's VERY convenient to write the eigenvectors as:<br /> <br /> [tex]\left(\begin{array}{c} 1 \\0 \\-i\end{array}\right)=\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)+i\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)[/tex]<br /> <br /> Alright, so let's compute the Real part and the Complex part for:<br /> <br /> [tex] \begin{align*}<br /> e^{(1+3i)t}\left(\begin{array}{c} 1 \\0 \\i\end{array}\right)<br /> &=e^{(1+3i)t}\left[\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)+i\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)\right] \\<br /> &=e^t\left[(Cos(3t)+iSin(3t))\left\{\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)+<br /> i\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)\right\}\right] \\<br /> &=e^t\left[Cos(3t)\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)+iCos(3t)\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)+iSin(3t)\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)-Sin(3t)\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)\right] \\<br /> &=e^t\left[C_1\left\{Cos(3t)\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)-Sin(3t)\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)\right\}+C_2\left\{Cos(3t)\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)+Sin(3t)\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)\right\}\right]<br /> \end{align}[/tex]<br /> <br /> This with the first solution then yields the general solution:<br /> <br /> [tex] \begin{align*}<br /> \mathbf{X}&=C_1e^{-t}\left(\begin{array}{c} 0 \\1 \\0\end{array}\right) \\<br /> &+e^t\left[C_2\left\{Cos(3t)\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)-Sin(3t)\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)\right\}+C_3\left\{Cos(3t)\left(\begin{array}{c} 0 \\0 \\1\end{array}\right)+Sin(3t)\left(\begin{array}{c} 1 \\0 \\0\end{array}\right)\right\}\right]<br /> \end{align}[/tex]<br /> <br /> That's read as:<br /> <br /> [tex]x(t)=C_2e^tCos(3t)+C_3e^tSin(3t)[/tex]<br /> <br /> [tex]y(t)=C_1e^{-t}[/tex]<br /> <br /> [tex]z(t)=-C_2e^tSin(3t)+C_3e^tCos(3t)[/tex][/itex][/itex]
 
Last edited:
[tex]avg=\frac{n_1+n_2+n_3+n_4+n_5}{n}[/tex][tex]dv_n=average-reading_n[/tex]
[tex]\delta=\sqrt{\frac{dv_1^2+dv_2^2+dv_3^2+dv_4^2+dv_5^2}{n-1}}[/tex]
[tex]\delta_{avg}=\frac{\delta}{\sqrt{n}}[/tex][tex]v_n=\frac{d_n}{t_n}[/tex][tex]\Delta v_n=v_n\sqrt{(\frac{\Delta L}{l})^2+(\frac{\Delta t_n}{t_n})^2}[/tex][tex]\rho_n=m_n(\frac{3}{4\pi})(\frac{d_n}{2})^{-3}[/tex][tex]\Delta\rho_n=6\sqrt{\frac{9\Delta d_n^2m_n^2}{\pi^2d_n^8}+\frac{\Delta m_n^2}{\pi^2d_n^6}}[/tex]

[tex]\eta_n=(\frac{2g}{9v_n})(\frac{d_n}{2})^2(\rho_n-\rho_l)[/tex]
 
Last edited by a moderator:
[tex]\Delta\eta_n=\frac{2}{9}\sqrt{\frac{4\Delta r_n^2g^2r_n^2(\rho_n-\rho_l)^2}{v_n^2}+\frac{\Delta\rho_n^2g^2r_n^4}{v_n^2}+\frac{\Delta g^2r_n^4(\rho_n-\rho_l)^2}{v_n^2}+\frac{\Delta\rho_l^2g^2r_n^4}{v_n^2}+\frac{\Delta v_n^2g^2r_n^4(\rho_n-\rho_l)^2}{v_n^4}}[/tex]
 
[tex]m_{avg}=\frac{19.837+19.839+19.840+19.841+19.840}{5}=19.84 g[/tex]

[tex]v_G=\frac{0.50}{24.25}= 0.0206 m/s[/tex]

[tex]\rho_G=0.01984(\frac{3}{4\pi})(\frac{0.02451}{2})^{-3}= 2573 kg/m^3[/tex]

[tex]\eta_G=(\frac{2g}{9(0.0206)})(\frac{0.02451}{2})^2(2573-1013)= 24.77 kg/ms[/tex]
 
Last edited by a moderator:
[tex]dv_1=24.51-24.70= -0.19 mm[/tex]
[tex]dv_2=24.51-24.40= 0.11 mm[/tex]
[tex]dv_3=24.51-24.44= 0.07 mm[/tex]
[tex]dv_4=24.51-24.32= 0.19 mm[/tex]
[tex]dv_5=24.51-24.68= -0.17mm[/tex]
[tex]\delta_G=\sqrt{\frac{-0.19^2+0.11^2+0.07^2+0.19^2+-0.17^2}{5-1}}= 0.1718 mm[/tex]
[tex]\delta_{avg of G}=\frac{0.1718}{\sqrt{5}}= 0.07683 mm[/tex]

[tex]\Delta v_G=0.0206\sqrt{(\frac{0.003}{0.50})^2+(\frac{\0.1566}{24.25})^2}= 0.00018 m/s[/tex]

[tex]\Delta\rho_G=6\sqrt{\frac{9(0.00007683)^2(0.01984)^2}{\pi^20.02451^8}+\frac{0.00005^2}{\pi^20.02451^6}}= 25.05 kg/m^3[/tex]
 
Last edited by a moderator:
[tex]\Delta\eta_G=\frac{2}{9}\sqrt{\frac{4(0.000038415)^2(9.8)^2(0.012255)^2(2573-1013)^2}{0.0206^2}+\frac{(25.05)^2(9.8)^2(0.012255)^4}{0.0206^2}[/tex]

[tex]\sqrt{adfsadfasf+\frac{(0.005)^2(0.012255)^4(2573-1013)^2}{0.0206^2}+\frac{(5)^2(9.8)^2(0.012255)^4}{0.0206^2}<br /> +\frac{(0.00018)^2(9.8)^2(0.012255)^4(2573-1013)^2}{0.0206^4}}= 0.4854 kg/ms}[/tex]
 
Last edited by a moderator:
First, let's clean up the functional: Really, if we want to minimize that integral, we can just move U across the integral sign right, and let's put the exponential in the numerator:

[tex]\mathbf{T}[y(x)]=\frac{1}{U}\int_{p_1}^{p_2} e^{y/h}\sqrt{1+(y^{'})^2}dx[/tex]


So:

[tex]F(x,y,y^{'})=e^{y/h}\sqrt{1+(y^{'})^2}[/tex]

and therefore:

[tex]\frac{\partial F}{\partial y}=\frac{e^{y/h}\sqrt{1+(y^{'})^2}}{h}[/tex]

and:

[tex]\frac{\partial F}{\partial y^{'}}=\frac{e^{y/h}y^{'}}{\sqrt{1+(y^{'})^2}}[/tex]

and so:

[tex] \begin{align*}<br /> \frac{d}{dx}\left(\frac{\partial F}{\partial y^{'}}\right)&=\frac{d}{dx}\left[\frac{e^{y/h}y^{'}}{\sqrt{1+(y^{'})^2}}\right] \\<br /> <br /> &=\frac{d}{dx}\left[e^{y/h}y^{'} \cdot \frac{1}{\sqrt{1+(y^{'})^2}}\right] \\<br /> <br /> &=\left[e^{y/h}y^{'}\cdot\frac{-1/2}{(1+(y^{'})^2)^{3/2}}\cdot 2 y^{'}y^{''} \\<br /> <br /> &+\frac{1}{\sqrt{1+(y^{'})^2}}\left(e^{y/h}y^{''}+y^{'}\frac{1}{h}y^{'}e^{y/h}\right) \\<br /> <br /> &=\frac{e^{y/h}y^{''}}{\sqrt{1+(y^{'})^2}}-\frac{e^{y/h}(y^{'})^2y^{''}}{(1+(y^{'})^2)^{3/2}}+<br /> \frac{e^{y/h}(y^{'})^2}{h\sqrt{1+(y^{'})^2}}<br /> <br /> <br /> <br /> <br /> <br /> <br /> <br /> \end{align}[/tex]

So, once we obtain the partials, then we substitute them into the Euler equation and equate the expression to zero. Now, can you please substitute these expressions into:

[tex]\frac{\partial F}{\partial y}-\frac{d}{dx}\left(\frac{\partial F}{\partial y^{'}}\right)=0[/tex]

and post the results?

Also, with regards to h=1: I just worked the problem with that value and obtained the results you indicated. Perhaps it works for any value of h. Not sure.
 
Last edited:
[tex]\Delta\eta_G=\frac{2}{9}\sqrt{\frac{4(0.000038415) ^2(9.8)^2(0.012255)^2(2573-1013)^2}{0.0206^2}+\frac{(25.05)^2(9.8)^2(0.012255 )^4}{0.0206^2}+\frac{(0.005)^2(0.012255)^4(2573-1013)^2}{0.0206^2}+\frac{(5)^2(9.8)^2(0.012255)^4} {0.0206^2}+\frac{(0.00018)^2(9.8)^2(0.012255)^4(2573-1013)^2}{0.0206^4}}= 0.4854 kg/ms}[/tex]
 
[tex]\overrightarrow{F_f}=-6\pi\eta(\frac{d}{2})\overrightarrow{v}[/tex]

[tex]\overrightarrow{mg}=\frac{4}{3}\pi(\frac{d}{2})^3\rho_s\overrightarrow g[/tex]

[tex]\overrightarrow{F_b}=-\frac{4}{3}\pi(\frac{d}{2})^3\rho_l\overrightarrow g[/tex]

[tex]\overrightarrow{F_b}+\overrightarrow{F_f}+\overrightarrow{mg} = 0[/tex]
 
Last edited by a moderator:
test

[tex] \int_{0}^{1} x dx = \left[ \frac{1}{2}x^2 \right]_{0}^{1} = \frac{1}{2}[/tex]
Let [tex]\theta[/tex] be a (p+1)-form defined on the ranges of all the cubes of a p-chain, where p>0. Then [tex]\int_{C} d \theta = \int_{\partial C} \theta[/tex]
[tex] Let\; \theta\; be\; a (p+1)-form\; defined\; on\; the\; ranges\; of\; all\; the\; cubes\; of\; a\; p-chain,\; where\; p>0.\; Then \int_{C} d \theta = \int_{\partial C} \theta.[/tex]
 
Last edited:
[tex]e=\frac{f}{a}[/tex]
[tex]e=\sqrt{1-\frac{b^2}{a^2}}[/tex]
 
[tex]F_b=m_ba_b[/tex]
[tex]F_e=m_ea_e[/tex]
[tex]F_b=F_e[/tex]
[tex]m_ba_b=m_ea_e[/tex]
[tex]\frac{m_b}{m_e}=\frac{a_e}{a_b}[/tex]
 
s

[tex]\nabla \cross E=0[/tex]
[tex]\nabla \cdot E =\frac {\rho}{\epsilon_0}[/tex]

abc
 
Last edited:
Taken [tex]\sum F = m_1 \cdot a[/tex]
We have [tex]\sum F_x = m_1 \cdot a_x[/tex] and [tex]\sum F_y = m_1 \cdot a_y[/tex]
With [tex]a_y = 0[/tex]
We We have [tex]\sum F_x = m_1 \cdot a_x[/tex] and [tex]\sum F_y = 0[/tex]
So [tex]T = m_1 \cdot a_x[/tex]
Taken [tex]\sum F = m_2 \cdot a[/tex]
We have [tex]\sum F_y = m_2 \cdot a_y[/tex]
Which is [tex]m_2 \cdot g - T = m_2 \cdot a[/tex]
Inserting [tex]T = m_1 \cdot a_x[/tex] into [tex]m_2 \cdot g - T = m_2 \cdot a[/tex] and solving for [tex]a[/tex],
We get [tex]a = \frac {m_2 \cdot g} {m_1 + m_2}[/tex]
Inserting this into [tex]T = m_1 \cdot a[/tex],
we get [tex]T = \frac {m_1 \cdot m_2 \cdot g} {m_1 = m_2}[/tex]Now, taken [tex]v^2 = v_0^2 + 2a\deltax[/tex]
 
Last edited:
Taken [tex]\sum F = m_1 \cdot a[/tex]
We have [tex]\sum F_x = m_1 \cdot a_x[/tex] and [tex]\sum F_y = m_1 \cdot a_y[/tex]
With [tex]a_y = 0[/tex]
We We have [tex]\sum F_x = m_1 \cdot a_x[/tex] and [tex]\sum F_y = 0[/tex]
So [tex]T = m_1 \cdot a_x[/tex]
Taken [tex]\sum F = m_2 \cdot a[/tex]
We have [tex]\sum F_y = m_2 \cdot a_y[/tex]
Which is [tex]m_2 \cdot g - T = m_2 \cdot a[/tex]
Inserting [tex]T = m_1 \cdot a_x[/tex] into [tex]m_2 \cdot g - T = m_2 \cdot a[/tex] and solving for [tex]a[/tex],
We get [tex]a = \frac {m_2 \cdot g} {m_1 + m_2}[/tex]
Inserting this into [tex]T = m_1 \cdot a[/tex],
we get [tex]T = \frac {m_1 \cdot m_2 \cdot g} {m_1 = m_2}[/tex]


Now, taken [tex]v^2 = v_0^2 + 2a\delta x[/tex]
 
Taken [tex]\sum F = m_1 \cdot a[/tex]
We have [tex]\sum F_x = m_1 \cdot a_x[/tex] and [tex]\sum F_y = m_1 \cdot a_y[/tex]
With [tex]a_y = 0[/tex]
We We have [tex]\sum F_x = m_1 \cdot a_x[/tex] and [tex]\sum F_y = 0[/tex]
So [tex]T = m_1 \cdot a_x[/tex]
Taken [tex]\sum F = m_2 \cdot a[/tex]
We have [tex]\sum F_y = m_2 \cdot a_y[/tex]
Which is [tex]m_2 \cdot g - T = m_2 \cdot a[/tex]
Inserting [tex]T = m_1 \cdot a_x[/tex] into [tex]m_2 \cdot g - T = m_2 \cdot a[/tex] and solving for [tex]a[/tex],
We get [tex]a = \frac {m_2 \cdot g} {m_1 + m_2}[/tex]
Inserting this into [tex]T = m_1 \cdot a[/tex],
we get [tex]T = \frac {m_1 \cdot m_2 \cdot g} {m_1 + m_2}[/tex]
Now, taken [tex]v^2 = v_0^2 + 2a\Delta x[/tex]
Solving for [tex]a[/tex] we get [tex]a = \frac {v^2 - V_0^2} {2\Delta x}[/tex]
Taking that [tex]\Delta x = .5[/tex] and [tex]v_0 = 0[/tex]
We have [tex]a = \frac {v^2} {meters*}[/tex]
Now [tex]a_exp[/tex] stands for experimental value of [tex]a[/tex]
and [tex]a_t [/tex[ stands for theoretical value of [tex]a[/tex]<br /> so [tex]a_exp = \frac {v^2} {meters*}[/tex]<br /> and [tex]a_t = \frac {m_2 \cdot g} {m_1 + m_2}[/tex]<br /> * meters is added so units sovle correctly.[/tex]
 
Taken [tex]\sum F = m_1 \cdot a[/tex]
We have [tex]\sum F_x = m_1 \cdot a_x[/tex] and [tex]\sum F_y = m_1 \cdot a_y[/tex]
With [tex]a_y = 0[/tex]
We We have [tex]\sum F_x = m_1 \cdot a_x[/tex] and [tex]\sum F_y = 0[/tex]
So [tex]T = m_1 \cdot a_x[/tex]
Taken [tex]\sum F = m_2 \cdot a[/tex]
We have [tex]\sum F_y = m_2 \cdot a_y[/tex]
Which is [tex]m_2 \cdot g - T = m_2 \cdot a[/tex]
Inserting [tex]T = m_1 \cdot a_x[/tex] into [tex]m_2 \cdot g - T = m_2 \cdot a[/tex] and solving for [tex]a[/tex],
We get [tex]a = \frac {m_2 \cdot g} {m_1 + m_2}[/tex]
Inserting this into [tex]T = m_1 \cdot a[/tex],
we get [tex]T = \frac {m_1 \cdot m_2 \cdot g} {m_1 + m_2}[/tex]
Now, taken [tex]v^2 = v_0^2 + 2a\Delta x[/tex]
Solving for [tex]a[/tex] we get [tex]a = \frac {v^2 - V_0^2} {2\Delta x}[/tex]
Taking that [tex]\Delta x = .5[/tex] and [tex]v_0 = 0[/tex]
We have [tex]a = \frac {v^2} {meters*}[/tex]
Now [tex]a_exp[/tex] stands for experimental value of [tex]a[/tex]
and [tex]a_t[/tex] stands for theoretical value of [tex]a[/tex]
so [tex]a_exp = \frac {v^2} {meters*}[/tex]
and [tex]a_t = \frac {m_2 \cdot g} {m_1 + m_2}[/tex]
* meters is added so units sovle correctly.
 
I've achieved oneness with the integral Tide . . . Here it is with your scalling factor in case other people are following this:

Above, after letting:

[tex]\sigma^2=s[/tex]

and completing the square, we obtain:

[tex]2e^{-\pi/t}\int e^{t(\sigma-\sqrt{\pi}/t)^2}\sigma d\sigma[/tex]

Now, in order to remove the t in the exponent, we let:

[tex]v=\sqrt{t}\sigma[/tex]

so that:

[tex]dv=\sqrt{t}d\sigma,\quad \sigma=\frac{v}{\sqrt{t}},\quad d\sigma=\frac{dv}{\sqrt{t}}[/tex]

Substituting this scalling factor into the exponent:

[tex] \begin{align*}<br /> t\left[\frac{v^2}{t}-\frac{2v\sqrt{\pi}}{t\sqrt{t}}+\frac{\pi}{t^2}\right]&=<br /> v^2-2v\sqrt{\pi/t}+\frac{\pi}{t} \\<br /> &=(v-\sqrt{\pi/t})^2<br /> \end{align}[/tex]

substituting this into the integral:

[tex]2e^{-\pi/t}\int e^{(v-\sqrt{\pi/t})^2}\left(\frac{v}{\sqrt{t}}\right)\frac{dv}{\sqrt{t}}[/tex]

Simplifying:

[tex]2\frac{e^{-\pi/t}}{t}\int e^{(v-\sqrt{\pi/t})^2}dv;\quad v=\sqrt{t}\sigma;\quad \sigma^2=s[/tex]
 
Last edited: