Undergrad Micromass' big October challenge

Click For Summary
The October challenge invites participants to solve a series of mathematical problems, with a focus on providing full proofs for their solutions. Participants are encouraged to submit new challenge ideas and adhere to specific rules, such as not using direct searches for problem solutions. Advanced challenges include topics like primitive recursive functions, generalized limits, and dynamic systems involving motion. Several problems have already been solved, showcasing a range of mathematical concepts and techniques. The challenge aims to foster engagement and collaboration among math enthusiasts.
  • #61
micromass said:
That is not what I wrote in the OP.
Sorry. I think what happened here is I read the initial problem and then began working on it, but I didn't have the coordinates in precise agreement because I didn't go back and re-read the question thoroughly. In my first couple of posts, I didn't think I was going to come close to solving it... editing... I think my slightly modified coordinates might even simply some of the vector algebra because it always has an x=0 for one of the ships...## \\ ## The substitution ## x=x_o-x' ## will give ## y ## vs. ## x' ## where the ## x' ## is the correct coordinate given in the OP. Putting this addition in the equation of post 57 hopefully gives the correct answer. My answer of post 57 for D is in agreement with @Chestermiller of post 47 for ## v_1=3 v_o ##.
 
Last edited:
Mathematics news on Phys.org
  • #62
mfb said:
You cannot link to your private mails. Even if you can see the image (I don't know), no one else can. You can use the "Upload" button to upload images here.

But it's size is greater that 3 MB , so what do I do ??
 
  • #63
You can also upload it to various websites that host images. A better solution would be to reduce its size, if that doesn't ruin the quality.
 
  • #64
Charles Link said:
I made too many mistakes in post #48. Here I will start with my equation from post #43 ## x/x_o=|sec(\phi)+tan(\phi)|^{-v_o/v_1} ## and proceed with ## dy/dx=tan(\phi) ## and ## sec(\phi)=-\sqrt{(dy/dx)^2+1} ##. A little algebra gives ## dy/dx=(1/2)[(x/x_o)^{-v_1/v_o}-(x/x_o)^{v_1/v_o}] ##. This integrates to give ## \\ ## ## y=(x_o/2)[-(v_o/(v_1-v_o))(x/x_o)^{-(v_1-v_o)/v_o}-(v_o/(v_1+v_o))(x/x_o)^{(v_1+v_o)/v_o)} ]+C ## ## \\ ## where ## C=x_o v_o v_1/(v_1^2-v_o)^2 ##. The constant of integration occurs because at ## x=x_o ## , ## y=0 ##. To find where this curve intercepts the good ship, we set x=0. This gives ## y_1=x_o v_o v_1/(v_1^2-v_o^2) ##. Distance D traveled by the pirate ship is ## D=(v_1/v_o)y_1=x_o v_1^2/(v_1^2-v_o^2) ##.
My solution to this problem for x vs y, expressed parametrically in terms of ##\theta## is as follows:
$$x=x_0\left[1-\frac{1}{(\sec{\theta}+\tan{\theta})^{V/v}}\right]$$
$$y=x_0\left[\int_0^{\theta}{\frac{\sec^2{\theta '}d\theta '}{(\sec{\theta '}+\tan{\theta '})^{V/v}}}-\frac{\tan {\theta}}{(\sec{\theta}+\tan{\theta})^{V/v}}\right]$$where ##\theta '## is a dummy variable of integration.

Charles, your analytic solution to this problem should match mine, and should thus somehow provide the result of correctly integrating of my "mystery integral" in the equation for y. Could you please see if you can back out the integral evaluation? Thanks.

Chet
 
  • #65
Charles Link said:
I made too many mistakes in post #48. Here I will start with my equation from post #43 ## \\ ## ## x/x_o=|sec(\phi)+tan(\phi)|^{-v_o/v_1} ## and proceed with ## dy/dx=tan(\phi) ## and ## sec(\phi)=-\sqrt{(dy/dx)^2+1} ##. A little algebra gives ## dy/dx=(1/2)[(x/x_o)^{-v_1/v_o}-(x/x_o)^{v_1/v_o}] ##. This integrates to give ## \\ ## ## y=(x_o/2)[-(v_o/(v_1-v_o))(x/x_o)^{-(v_1-v_o)/v_o}-(v_o/(v_1+v_o))(x/x_o)^{(v_1+v_o)/v_o)} ]+C ## ## \\ ## where ## C=x_o v_o v_1/(v_1^2-v_o)^2 ##. The constant of integration occurs because at ## x=x_o ## , ## y=0 ##. To find where this curve intercepts the good ship, we set x=0. This gives ## y_1=x_o v_o v_1/(v_1^2-v_o^2) ##. Distance D traveled by the pirate ship is ## D=(v_1/v_o)y_1=x_o v_1^2/(v_1^2-v_o^2) ##.
@micromass With a little effort I found a couple of algebraic errors in my previous posts. Hopefully this one is now correct. ## \\ ## Adding on to this and making a couple of corrections, thanks to @Chestermiller post #47, I see I have an error in post 43, and my equation should read (instead of ## v_o/v_1 ##) ## \\ ## ## x=x_o|sec(\phi)+tan(\phi)|^{-v_1/v_o} ##. ## \\ ## My equation above then becomes ## \\ ## ## dy/dx=-(1/2)[(x/x_o)^{-v_o/v_1}-(x/x_o)^{v_o/v_1}] ## ## \\ ## where the addition of the minus sign (missing in the above quote) in front of the (1/2) comes as a result of the absolute value sign putting a minus sign on the dy/dx which for my choice of coordinates is always negative. Solving for y (integrating), the result is ## \\ ## ## y= (1/2) x_o [(v_1/(v_1+v_o))(x/x_o)^{(v_1+v_o)/v_o}-(v_1/(v_1-v_o))(x/x_o)^{(v_1-v_o)/v_1}]+C ## ## \\ ## where ## C= x_o v_o v_1/(v_1^2-v_o^2) ##. ## \\ ## The substitution ## x=x_o-x' ## is then needed to get in the form of the OP because I did a coordinate transformation where I had the good ship sailing up the y-axis and the pirate ship beginning at ## (x_o,0) ##... Making these couple of corrections did not affect the distance ## D ## that the pirate ship travels which is ## D=x_o v_1^2/(v_1^2-v_o^2) ##. Note: ## v_1=V ## and ## v_o=v ## of the OP.
 
Last edited:
  • #66
Charles Link said:
Adding on to this and making a couple of corrections, thanks to @Chestermiller post #47, I see I have an error in post 43, and my equation should read (instead of ## v_o/v_1 ##) ## \\ ## ## x=x_o|sec(\phi)+tan(\phi)|^{-v_1/v_o} ##. ## \\ ## My equation above then becomes ## \\ ## ## dy/dx=-(1/2)[(x/x_o)^{-v_o/v_1}-(x/x_o)^{v_o/v_1}] ## ## \\ ## where the addition of the minus sign (missing in the above quote) in front of the (1/2) comes as a result of the absolute value sign putting a minus sign on the dy/dx which for my choice of coordinates is always negative. Solving for y (integrating), the result is ## \\ ## ## y= (1/2) x_o [(v_1/(v_1+v_o))(x/x_o)^{(v_1+v_o)/v_o}-(v_1/(v_1-v_o))(x/x_o)^{(v_1-v_o)/v_1}+C ## ## \\ ## where ## C= x_o v_o v_1/(v_1^2-v_o^2) ##. ## \\ ## The substitution ## x=x_o-x' ## is then needed to get in the form of the OP because I did a coordinate transformation where I had the good ship sailing up the y-axis and the pirate ship beginning at ## (x_o,0) ##... Making these couple of corrections did not affect the distance ## D ## that the pirate ship travels which is ## D=x_o v_1^2/(v_1^2-v_o^2) ##. Note: ## v_1=V ## and ## v_o=v ## of the OP.

It still doesn't seem to be correct. If I graph the curve you posted, then I notice that if ##v1>v0##, then the ##y## becomes negative, meaning that the the pirate ships goes in the negative direction first. This seems wrong.
 
  • #67
micromass said:
It still doesn't seem to be correct. If I graph the curve you posted, then I notice that if ##v1>v0##, then the ##y## becomes negative, meaning that the the pirate ships goes in the negative direction first. This seems wrong.
I will need to take another look at it but my x goes from ## x_o ## to 0 as the ## x' ## (which is the coordinate I should be using) goes from 0 to ## x_o ##. Did you include the constant of integration C in the graph ? I also missed a " ]" which I will now include.editing... I found the typo: THe denominator of the first exponent should be ## v_1 ##. Hopefully that makes it work.
 
Last edited:
  • #68
@micromass I found a typo: The denominator in the first exponent should be ## v_1 ##. It correctly reads ## \\ ## ## y= (1/2) x_o [(v_1/(v_1+v_o))(x/x_o)^{(v_1+v_o)/v_1}-(v_1/(v_1-v_o))(x/x_o)^{(v_1-v_o)/v_1}]+C ## ## \\ ## where ## C= x_o v_o v_1/(v_1^2-v_o^2) ##. ## \\ ## Hopefully this is now finally correct.
 
  • Like
Likes micromass
  • #69
Charles Link said:
@micromass I found a typo: The denominator in the first exponent should be ## v_1 ##. It correctly reads ## \\ ## ## y= (1/2) x_o [(v_1/(v_1+v_o))(x/x_o)^{(v_1+v_o)/v_1}-(v_1/(v_1-v_o))(x/x_o)^{(v_1-v_o)/v_1}]+C ## ## \\ ## where ## C= x_o v_o v_1/(v_1^2-v_o^2) ##. ## \\ ## Hopefully this is now finally correct.

That seems ok, but the solution is invalid if ##v_1=v_0##.
 
  • #70
micromass said:
That seems ok, but the solution is invalid if ##v_1=v_0##.
The ## v_1=v_o ## is a special case that has the result ## \\ ## ## y=(1/2)x^2/x_0-(1/2)x_o ln(x/x_o)-(1/2) x_o ## where again the substitution ## x=x_o-x' ## is necessary to get it into the form of the original post.
 
  • Like
Likes micromass
  • #71
Awesome, it appears (a) and (c) are solved then.
 
  • Like
Likes Charles Link
  • #72
And I saw that (b) was solved too!
 
  • Like
Likes Charles Link
  • #73
Ok, I'll try advanced problem number 7. First, we note that ##X=b\tan{\theta}##. To get the expected value of ##X##, we use the law of the unconcscious statistician:
$$E[X(\theta)] = \int^{\frac{\pi}{2}}_{-\frac{\pi}{2}} X(\theta) f(\theta)d\theta$$
where ##f(\theta)## is the probability density function of ##\theta##. We know that ##\theta## is uniformly distributed over ##(-\pi/2,\pi/2)##, so ##f(\theta) = 1/\pi##. Thus the expected value of ##X## is
$$E[X(\theta)] =\frac{b}{\pi} \int^{\frac{\pi}{2}}_{-\frac{\pi}{2}} \tan{\theta} d\theta$$
Since ##\tan \theta## is an odd function and the integration is symmetric about ##x=0##, the expected value is ##E[X] = 0##.

For the variance, we have:
$$\sigma^2 = E([X^2])-(E[X])^2$$
so we need to evaluate ##E([X^2])##. We use the unconscious statistician again, which gives us the integral:
$$E[X^2] = \frac{b^2}{\pi} \int^{\frac{\pi}{2}}_{-\frac{\pi}{2}} \tan^2{\theta} d\theta$$
This integral diverges to infinity as we take the bounds of integration out to ##\pm \pi/2##. This would imply that the variance is infinite ##(\sigma^2 = \infty)##. I'm not sure if this is right, but my intuition says this makes sense, since there's a finite probability of ##|X|## being arbitrarily large. But I don't know that for certain, and I'm not sure how to make it more mathy.
 
  • #74
TeethWhitener said:
Ok, I'll try advanced problem number 7. First, we note that ##X=b\tan{\theta}##. To get the expected value of ##X##, we use the law of the unconcscious statistician:
$$E[X(\theta)] = \int^{\frac{\pi}{2}}_{-\frac{\pi}{2}} X(\theta) f(\theta)d\theta$$
where ##f(\theta)## is the probability density function of ##\theta##. We know that ##\theta## is uniformly distributed over ##(-\pi/2,\pi/2)##, so ##f(\theta) = 1/\pi##. Thus the expected value of ##X## is
$$E[X(\theta)] =\frac{b}{\pi} \int^{\frac{\pi}{2}}_{-\frac{\pi}{2}} \tan{\theta} d\theta$$
Since ##\tan \theta## is an odd function and the integration is symmetric about ##x=0##, the expected value is ##E[X] = 0##.

For the variance, we have:
$$\sigma^2 = E([X^2])-(E[X])^2$$
so we need to evaluate ##E([X^2])##. We use the unconscious statistician again, which gives us the integral:
$$E[X^2] = \frac{b^2}{\pi} \int^{\frac{\pi}{2}}_{-\frac{\pi}{2}} \tan^2{\theta} d\theta$$
This integral diverges to infinity as we take the bounds of integration out to ##\pm \pi/2##. This would imply that the variance is infinite ##(\sigma^2 = \infty)##. I'm not sure if this is right, but my intuition says this makes sense, since there's a finite probability of ##|X|## being arbitrarily large. But I don't know that for certain, and I'm not sure how to make it more mathy.

Your expected value is wrong.
 
  • #75
micromass said:
Your expected value is wrong.
Did I cut too many corners? Upon closer inspection, I get
$$\int^{\frac{\pi}{2}}_{-\frac{\pi}{2}} \tan{\theta} d \theta = \infty - \infty$$
so it's an indeterminate form. I have no idea how to handle this if not by the symmetry of the function.
 
  • #76
micromass said:
Your expected value is wrong.
Really? As far as I can see, the symmetry of the problem implies that the expectation is ##0##. If we take two infinitesimal intervals ##[-x-dx, x]## and ##[x,x+dx]## on the wall, corresponding to the angle intervals ##[-\theta-d\theta, \theta]## and ##[\theta,\theta+d\theta]##, respectively, both having probability ##d\theta/\pi##, their contributions to the expectation are ##-xd\theta/\pi## and ##xd\theta/\pi##, respectively, so they cancel each other out. Every infinitesimal interval has such a "mirror" interval on the other side of ##a##, so the expectation must be ##0##.

How could it be in any other way?
 
  • #77
Erland said:
Really? As far as I can see, the symmetry of the problem implies that the expectation is ##0##. If we take two infinitesimal intervals ##[-x-dx, x]## and ##[x,x+dx]## on the wall, corresponding to the angle intervals ##[-\theta-d\theta, \theta]## and ##[\theta,\theta+d\theta]##, respectively, both having probability ##d\theta/\pi##, their contributions to the expectation are ##-xd\theta/\pi## and ##xd\theta/\pi##, respectively, so they cancel each other out. Every infinitesimal interval has such a "mirror" interval on the other side of ##a##, so the expectation must be ##0##.

How could it be in any other way?
Then by that same reasoning, would you say that ##\int_{-\infty}^{+\infty} xdx = 0## too?
 
  • #78
micromass said:
Then by that same reasoning, would you say that ##\int_{-\infty}^{+\infty} xdx = 0## too?
So is it just an indeterminate form, like in my post #75 above?
 
  • #79
TeethWhitener said:
So is it just an indeterminate form, like in my post #75 above?

Yes, the expectation value doesn't exist.

You can do the test here: http://www.math.uah.edu/stat/apps/CauchyExperiment.html Run it for 1000 or 10000 turns. You'll see it getting close to 0, but then suddenly it'll jump away from 0 again.
 
  • #80
micromass said:
Yes, the expectation value doesn't exist.
Yes, Ok, I agree. The integral ##\int_{-\pi/2}^{\pi/2}\tan\theta d\theta## diverges. Sorry...

But then, the variance is undefined also, since it is defined in terms of the expectation: ##V(X)=E((X-E(X))^2)##.

So the entire problem was a poser... :wink:
 
  • Like
Likes TeethWhitener
  • #81
Erland said:
So the entire problem was a poser... :wink:

Haha, yes indeed!
 
  • #82
Erland said:
But then, the variance is undefined also
Good point. I was assuming that ##E[X]=0##. Call it the "renormalized" answer :smile:
 
  • #83
Ok, I'll try the unsolved advanced number 2. Last month, I was somewhat unadvisedly trying to show that ##dM_p/dp## was positive everywhere, but the hint this month helps out a whole bunch. We want to know the nature of the relationship:
$$\left(\sum_i {\frac{x_i^p}{n}}\right)^\frac{1}{p} \sim \left(\sum_i {\frac{x_i^q}{n}}\right)^\frac{1}{q}$$
For the moment, let us consider both terms raised to the power of ##p##. We have:
$$\sum_i {\frac{x_i^p}{n}}$$
and
$$\left(\sum_i {\frac{x_i^q}{n}}\right)^\frac{p}{q}$$
We notice that
$$\sum_i {\frac{\left(x_i^q\right)^{\frac{p}{q}}}{n}} = \sum_i {\frac{x_i^p}{n}}$$
which allows us to apply Jensen's inequality:
$$\left(\sum_i {\frac{x_i^q}{n}}\right)^\frac{p}{q} \leq \sum_i {\frac{\left(x_i^q\right)^{\frac{p}{q}}}{n}} = \sum_i {\frac{x_i^p}{n}}$$
for ##x^{p/q}## convex (2nd derivative ##> 0##) and
$$\left(\sum_i {\frac{x_i^q}{n}}\right)^\frac{p}{q} \geq \sum_i {\frac{x_i^p}{n}}$$
for ##x^{p/q}## concave (2nd derivative ##< 0##). For positive ##x##, ##x^{p/q}## is convex when ##p/q > 1## and ##p/q < 0## and concave otherwise. (This is true because the second derivative of ##x^a## is ##a(a-1)x^{a-2}## and the function ##a(a-1)## is a parabola which is only negative between 0 and 1).

We want to know what happens for ##p>q##. This condition breaks the problem into 3 cases: ##p## and ##q## both positive (convex), ##p## positive and ##q## negative (convex), and ##p## and ##q## both negative (concave). For the convex case, we have
$$(M_p)^p \geq (M_q)^p$$
Since ##p## is positive whenever ##x^{p/q}## is convex, this result implies that ##M_p \geq M_q##. For the concave case, we have
$$(M_p)^p \leq (M_q)^p$$
Since ##p## is negative whenever ##x^{p/q}## is concave, this result also implies that ##M_p \geq M_q##, which completes the proof.
 
  • Like
Likes micromass
  • #84
Shreyas Samudra said:
https://mail.google.com/mail/u/0/?u...68608010&rm=15794b3f1927dbe0&zw&sz=w1366-h662

https://mail.google.com/mail/u/0/?ui=2&ik=8afd5aa52a&view=fimg&th=15794b3f1927dbe0&attid=0.1&disp=inline&realattid=1547350341991792640-local0&safe=1&attbid=ANGjdJ9AL_gRE5sXmUYaevAhKSEVWeCDGMrM2iIlWdWoI_HQXuUtcdRj5gRNTsNVhx1fD7uUlFXVMOekFVnhgzPXCNVO8-rWTr1U5G2r1pmTcTPjii6g4adAZZBM6Ic&ats=1475668608010&rm=15794b3f1927dbe0&zw&sz=w1366-h662
Is it right ??
20161002_235107_1.jpg
free adult image hosting
 
  • #86
Attempting Advanced Problem 6:

##X## and ##Y## are independent stochastic variables which both have uniform distrubution on ##[0,1]##. This means that for every Lebesgue mesaurable set ##A\subseteq \mathbb R^2##, ##P((X,Y)\in A)=m(A\cap[0,1]^2)##, where ##m## is the Lebesgue measure on ##\mathbb R^2##. Also, for every Lebesgue measurable set ##B\subseteq \mathbb R##, ##P(X^Y\in B)=m(\{(x,y)\in[0,1]^2\,|\,x^y\in B\})##.

The distribution of ##X^Y## is determined by its cumulative distribution function ##F:\mathbb R \to[0,1]##, given by ##F(t)=P(X^Y\le t)=m(\{(x,y)\in[0,1]^2\,|\,x^y\le t\})##, for all ##t\in\mathbb R##. To find the distribution, it should suffice to find ##F##.

The function ##f(x,y)=x^y##, defined on ##[0,1]^2## (for definiteness, we define ##f(0,0)=0^0=1##, but this does not really matter) is increasing in ##x## and decreasing in ##y##. ##f## is continuous on all ##[0,1]^2## except at ##(0,0)##, so it is Lebesgue measurable. We have ##f(x,0)=1## and ##f(x,1)=x##, for all ##x\in[0,1]##, and ##f(0,y)=0## and ##f(1,y)=1## for all ##y\in\,]0,1]##, so the range of ##f## is ##[0,1]##.

Hence, ##F(t)=0## for ##t<0## and ##F(t)=1## for ##t>1##.
For ##t\in[0,1]## and ##(x,y)\in[0,1]\times\,]0,1]##, we have ##x^y\le t\Leftrightarrow x\le t^{\frac1y}(\le 1)##, since ##\frac1y\ge 1##. For each ##\epsilon\in\,]0,1[\,##, we then have ##m(\{(x,y)\in[0,1]\times[\epsilon,1]\,|\,x^y\le t\})=\int_\epsilon^1 t^{\frac1y}dy##. This tends to both ##F(t)## (since ##m(\{(x,0)\,|\,0\le x\le 1\})=0)## and ##\int_0^1t^{\frac1y}## (since this improper integral converges, because the integrand is continuous and bounded on ##]0,1]\,##) as ##\epsilon\to 0##. Hence ##F(t)=\int_0^1 t^{\frac1y} dy##.
I don't know of any closed expression for this integral, and I doubt that such an expression exists or is known.

So my answer is that the desired distribution is given by the cumulative distribution function:

F(t)=\begin{cases} 0 &amp; t&lt;0,\\ \int_0^1 t^{\frac1y}dy,&amp; 0\le t \le 1, \\1,&amp;t&gt;1.\end{cases}
 
  • Like
Likes micromass
  • #87
I just posted one answer, and only after replying I noticed that it said this:

micromass said:
CHALLENGES FOR HIGH SCHOOL AND FIRST YEAR UNIVERSITY:

before that problem. Whoops... Did that mean that I should not have posted that answer?
 
  • #88
jostpuur said:
Whoops... Did that mean that I should not have posted that answer?
Yes, and we've removed it.
 
  • #89
Chestermiller said:
My solution to this problem for x vs y, expressed parametrically in terms of ##\theta## is as follows:
$$x=x_0\left[1-\frac{1}{(\sec{\theta}+\tan{\theta})^{V/v}}\right]$$
$$y=x_0\left[\int_0^{\theta}{\frac{\sec^2{\theta '}d\theta '}{(\sec{\theta '}+\tan{\theta '})^{V/v}}}-\frac{\tan {\theta}}{(\sec{\theta}+\tan{\theta})^{V/v}}\right]$$where ##\theta '## is a dummy variable of integration.

Charles, your analytic solution to this problem should match mine, and should thus somehow provide the result of correctly integrating of my "mystery integral" in the equation for y. Could you please see if you can back out the integral evaluation? Thanks.

Chet
I finally got the integral in the parameteric equation for y. @fresh_42 submitted the integrand to Wolfram Alpha, and it provided the closed-form result. Thanks you so much @fresh_42. Here is the desired result:
$$y=\frac{x_0(V/v)}{\left(\frac{V}{v}\right)^2-1}\left[1-\frac{(\sec{\theta}+(V/v)\tan{\theta})}{(\sec{\theta}+\tan{\theta})^{V/v}}\right]$$
I hope I did the "arithmetic" correctly and that this result agrees with Charles Links' results.
 
Last edited:
  • Like
Likes Charles Link
  • #90
problem 4 highshool:
x coordinate = ##-1+\frac{1}{2}-\frac{1}{4}... = -1 + \frac{1}{2}(1-\frac{1}{2}+\frac{1}{3}...) = -1 + \frac{log(2)}{2}##
y coordinate = ##1-\frac{1}{3}+\frac{1}{5}..=\frac{\pi }{4}##
converges to the point ##\left ( -1+\frac{log(2)}{2},\frac{\pi }{4} \right )##
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
918
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 5 ·
Replies
5
Views
1K
  • · Replies 16 ·
Replies
16
Views
909
  • · Replies 21 ·
Replies
21
Views
1K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
10K
  • · Replies 2 ·
Replies
2
Views
2K