[ home ] [ math / cs / ai / phy / as / chem / bio / geo ] [ civ / aero / mech / ee / hdl / os / dev / web / app / sys / net / sec ] [ med / fin / psy / soc / his / lit / lin / phi / arch ] [ off / vg / jp / 2hu / tc / ts / adv / hr / meta / tex ] [ chat ] [ wiki ]

/math/ - Mathematics


Name
Email
Subject
Comment
Verification
Instructions:
  • Press the Get Captcha button to get a new captcha
  • Find the correct answer and type the key in TYPE CAPTCHA HERE
  • Press the Publish button to make a post
  • Incorrect answer to the captcha will result in an immediate ban.
File
Password (For file deletion.)

25 Dec 2021Mathchan is launched into public


File: Saudi combined e.png ( 334.65 KB , 2033x1433 , 1663822754492.png )

Image
Peak oil.

It's mathematically impossible for you to not die within a year.

I posted with the alias mustang on pol, you can see in archive.


File: r_triangle_1.jpg ( 4.38 KB , 242x208 , 1655398539385.jpg )

Image
Trigonometric functions are used to relate the angles of a right triangle to its sides.

sin(angle)=oppositehypotenusecos(angle)=adjacenthypotenusetan(angle)=sin(angle)cos(angle)=oppositeadjacent\qquad \sin(\text{angle}) = \frac{\text{opposite}}{\text{hypotenuse}}\quad\cos(\text{angle}) =\frac{\text{adjacent}}{\text{hypotenuse}}\quad\tan(\text{angle}) = \frac{\sin(\text{angle})}{\cos(\text{angle})} = \frac{\text{opposite}}{\text{adjacent}}\quad


With respect to any of the two angles
α\alpha
or
β\beta
(the third angle of a right triangle is, of course,
90o90^\text{o}
thereofre it's irrelevant) three sides can be identified with the following names: hypotenuse, opposite and adjacent. The hypotenuse is easy to identify and it's the longest i.e.
cc
. The adjacent side is the one that's "touching" the angle; for angle
α\alpha
that would be
bb
, while for angle
β\beta
that would be
aa
. The opposite side is the one that's NOT "touching" the angle; for angle
α\alpha
that would be
aa
, while for angle
β\beta
that would be
bb
.

Therefore, the complete set of trigonometric functions for the triangle in the pic related is:

sin(α)=accos(α)=bctan(α)=ab\qquad \sin(\alpha) = \frac{a}{c}\quad \cos(\alpha) = \frac{b}{c}\quad \tan(\alpha) = \frac{a}{b}


sin(β)=bccos(β)=actan(β)=ba\qquad \sin(\beta) = \frac{b}{c}\quad \cos(\beta) = \frac{a}{c}\quad \tan(\beta) = \frac{b}{a}



Thus, just knowing the angle
α\alpha
allows conversion of any side to any other side:
aa
bb
cc
a=a=
aa
btan(α)b\tan(\alpha)
csin(α)c\sin(\alpha)
b=b=
atan(α)\frac{a}{\tan(\alpha)}
bb
ccos(α)c\cos(\alpha)
c=c=
asin(α)\frac{a}{\sin(\alpha)}
bcos(α)\frac{b}{\cos(\alpha)}
cc

And similarly, just knowing the angle
β\beta
also allows conversion of any side to any other side:
aa
bb
cc
a=a=
aa
btan(β)\frac{b}{\tan(\beta)}
ccos(β)c\cos(\beta)
b=b=
atan(β)a\tan(\beta)
bb
csin(β)c\sin(\beta)
c=c=
acos(β)\frac{a}{\cos(\beta)}
bsin(β)\frac{b}{\sin(\beta)}
cc

Intuition behind this is straightforward: For angles of a triangle, sine and cosine are bounded functions between zero and one.
(0sin(x)1,  0cos(x)1)(0 \leq \sin(x)\leq 1,\; 0\leq \cos(x) \leq1)
. Multiplying something with these functions should produce a smaller or equal value, while dividing something with these functions should produce a larger or equal value. Since hypotenuse is the longest, multiplying it with sine or cosine will "reduce" it to a leg. Similarly, dividing a leg by sine or cosine will "grow" it into a hypotenuse. A leg can be converted to another leg by "growing" it into a hypotenuse first, then "reduciPost too long. Click here to view the full text.
7 posts omitted. Click here to view.
>>
Taylor series of trigonometric functions

Sine and cosine have the following Maclaurin series (which are Taylor series around
x0=0)x_0 = 0)
:

sin(x)=n=0(1)n(2n+1)!x2n+1\qquad \sin(x) = \sum_{n=0}^\infty\frac{(-1)^n}{(2n + 1)!}x^{2n + 1}

cos(x)=n=0(1)n(2n)!x2n\qquad \cos(x) = \sum_{n=0}^\infty \frac{(-1)^n}{(2n)!}x^{2n}


This is the same as writing:

sin(x)=xx33!+x55!x77!+x99!x1111!+\qquad \sin(x) = x -\frac{x^3}{3!} + \frac{x^5}{5!} - \frac{x^7}{7!} + \frac{x^9}{9!} -\frac{x^{11}}{11!} +\dots

cos(x)=1x22!+x44!x66!+x88!x1010+\qquad \cos(x) = 1 - \frac{x^2}{2!} + \frac{x^4}{4!} - \frac{x^6}{6!} + \frac{x^8}{8!} - \frac{x^{10}}{10}+\dots


By summing the Maclaurin series of
cos(x)\cos(x)
and
isin(x)i\sin(x)
it's possible to prove Euler's formula
eix=cos(x)+isin(x)e^{ix} = \cos(x) + i\sin(x)


cos(x)+isin(x)=1+ix+(ix)22!+(ix)33!+(ix)44!+=eix\qquad \cos(x) + i\sin(x) = 1 + ix + \frac{(ix)^2}{2!} + \frac{(ix)^3}{3!} + \frac{(ix)^4}{4!} + \dots = e^{ix}


Maclaurin series can be used to approximate trigonometric functions to an arbitrary degree
For small angles (i.e.
α0\alpha\approx 0
), good approximations (often used in physics) are:

sin(α)x\qquad \sin(\alpha)\approx x

cos(α)1\qquad \cos(\alpha)\approx 1


Which are just Maclaurin series with only 1 term.
>>
Computation of trigonometric functions

Maclaurin series lend an easy way to numerically compute.
We need a helper function
ifac(n)
that computes
n!,  nNn!,\; n\in\mathbb{N}
and a helper function
fpow(x, n)
that computes
xn,  xR,nNx^n,\;x\in\mathbb{R}, n\in\mathbb{N}
:

int ifac(int n) {

return (n > 0) ? (n * ifac(n - 1)) : 1;
}
float fpow(float x, int n) {
return (n < 0) ? fpow(1/x, -n)
: (n == 0) ? 1
: (n == 1) ? x
: (n > 1) ? x * fpow(x, n - 1)
: 1;
}

Now
sin(x)\sin(x)
can be computed with
fsin(x, prec)
:
float fsin(float x, int termCount) {

int n, sign = 1;
float numer = x,
denom = 1,
result = sign * (numer / denom);

for (n = 1; n < termCount; n++) {
sign *= -1;
numer = fpow(x, 2*n+1);
denom = ifac(2*n+1);
result += sign * (numer / denom);
}

return result;
}

While
cos(x)\cos(x)
can be computed with
fcos(x, prec)
:

float fcos(float x, int termCount) {

int n, sign = 1;
float numer = 1,
denom = 1,
result = sign * (numer / denom);
for (n = 1; n < termCount; n++) {
sign *= -1;
numer = fpow(x, 2*n);
denom = ifac(2*n);
result += sign * (numer / denom);
}
return result;
}


While Maclaurin series for
sin(x)\sin(x)
and
cos(x)\cos(x)
converge on entire
R\mathbb{R}
, they converge the fastest when
x(π,π]x\in(-\pi,\pi]
so its good to bound
xx
to that interval considering trigonometric functions are periodic. In fact, only 5-15 are iterations are usually good enough in double-precision floating point precision.
>>

File: 8ce0ff601f99fdbf275b63c9b3b1944e.jpg ( 172.93 KB , 1040x1136 , 1656009606752.jpg )

Image
>>198
You touched on what nobody actually mentions, how these functions are computed. Excellent exposition.
>>

File: 1648530726809.png ( 25.95 KB , 644x800 , 1662284901464.png )

Image
gemmy though, thanks for the full info
>>
great thread


File: kotberet.png ( 154.62 KB , 393x403 , 1642718192578.png )

Image
As well as "Jew" or "gay" id on /pol/. I mean it's a three-letter word where second letter can be written as a number (k0t) and first and last letter aren't the same.

Now let's talk about IDs. They are 8-character sequences consisting of base64 characters (numbers, lowercase, uppercase letters, / and +). So, there can be

10+26+26+2=6410 + 26 + 26 + 2 = 64


possible characters for one position in id and thus

K0=648=248=281 474 976 710 656K_0 = 64^8 = 2^{48} = 281\ 474\ 976\ 710\ 656


possible IDs (BIG number :o).

Now back, to KOT.
You can write letters k, o and t in uppercase or lowercase and you can also write o as zero. That means you have
2×3×2=122\times3\times2 = 12
combinations of base64 characters that can make up the word "k0T".

Suppose that first three letters of your ID is one of 12 possible K0Ts. We have 5 letters of ID left and they can be anything from base64 alphabet.
So, there are

K1=12×645=12×230=12 884 901 888K_1 = 12 \times 64^5 = 12 \times 2^{30} = 12\ 884\ 901\ 888


So, it's also many, many ids. More than all cats, all humans, but less than all chickens in the world. So, maybe kot ID isn't really that rare.

Also there are 6 positions on which KoT id can happen:
KOTxxxxx
xKOTxxxx
xxKOTxxx
xxxKOTxx
xxxxKOTx
xxxxxKOT

We can assume that there are the same number of ids with KoT on n-th position where
n=1...6n = 1 ... 6
.

Remember, that KOT id can happen twice, so we must exclude duplicates:
KOTKOTxx
KOTxKOTx
KOTxxKOT
xKOTKOTx
xKOTxKOT
xxKOTKOT

Each of these duplicate positions exists in
K2=12×12×642=144×212=589824K_2 = 12 \times 12 \times 64^2 = 144 \times 2^{12} = 589824

configurations.

Because one duplicate can belong to two sets of IDs with at least one kot id (KOTxxKOT can belong to KOTxxxxx and xxxxxKOT), our equation for number of all KOT ids must be:

6×K16×K26 \times K_1 - 6 \times K_2


Let's calculate it:

6×K16×K2=6 \times K_1 - 6 \times K_2 =

6×12×2306×144×212=6 \times 12 \times 2^{30} - 6 \times 144 \times 2^{12} =

72×230864×212=72 \times 2^{30} - 864 \times 2^{12} =

77 305 872 38477\ 305\ 872\ 384


For every human on this planet there are approximately 10 KOT ids. :o

Now we can calculate probablility of having koT id:

77305872384K0=773058723842814749767106560.0002746456313641\frac{77305872384}{K_0} = \frac{77305872384}{281474976710656} \approx 0.00027464563 \approx \frac{1}{3641}


Thus, k0t id happens approximately one time for every 3641 ids.
Now you can calculate daily probability, because I don't know if ids on /bant/ are boardwise or just one for every poster and thread. x)
2 posts omitted. Click here to view.
>>
>>110
can you calculate the probability of "jak" id on /qa/ though
>>
>>126

Assuming that you can write second 'a' of 'jak' as 4 it's the same probability.

Thoughever, you can write it also as 'jaq' or 'jac', so there are 6 possibilities of writing the last letter there are:

2×3×6=362 \times 3 \times 6 = 36


Possible variation of previous constants for JaK:

K0=the sameK_0 = the\ same


K1=36×645=36×230=38 654 705 664K_1 = 36 \times 64^5 = 36 \times 2^{30} = 38\ 654\ 705\ 664


K2=36×36×642=1296×212=5 308 416K_2 = 36 \times 36 \times 64^2 = 1296 \times 2^{12} = 5\ 308\ 416



The final equation for the number of jAc IDs:

6×K16×K2=6 \times K_1 - 6 \times K_2 =

6×36×2306×1296×212=6 \times 36 \times 2^{30} - 6 \times 1296 \times 2^{12} =

216×2307776×212=216 \times 2^{30} - 7776 \times 2^{12} =

231 896 383 488231\ 896\ 383\ 488


And the probability:

231896383488K0=2318963834882814749767106560.00082386145311214\frac{231896383488}{K_0} = \frac{231896383488}{281474976710656} \approx 0.000823861453 \approx \frac{1}{1214}


So it's just pretty much 3 times bigger than "KoT". :o


/qa/ is dead tohough :(
>>
>>106
Let's eneralize this to id of length
nn
(for 4chan:
n=8n=8
) with set of
qq
different possible characers (for base64:
q=64q=64
) and substring of interest of length
 (n)\ell\ (\ell\leq n)
(for k0t:
=3\ell=3
) with
pp
different variations of the substring (for k0t:
p=12p = 12
)

The number of all possible IDs is:

K0=qn\qquad K_0 = q^n


If
k0t
appears in at least one place, we must consume
=3\ell = 3
characters to assemble them into
k0t
, then we return the
k0t
into the bunch (so we now have
n+1n -\ell + 1
items in the bunch) and then we arrange the
k0t
in
(n+11)n - \ell + 1\choose 1
ways. We have
pp
variations of k0t, plus
nn-\ell
other characters which can each be
qq
different vays (so
qnq^{n-\ell}
) (even tho I said "plus" we use
×\times
bcoz of thing called multiplication principle but w/e)

K1=(n+11)pqn\qquad K_1 = {n - \ell + 1 \choose 1}\cdot p \cdot q^{n - \ell}


If there are two
k0t
s somewhere, we must consume
2=62\ell = 6
characters to assemble 2
k0t
s (both with
p=12p=12
variations so
p2p^2
), then we return the k0t in the bunch (we now have
n2+2n -2\ell + 2
items) and we arrange the two kots in
(n2+2)2)n-2\ell + 2)\choose 2
ways. We have two kots which can each be
pp
different ways (so
p2p^2
) and
n2n - 2\ell
unconsoomed characters which can each be
qq
different ways (
qn2q^{n - 2\ell}
)

K2=(n2+2)2)p2qn2\qquad K_2 = {n - 2\ell + 2)\choose 2}\cdot p^2\cdot q^{n - 2\ell}


Therefore if we have
kk
k0t
s, the fomula is the following: (and I have a truly marevolous way to prove this formula by induction but I am also truly a lazy fuck.)

Kk=(nk(1)k)pkqnk\qquad K_k = {n - k(\ell - 1)\choose k}\cdot p^k\cdot q^{n-k\ell}


Now if we want to count all the
k0t
s without duplicates, we must use inclusion-exclusion principle:

KΣ=k=1nk(1)>0(1)k1Kk\qquad K_{\Sigma} = \sum_{k = 1}^{n - k(\ell - 1) > 0}(-1)^{k - 1}\cdot K_k


Then the probability of
k0t
appearing is

P=KΣK0\qquad P = \frac{K_\Sigma}{K_0}



So if I plug in the numbers for
k0t
on /pol/ slash /biz/:

n=8q=64=3p=12\begin{aligned} \qquad n &= 8\\q &= 64\\\ell&=3\\p&=12 \end{aligned}


We get the following result
P=KΣK0=1qnk=1nk(1)>0(1)k1Kk=1648((83+11)12645(86+22)122642)=12642648(6643612)=773058723842814749767106560.02746%\qquad P = \frac{K_\Sigma}{K_0} = \frac{1}{q^n}\sum_{k = 1}^{n - k(\ell - 1) > 0}(-1)^{k - 1}\cdot K_k = \frac{1}{64^8}\cdot\left({8 - 3 + 1\choose 1}\cdot 12\cdot 64^5 - {8 - 6 + 2\choose 2}\cdot 12^2\cdot 64^2\right) = \frac{12\cdot 64^2}{64^8}\cdot(6\cdot 64^3 - 6\cdot 12) = \frac{77305872384}{281474976710656} \approx 0.02746\%


Which matches OP's numbers!
>>
>>129
Nice, I treat the fact that our results are the same as a miracle because I always get something wrong during calculations. x)

Also I wonder if determining number for arbitrary word and
nn
(or just equation for
Kk,k1,...,nK_k, k \in 1,...,n
) can have an elegant equation or generic brute-force algorithm is the best way.

For example the word SUS is kinda sus, because SUSUSUSx is a valid id with three occurences of the searched word.

We can assume that there are problems when
kk
-long suffix of the word is
kk
-long prefix like when 'S' is 1-long suffix and prefix of 'SuS' and I see hope there, but words like 'SSUSSUSS' have the same words of length 1,2 and 5 as prefixes and suffixes of respective lengths.
>>
test


File: 1638294292313.png ( 363.02 KB , 1195x1240 , 1640551533576.png )

Image
Supose you have two baskets with one ball each. Probability that one basket has k ball in 0-th time is

p0,kp_{0,k}


therefore

p0,1=1p_{0,1} = 1


and for k != 1

p0,k=0p_{0,k} = 0
.

Now (You) choose random ball from any basket (every ball have the same probability of being chosen, not dependent on which basket is it placed in). You put new ball in a basket where this choosen ball was. :0 You repeat this t times.

In math language it is probably:

pt,k=pt1,k1k1t+1+pt1,ktkt+1p_{t,k} = p_{t-1,k-1}\frac{k-1}{t+1} + p_{t-1,k}\frac{t-k}{t+1}


Now

The most bone-chilling, slow-burn, atmosphere-oozing thing I discovered about it:

pt,k=1t+1p_{t,k} = \frac{1}{t+1}


HOLY SCIENCE!!!!!!!!1

Probability of having k balls after t time is always the same!!!!!
5 posts omitted. Click here to view.
>>
>Probability of having k balls after t time is always the same!!!!!
I agree. And this probability is 1/2. It either happens or it doesn't.
>>
>>254
and now consider the conditional probability of it happening or not.
>>
>>73
op is a faggot (test)
>>
>>269
Your containment board >>>/test/
>>
>>73
Your post was kind of difficult to read so I will try to summarise the problem and solution in my own words. I also think you made a typo, so my reccurance relation looks very slightly different to yours, but the solution is the same.


Problem

You have two boxes, box #1 and box #2. You also have a bag with an unlimited number of balls. At time 0, box #1 has one ball and box #2 has one ball.

In every iteration:
- Choose a ball at random from the set of balls that are in the boxes.
- Take a new ball from the bag and put it into the box where the chosen ball was.

So at time
tt
, there are
t+2t+2
balls in both boxes. Each box has between
11
and
t+1t+1
balls (inclusive).

Let
pt,kp_{t,k}
denote the probability that box #1 has exactly k balls at time t.

So by the initial conditions:
-
p0,k=1p_{0,k}=1
if
k=1k=1

-
p0,k=0p_{0,k}=0
otherwise

Derive a closed form expression for
pt,kp_{t,k}
.


Solution

We know that at time
t1t-1
there were
(t1)+2=t+1(t-1)+2=t+1
balls in both boxes.

At time
tt
, how can box #1 come to contain
kk
balls? There are two cases:
- at time
t1t-1
it had
k1k-1
balls and was chosen
- at time
t1t-1
it had
kk
balls and was not chosen

This leads to the following reccurence relation

pt,k=pt1,k1k1t+1+pt1,kt+1kt+1p_{t,k} = p_{t-1,k-1}\frac{k-1}{t+1} + p_{t-1,k}\frac{t+1-k}{t+1}


which has the solution by induction on
tt


pt,k=1t+1p_{t,k} = \frac{1}{t+1}


where
1<=k<=t+11<=k<=t+1
.

Base case for
t=1t=1
is trivial.

pt,kp_{t,k}


=pt1,k1k1t+1+pt1,kt+1kt+1 = p_{t-1,k-1}\frac{k-1}{t+1} + p_{t-1,k}\frac{t+1-k}{t+1}


=1(t1)+1k1t+1+1(t1)+1t+1kt+1 = \frac{1}{(t-1)+1}\frac{k-1}{t+1} + \frac{1}{(t-1)+1}\frac{t+1-k}{t+1}


by inductive assumption and because
t>=1t>=1
so no problems with dividing by 0

=k1t(t+1)+t+1kt(t+1) = \frac{k-1}{t(t+1)} + \frac{t+1-k}{t(t+1)}


=k1+t+1kt(t+1) = \frac{k-1+t+1-k}{t(t+1)}


=tt(t+1) = \frac{t}{t(t+1)}


=1t+1 = \frac{1}{t+1}


which is what we wanted.


Now for some values of
kk
, namely
11
and
t+1t+1
this is somewhat an obvious result because in these cases, we have to look at the probability of choosing one box repeatedly which is just
t!(t+1)!=1t+1\frac{t!}{(t+1)!}=\frac{1}{t+1}
, but for the other cases it is not so obvious. Drawing a tree diagram helps verify the result, I did it up to
t=3t=3
.

Cool problem overall anon :)


File: Oyu.jpg ( 52.67 KB , 960x720 , 1657222704337.jpg )

Image
What are some good books to learn symplectic geometry with?
>>
Have you had a look at Abraham and Marsden? Alternatively, Berndt?
>>
could symplectic geometry replace a classical mechanics course?
>>
>>252
No. You would lack almost all physical understanding of what is happening. It would likely seem very poorly motivated too, if you didn't know where the idea of a symplectic manifold came from and why. I'm sure some people do study symplectic geometry without studying classical mechanics, but it seems pointless to me. And to actually replace a course on classical mechanics (the implication being that this is someone studying physics), I think would be a bad idea.


File: dbqOJ.gif ( 239.17 KB , 1838x1304 , 1658138309873.gif )

Image
The picture below shows (in red) a sum of Gaussian kernels with different means, but all having variance 1, i.e.
f(x)=12π−−√∑i=1ne−12(x−μi)2 .
. The green curves show what I get if I select a point x0, and then "fit" a scaled Gaussian distribution g to f about x0, in the sense that I determine c and μ0 in the function
g(x):=c2π−−√e−12(x−μ0)2
so that g(x0)=f(x0) and g′(x0)=f′(x0). From numerical experimentation, it seems to always be the case that g(x)≤f(x) for all x. Does anyone know why this is so?


File: IMG_20220519_220921.jpg ( 3.11 MB , 4000x3000 , 1653047795771.jpg )

Image
Are there mathemathics i can use to fix the silhouette of this image? like the back tip of the helmet?
to make it look "correct"?


File: image.png ( 126.74 KB , 2048x2048 , 1643699011807.png )

Image
does anyone have good resources to learn the notation? i could figure out basic things like
((λx.x)(λa.a))((λ x . x) (λ a . a))

but got filtered by
(((λf.(λx.(fx)))(λa.a))(λb.b))(((λ f . (λ x . (f x))) (λ a . a)) (λ b . b))
>>
>>134
What's your background
>>
>>138
mathematics, i managed to figure out the notation by fucking around with emacs lisp
>>

File: da9ffae136e360d515303cbbb8af99e1.jpg ( 350.26 KB , 744x1052 , 1644299884658.jpg )

Image
>>139
https://www.mscs.dal.ca/~selinger/papers/papers/lambdanotes.pdf
>>
>>134
Usually repeated abstraction is notated as just
lam f x. ...
instead of
lam f. lam x. ...
Lambda calculus isn't usually used in practice.
Not sure what you want to use it for.


File: 00000.jpg ( 81.33 KB , 289x293 , 1649020837918.jpg )

Image
how do i start learning advanced math, ahead of my classes? any good textbooks/sources? what is your advice?
>>

File: orthocenter.png ( 329.63 KB , 1521x840 , 1649977152884.png )

Image
>>159
Mathematics at an advanced level diverges into many different topics, each something you could study for years. What do you want to know more about?
>>

File: 1636870494174.jpg ( 202.36 KB , 2048x2048 , 1650886614985.jpg )

Image
If you are just saying something as vague as "learning advanced math", then we can only tell you to google good introductory textbooks on various subjects. Maybe first open up Wikipedia and just check different genres there and see if something clicks and then find intro books on that!
>>
A graphical approach to Algebra and Trigonometry basically gets you into Calculus


File: Tensegrity_simple_3_RL.png ( 10.89 KB , 384x256 , 1647295247452.png )

Image
The brain dump of a managerie of often contracdictory and fantastic characters, from a mossad member turned professor, to an indian student turned part time street criminal, a homeless poet in philly, an ex-skin head, a former chinese investment banker, the owner of a multimillion dollar startup, two military veterans, and a cast of other folks.

The result of a a simple discovery that a rich set of algebraic identities underlay the product of large primes.

For example, a series of variables related to n unknown value, d4, easily derived from known variables. d4a, d4u, d4z, d4H, etc.
Or c/d4, neither known by themselves, but the ratio of which is easily found, and the *product* of which, cd4, is the ratio of our product's factors, b/a.

Be warned the code is thick, mostly written by an indian guy, and translated by the rest of us over time. It is 198 pages of dense work, defining so many variables the original authors resorted to greek letters and elements from the periodic table for naming.

Good luck.

https://pastebin.com/Ad46Awcp
>>
>list of primes
>Miller test (deterministic Miller-Rabin)
>a bunch of tests and variables I have no reason to memorize, followed by a schizo string list including elements of the periodic table..
>...
I'd like to converse with you about your findings, but I have no paper, no lemmas and no theorems.
>>
"I'd like to converse with you about your findings, but I have no paper, no lemmas and no theorems. "

It appears there is no lemma or theorem and OP is working at best in 8th grade math.

I'm just a lowly front end developer but I'll take a crack at it.

All the key variables are on two apparent slopes, d4, and 'gold'.

I could see getting estimates on unknown factors by simply taking known products with common variables (d4a, d4u, d4z, d4H, etc), and then finding the smallest of the set, on the basis that would tightly constrain the upper and lower bound of the factors of those variables, but that assumes none of the unknown factors like u, d4, c, z, or H, are less than 1.

I also wrote some functions for converting variables between the gold line, the _gold line, and the d4 line.

def _gold_to_d4(var):
return (var**2)*Hc


def d4_to__gold(var):
return (var/Hc).sqrt()

This followed from the observation that

d4u/Hc/G
Decimal('1.332259228726959456136695243518871469428174838645985237776420270')

(_goldc**2/G)
Decimal('1.332259228726959456136695243518871469428174838645985237776420270')

meaning

d4u/Hc == some _gold**2

Shows d4 and _gold series are interconnected.

Taking any _gold, multiplying by the power of 2, and then multiplying by Hc should give
a d4 variable.

While dividing a d4 variable by Hc, and then rooting it, should give a _gold series variable

if d4 is inverted (less than 1), then (((_goldc)**2)*Hc)/d4 > d4u
To test this I ran the script a few times with various factors.

a
Decimal('84780491')
b
Decimal('7022437637')

>>> (((_goldc)**2)*Hc)/d4 > d4u
True

>>> a
Decimal('247017779')
>>> b
Decimal('8552169383')

(((_goldc)**2)*Hc)/d4 > d4u
True

>>> a
Decimal('39420301')
>>> b
Decimal('8490912585839')
(((_goldc)**2)*Hc)/d4 > d4u
False

As you can see with the above example, both c and d4 are normal, but c is less than d4,
so _cd4 is less than 1, and the outcome is of course false in this case.
Good to know.

And what this means is that when _cd4 is less than 1, Post too long. Click here to view the full text.
>>
Making sense of all of what looks like dead code, the elListI() function appears to list elements based on some mysterious "I" value.
The elements that have matching I values tend to share lots of identities. Thats all I could tell about it though.
>>
i love primes bro