1

I am trying to find the optimum solution for calculating and visualizing the critical point of an inverse power curve.

As an example, consider the curve of: $150x^{-0.4}$

link to image of curve

I modeled the function using https://www.geogebra.org/graphing (type in $150x^-0.4$ as the function). A power curve is unlike a polynomial because it does not have a sudden vertex where the rate of change $=0$. Therefore I need to use a different method of calculating the maximum or minimum (depending on orientation) for the curve.

Visual inspection of this curve shows that maximum curvature in the 'elbow' is likely somewhere between $10<x<60$. I used recommendations for finding maximum curvature based on protocol provided in other forum threads:

The answer I calculated as the point of maximum curvature is (50.7, 31.2). However, I am not confident about this answer because I calculated curvature at each point and the answer is between 16 and 18 based on my calculation of curvature in excel. I don't know how to validate it though.

table of curve values. note change at 16-18

I have a couple inquiries:

  • Is this the equation I should be using to find the apex or maximum of this curve?
  • What is recommended software or programming that would aid in quantifying and visualizing this maximum curvature. It would be nice to have software like what is seen 52 seconds into this video: https://youtu.be/wyPXbvsd9nI?t=52 (other examples here, but none of them have been satisfying after a quick test. I'm still exploring)
wayne r
  • 13

1 Answers1

1

Here are my calculations, following the method given in the linked posts.

General power function: $$y=ax^p$$ Derivative: $$dy/dx=apx^{p-1}$$ Second derivative: $$d^2y/dx^2=ap(p-1)x^{p-2}$$ Curvature, squared: $$k=\frac{d^2y/dx^2}{\sqrt{1+(dy/dx)^2\,}^3}$$ $$k^2=\frac{(d^2y/dx^2)^2}{(1+(dy/dx)^2)^3}$$ $$=\frac{a^2p^2(p-1)^2x^{2p-4}}{(1+a^2p^2x^{2p-2})^3}$$ Derivative of curvature vanishes at maximum: $$0=d(k^2)/dx=\frac{a^2p^2(p-1)^2(2p-4)x^{2p-5}}{(1+a^2p^2x^{2p-2})^3}-3\frac{a^2p^2(p-1)^2x^{2p-4}}{(1+a^2p^2x^{2p-2})^4}a^2p^2(2p-2)x^{2p-3}$$ Dividing by $2a^2p^2(p-1)^2$: $$0=\frac{(p-2)x^{2p-5}}{(1+a^2p^2x^{2p-2})^3}-3\frac{x^{2p-4}}{(1+a^2p^2x^{2p-2})^4}a^2p^2(p-1)x^{2p-3}$$ Dividing by $x^{2p-5}$ and multiplying by $(1+a^2p^2x^{2p-2})^4$: $$0=(p-2)(1+a^2p^2x^{2p-2})-3(x^{2p-2})a^2p^2(p-1)$$ $$0=(p-2)+a^2p^2(p-2-3p+3)x^{2p-2}$$ $$x^{2p-2}=\frac{(p-2)}{a^2p^2(2p-1)}$$ $$x=\left(\frac{(p-2)}{a^2p^2(2p-1)}\right)^{1/(2p-2)}$$ Here's our general solution. Plugging in the values $a=150$ and $p=-0.4$, I get $$x=2700^{5/14}\approx16.8069,\quad y=150\cdot2700^{-1/7}\approx48.5173$$

mr_e_man
  • 5,986
  • Are you willing to explain your answer a bit more? Specifically where you calculated the derivative of the curvature equation? I eventually figured out that the answer was somewhere between 16 and 17 by calculating curvature at each point, but I wasn't able to find the correct general solution. I think the fault in my answer is with finding the derivative of the curvature equation. – wayne r Jun 30 '22 at 16:02
  • The curvature has the general form $k^2=u(1+v)^{-3}$. By the usual rules for differentiation, $$d(k^2)/dx=(du/dx)(1+v)^{-3}+u(-3)(1+v)^{-4}(dv/dx)$$ – mr_e_man Jun 30 '22 at 16:05
  • $$u=a^2p^2(p-1)^2x^{2p-4},\quad v=a^2p^2x^{2p-2}$$ – mr_e_man Jun 30 '22 at 16:07
  • Thank you. I used the wrong derivative form. I mixed the dividing rule with the multiplication rule. – wayne r Jun 30 '22 at 16:20