0
|
20636 Views
|
22 Replies
|
8 Total Likes
View groups...
Share
GROUPS:

# How do I find x in a polynomial when x is approximate?

Posted 10 years ago
 I have difficulty programming to find x in a polynomial. I try NSolve to set the equation to zero and solve for x.If N = 85, x= 5, and y = 17I have an equation:y = sqrt((N * y  x^2)/x)as an approximation of the larger Prime product (or any possible product)My equation as it stands does not find x. I know there is a margin of error in the equation. This error can be subtracted to get a closer approximation. For example:In the following equations: y = 17.2941176471 = 16.8522995464 This is the answer when subtracting the 2 equations below to find x knowing only N. Look at the Yellow.So as the equations approach 0 (plus the margin of error). So since it is an approximation I need to find those values that are less than 1 and greater than 0. But 16.9082 is the value of the equation and not the value of x which should be found.I need help programming this. So if anyone can help it would be much appreciated.y = sqrt[(((85/x) * 85 - x^2)/ x) ]= ((85^2/x) + x^2)/ 85p = ((((85^2/x) + x^2)/ 85 ) ^2)- (((85/x) * 85 - x^2)/ x )^2 - 0.4418181007  ;sol = NSolve[p?0, x] {{x?-86.8953},{x?82.8717},{x?23.3122},{x?-7.50642+19.0312 ?},{x?-7.50642-19.0312 ?},{x?-10.592+14.3276 ?},{x?-10.592-14.3276 ?},{x?16.9082}}y = sqrt[(N * y - x^2) / x] y = sqrt[(7872197 * 3191 - 2467^2) / 2467]You can see this work at http://www.3dbuzz.com/forum/threads/200441-New-One-Way-Function/page2Look at the last page first.Also my website is www.constructorscorner.net
22 Replies
Sort By:
Posted 10 years ago
Posted 10 years ago
Posted 10 years ago
 I know what you are thinking after the last equation: "You call yourself Trurl the Constructor and  you throw a junk equation into an already hard to follow problem?"Well I will show you a relationship between Prime numbers that doesn't take Mathematica to calculate some high degree polyonmial. But you may say, "There are infinitely many of such relationships. What makes this one significant?"Nothing. An equation is only as significant as you make it. I just saying we are taught that Prime numbers are impossible to work with and I say they are not.Where x = 571 and y = 13811381 / 2\ = 219.79297642\ * (remainder(1381/2\) * 2\) * 571 = 13208.65186rewritten: 2\ * ((1381/2\-whole number part) * 2\) * 571 = 13208.6518613208.65186 / (1381/2\) = 60.0913208.65186 / 60.09 = 219.7929764219 = 219So judge for yourself if this is significant. This is very confusing so I hope I wrote this down right. But as I once invented a machine that made all things "N'. Numbers. I am in the business of constructing numbers.If the last equation didn't work logically. But sometimes you plug numbers into a calculator along the lines of your equation and it gives you a "feel" for what is happening with those numbers.________________________________________________________________________________________________________________________________New equations:2\ * (remainder (17 / 2\)) = 4.424.42 * 5 = 22.122.1 * (5 / (17* 2\)) = 0.8840.884 * 5 = 4.424.42 = 4.42
Posted 10 years ago
 only derived by logic not tested:    ((85/\) / 2\) = x*(remander of (y/2\)) where N= 85
Posted 10 years ago
 A simple trigonometric definition:s = R\A vector where x + y = NBut this isn't useful in this form so we solve with a triangle that is similar that we do know some values."x" is the value we are solving for. We know N - ((y/2\). The obtuse angle between the 2 "x" radiuses is (y/2\). This is if we hadtaken the angle of y and rotated it around the cirlce with radius "x". So "x" * "y" equals N. And as in the trigonometric definition s = R\.But we know that the smaller, similar triangle n = (N/x) / 2\) by inspection. (Where N is the given product of 2 Prime numbersand the angle wraps around the circle is divided by 2\. And we know that y = (N/x).So by vector division u + v = the longest chord or N. So we form the following equation:|tan(y/(2\))| * x + x = (y/(2\))substitute (N/x) for y|tan((N/x)/(2\))| * x + x = ((N/x)/(2\))multiply both sides by x |tan(N//(2\))| * x + x^2 = (N/(2\)) example where N = 85 |tan(85//(2\))| * x + x^2 = (85/(2\))
Posted 10 years ago
 "Do you believe that the last equation solves for x knowing only N?" The last few equations in the thread do not have an N.  Find the equation you mean. If it has only x's, N's, and numbers, then there could probably be at least an implicit numerical solution. I don't know now if a symbolic solution would exist.
Posted 10 years ago
 I truly believe this problem has merit. It took 2 seconds to dream it up as Prime numbers being wrapped around a logarithmic spiral and 5+ years to write it all down. I do not want to dwell on this problem. I have spent a lot of time on it. My next step is to show how the logarithmic spiral should work.I would not have an estimate of which numbers to try as the Prime products if it were not for Wolfram Alpha (Mathematica).I sent it to a library researches and they basically said to call them back when I can factor a 256 bit number. So there is still work to test to see if it works. But I see something here and hope others see it too.I was hoping this had merit and went viral. So far I estimate over 10,000 views. But it wouldnt matter if it doesnt work.@ Bruce MillerYou showed me how to find x knowing y in Mathematica code. Do you believe that the last equation solves for x knowing only N?
Posted 10 years ago
 Sqrt[(85*(85/x)  x^2)/ x] = yAbs(tan(85 / (2*Pi))) * x + x = 85/(2*Pi) tan in radians of course.((85^4/x + 2*(85^2 * x^2) + x^5) / 85^3 ) * x  85 = 0________ These three equations solve show how two Prime numbers x and y equal 85.This is important because in asymmetric cryptography such as RSA, the public key is the product of 2 unknown Prime numbers. If we can set a relationship to find x or y when given N, which in this case is 85, the one way function that RSA is based on is mathematically defeated!!!_______ So if we set equation 1 equal to equation 3 the polynomial produced, if solvable, will solve for x.Equation used to find 5 knowing 85:((85^4/x + 2*(85^2 * x^2) + x^5) / 85^3 )  0.6 = Sqrt[(85*(85/x)  x^2)/ x] + 0.15= approxiamtely 5.02 when solved by Wolfram Alpha!I believe this will defeat RSA. However no one believes me or just doesnt care. To me it is interesting, but when I ask people about it no response. If it worked I think it would go viral. What do you think?
Posted 10 years ago
Posted 10 years ago
 I don't follow most of this but I can answer the last question. No, it (NSolve-based methods) cannot beprogrammed to factor the product of two large primes.
Posted 10 years ago
 Daniel,This solution is deceiving. The angle is in radians of course. But does it work for other values of N?The answer is not exactly. There is some range of error. Also I made a algebraic mistake in the simplification.But if this is fixed and there are test values it does seem to have some use. But solving for N proves difficult.My question is does the following show any benefit from my effort? y = sqrt[(((85/x) * 85 - x^2)/ x) ] = ((85^2/x) + x^2)/ 85  In:=  p = ((85/x) * 85 - x^2)/ x - (( 85^2/((85^2/x) + x^2) ) ^2);  sol = NSolve[p == 0, x]  Out= {{x -> -36.2894}, {x -> 27.7376 + 21.7226 I}, {x ->     27.7376 - 21.7226 I}, {x -> -10.093 + 32.5167 I}, {x -> -10.093 -     32.5167 I}, {x -> 9.44493}, {x -> 0.257048 + 9.23565 I}, {x ->    0.257048 - 9.23565 I}, {x -> -8.95875}}This is still an approximation. In the above example 9 solves the equation better than 5. That is where 5 should be the desired answer. N = 85; x = 5; and y = 17 . So there still isn't a perfect solution just estimated. But take:In:= y = sqrt[(N * y - x^2) / x]y = sqrt[(7872197 * 3191 - 2467^2) / 2467]Out= sqrt[(-x^2 + N sqrt[(-x^2 + N sqrt[(-x^2 + 1/85 N (7225/x + x^2))/x])/x])/x]Out= sqrtSo sqrt of 10180014 = 3190.613421 which is approx y or 3191.The above substitution shows it works for larger values of N. Can it be programmed to factor the product of 2 large Prime Numbers?
Posted 10 years ago
 Looks like a reasonable approximation to a root.In:= FindRoot[Abs[Tan[(85/(2*Pi))]]*x + x - 85/(2*Pi) == 0, {x, 2}]Out= {x -> 5.55835870667}
Posted 10 years ago
 http://www.wolframalpha.com/input/?i=Abs%28tan%2885+%2F+%282*Pi%29%29%29+*++x+%2B+x+%3D+85%2F%282*Pi%29Abs(tan(85 / (2*Pi))) *  x + x = 85/(2*Pi)comes to 5.5 in Wolframalpha.Is this anything?
Posted 10 years ago
 Using vector addition I have another equation. I can keep making these and maybe I will find a helpful relationship. But it does no good to find an equation in polynomial form if I cannot solve for x in the polynomial.This time I keep it simple. My question is does this accomplish anything?Abs[tan(y / (2*Pi))]*x+x=(y ? ((2*Pi)))
Posted 10 years ago
 If the goal really is to find the integer factors, perhaps one could just use Solve or FactorInteger. In:= Solve[x y == 85 && 1 < x < y, {x, y}, Integers]Out= {{x -> 5, y -> 17}}In:= FactorIntegerOut= {{2467, 1}, {3191, 1}}
Posted 10 years ago
 Bruce can you work your Mathematica magic on this?Sqrt[((85*(85/x)) - x^2)/x] - (( 85^2/((85^2/x) + x^2) ) ^2) = 0Plugged into Wolfram Alpha this yields for x :  x is approximately 4.42297Getting pretty close to 5 not knowing y!This is a margin of error in the 2 equations that are subtracted about 0.599Is this any closer to solving knowing only N?
Posted 10 years ago
 To get a value for x, you need y.  With just n, you can get a relationship between x and y.  In:= soln = Solve[y == Sqrt[(n*y - x^2)/x], x] Out=         1    2                 4           1    2                 4 {{x -> - (-y  - Sqrt[4 n y + y ])}, {x -> - (-y  + Sqrt[4 n y + y ])}}        2                                  2  In:= soln /. {n -> 85, y -> 16.85229955} Out= {{x -> -288.957}, {x -> 4.95729}} In:= soln /. n -> 85Out=        1    2                 4           1    2                 4{{x -> - (-y  - Sqrt[340 y + y ])}, {x -> - (-y  + Sqrt[340 y + y ])}}       2                                  2In:= \$VersionOut= "9.0 for Mac OS X x86 (64-bit) (January 24, 2013)"In this plot, y is the horizontal axis and x is the vertical axis. In:= Plot[Evaluate[x /. soln /. n -> 85], {y, -30, 30},     PlotStyle -> {{Blue, Thick}, {Red, Thick}}, Frame -> True, GridLines -> {None, {0}}, FrameTicks -> All] Posted 10 years ago
 In:= soln = Solve[y == Sqrt[(n*y - x^2)/x], x]soln //. {n -> 85, y -> 16.85229955} // NOut= {{x -> 1/2 (-y^2 - Sqrt[4 n y + y^4])}, {x ->    1/2 (-y^2 + Sqrt[4 n y + y^4])}}Out= {{x -> -288.957}, {x -> 4.95729}}Bruce I am in your debt. You helped show the equation works. I moved a parenthesis to correct.But the equation works. Unfortunately there is still the difficult task of solving the polynomial only knowing N. I isolated x, but does that help knowing only N?The problem is this. I know N, I dont know that we can solve polynomials. Because I know N and say if 17 where tried it would equal 5.I have an equation with x unknown, but the only variable in the equation is x. So interval proves those values are true. But can we solve the polynomial for x only knowing N?
Posted 10 years ago
 y = sqrt[(((85/x) * 85 - x^2)/ x) ] = ((85^2/x) + x^2)/ 85p = ((((85^2/x) + x^2)/ 85 ) ^2) - (((85/x) * 85 - x^2)/ x )^2 -    0.4418181007  ;sol = NSolve[p == 0, x]{{x -> -86.8953}, {x -> 82.8717}, {x ->    23.3122}, {x -> -7.50642 + 19.0312 I}, {x -> -7.50642 -     19.0312 I}, {x -> -10.592 + 14.3276 I}, {x -> -10.592 - 14.3276 I}, {x ->    16.9082}}This is still an approximation. In the above example 9 solves the equation better than 5. That is where 5 should be the desired answer. N = 85; x = 5; and y = 17 . So there still isn't a perfect solution just estimated. But take:y = sqrt[(N * y - x^2) / x]y = sqrt[(7872197 * 3191 - 2467/^2) / 2467]
Posted 10 years ago
 Sorry for the delay. I much value your opinion. I have been trying for 6+ years in my spare time to find a pattern of Prime numbers in a logarithmic spiral. This here is not a pattern in Prime numbers, but a description of multiplication.The N = p *eI am saying knowing N I can find p the smaller product. My example is 85 = 5 * 17The trouble is it is a polynomial which means that the equation is difficult if not impossible to solve. Better methods of solving polynomials need to be found. I have not learned higher mathematics  beyond linear algebra, so I rely on Mathematica which crunches numbers.The trouble is both sides of the equation are equal. That is, they both equal  17. But my equation when solving for x = 17 for the single case of 85 = 5 * 17. What I dont know is why I plug 5 in for x and 85 for N it equals (with a margin of error) 17.Is my method of programming Mathematica wrong? Why cant I solve for x?I hope this clears up my question, but if it is still unclear let me know. The code is posted below.BTW I am looking at interval now. Thanks for the direction.
Posted 10 years ago
 Have you tried Interval?  In:= soln = Solve[ y == Sqrt[(n*y - x^2)]/x, x]                                                         Sqrt[n] Sqrt[y]          Sqrt[n] Sqrt[y] Out= {{x -> -(---------------)}, {x -> ---------------}}                             2                        2                   Sqrt[1 + y ]             Sqrt[1 + y ]  In:= soln //. {n -> 85, y -> 17} // N                                                 Out= {{x -> -2.23221}, {x -> 2.23221}}In:= ?Interval                                                                       Interval[{min, max}] represents the range of values between min and max. In:= soln  //. {n -> 85, y -> Interval[{16.8523, 17.294}]}                           Out= {{x -> Interval[{-2.27109, -2.18484}]}, {x -> Interval[{2.18484, 2.27109}]}}
Posted 10 years ago
 The code as presented in the first web page you  link to is easier to read. The -> for Rule and I for Sqrt[-1] come across better than special characters. Using a Mathematica built-in function name (N) as a user variable name is asking for trouble. It wasn't clear from the code above or the link if "sqrt" is your own function or a mis-capitalized Sqrt.