- Thread starter
- #1

- Thread starter Wilmer
- Start date

- Thread starter
- #1

- Admin
- #2

- Jan 26, 2012

- 4,198

Several thoughts:

1. You could quite easily eliminate one variable, as all three equations are explicit. [EDIT]: See Bacterius below for a better answer.

2. You're most definitely not guaranteed a unique solution. The floor function is not 1-1, so inverting it can be problematic. You're not guaranteed that a solution exists, and if one does exist, you're not guaranteed that it is unique.

3. For a highly nonlinear system like this, I'm thinking numerical is your only option. The Excel Solver routine might give you good results, or maybe MATLAB or a MATLAB clone like GNU Octave if you can't afford MATLAB.

1. You could quite easily eliminate one variable, as all three equations are explicit. [EDIT]: See Bacterius below for a better answer.

2. You're most definitely not guaranteed a unique solution. The floor function is not 1-1, so inverting it can be problematic. You're not guaranteed that a solution exists, and if one does exist, you're not guaranteed that it is unique.

3. For a highly nonlinear system like this, I'm thinking numerical is your only option. The Excel Solver routine might give you good results, or maybe MATLAB or a MATLAB clone like GNU Octave if you can't afford MATLAB.

Last edited:

- Jan 26, 2012

- 644

- Admin
- #4

- Jan 26, 2012

- 4,198

Good catch!

- Jan 26, 2012

- 890

Change the variables to a=1/A, B and C and the solution is a=0, B=-2, C=-2

(or my algebra could have gone horribly wrong, you could also help yourself by simplifying the equations)

CB

- Thread starter
- #6

I tried it, REMOVING the floor function, this way:

Let u = 1.25 * 10^(-9) / (27/2) , v = 1.29 , w = 4.476

Then the equations can be rewritten:

A = u * [(B + 2)(C + 2)]^2

B = v / A^(1/4) - 2

C = w / A^(1/4) * SQRT[(B + 2) / (11B + 43)] - 2

Manipulations/contortions(!) of above lead to:

11vA^2 + 21A^(9/4) - uv^3w^2 = 0 ; (I'm 99.9% sure that's correct!)

And that also seems to require numeric solving; you were correct in "many solutions";

Wolfram hands out 8 of them; with one real where A = .0000160191 :

http://www.wolframalpha.com/input/?i=11*1.29*a%5E2%2B21*a%5E%289%2F4%29-%281.25*10%5E%28-9%29%2F13.5%29*1.29%5E3*4.476%5E2%3D0

These equations I'm told are for buckling of columns, in structural engineering!

Further comments appreciated!

- Admin
- #7

- Jan 26, 2012

- 4,198

We can simplify the expressions as follows:

\begin{align*}

A&=\frac{5}{54}\,\times 10^{-9}(B+2)^{2}(C+2)^{2}\\

B&=\left\lfloor \frac{1.29}{\sqrt[4]{A}}\right\rfloor-2\\

C&=\left\lfloor\frac{4.476}{\sqrt[4]{A}}\sqrt{\frac{B+2}{11B+43}}\right\rfloor

-2

\end{align*}

Now we follow Bacterius's suggestion to obtain the following single equation

for $A$:

$$A=\frac{5\times 10^{-9}}{54}\left\lfloor\frac{1.29}{\sqrt[4]{A}}\right\rfloor^{2}

\times\left\lfloor\frac{4.476}{\sqrt[4]{A}}\sqrt{\frac{\left\lfloor\frac{1.29}{\sqrt[4]{A}}\right\rfloor}{11\left\lfloor\frac{1.29}{\sqrt[4]{A}}\right\rfloor+21}}\right\rfloor^{2}.$$

If you plot the LHS and the RHS on the same axis, you do get an intersection

on a horizontal portion of the RHS's graph at around $A=0.000015$. Using the

FindRoot command in Mathematica yielded the solution

$$A\to 0.0000150371,\quad B\to 18.,\quad C\to 18.$$

Naturally, $B$ and $C$ must be integers.

If I use the FindRoot command just on the $A$ equation, I get a more refined value of $A\to 0.0000149558.$

\begin{align*}

A&=\frac{5}{54}\,\times 10^{-9}(B+2)^{2}(C+2)^{2}\\

B&=\left\lfloor \frac{1.29}{\sqrt[4]{A}}\right\rfloor-2\\

C&=\left\lfloor\frac{4.476}{\sqrt[4]{A}}\sqrt{\frac{B+2}{11B+43}}\right\rfloor

-2

\end{align*}

Now we follow Bacterius's suggestion to obtain the following single equation

for $A$:

$$A=\frac{5\times 10^{-9}}{54}\left\lfloor\frac{1.29}{\sqrt[4]{A}}\right\rfloor^{2}

\times\left\lfloor\frac{4.476}{\sqrt[4]{A}}\sqrt{\frac{\left\lfloor\frac{1.29}{\sqrt[4]{A}}\right\rfloor}{11\left\lfloor\frac{1.29}{\sqrt[4]{A}}\right\rfloor+21}}\right\rfloor^{2}.$$

If you plot the LHS and the RHS on the same axis, you do get an intersection

on a horizontal portion of the RHS's graph at around $A=0.000015$. Using the

FindRoot command in Mathematica yielded the solution

$$A\to 0.0000150371,\quad B\to 18.,\quad C\to 18.$$

Naturally, $B$ and $C$ must be integers.

If I use the FindRoot command just on the $A$ equation, I get a more refined value of $A\to 0.0000149558.$

Last edited:

- Admin
- #8

- Jan 26, 2012

- 4,198

- Thread starter
- #9

- Admin
- #10

- Jan 26, 2012

- 4,198

You get an exact value from plugging $B=18$ and $C=18$ into the first equation. There's no range. I get $A=1/67500$. How are you getting your range?Thanks.

However, B=18 and C=18 results in range A = 1.43000 * 10^(-5) to 1.72000 * 10^(-5)

Agree?

- Thread starter
- #11

Sorry...was doing something different...not applicable here...How are you getting your range?

Thanks everyone for the responses. All's fine now.