PyMath #1 - When Dividorial and π Meet - Their Beautiful Relationship

[@NPCooking69]

Hey you! Thanks for watching ! :D If you like what you saw, please share the video around :3 Love ya

[@integralboi2900]

You forgot to add the program in the description.

[@PapaFlammy69]

@@integralboi2900 Damn, thank you! It's been fixed now! =)

[@MidnightStorm4990]

@@PapaFlammy69 I'm 16 and I'm really interested in math and learning python. Do you think I should attempt to learn this if I have barely any knowledge on python.

[@Ou_dembele]

@@MidnightStorm4990 yes start learning

[@AndrewDotsonvideos]

Hey last year you left your used napkin on my kitchen table. I still have it in case you want it back.

[@PapaFlammy69]

yes, plis.

[@Hi_Brien]

Ahahaha

[@nanare137]

Introducing python coding and combining it with the maths on the board sounds good, I would love more of that. :)

[@dr.peyamsfan7158]

Mee too

[@PapaFlammy69]

Definitely more to come soon! :) Thank you for your feedback! :)

[@osman9750]

Exactly my thoughts 👍

[@JB-ym4up]

Python and n¡? Sounds like a skit about mathematicians that say n¡.

[@VerrouSuo]

Dividorial is one of the silliest sounding terms I’ve ever heard of... but it’s starting to grow on me.

[@PapaFlammy69]

Same! :DDD

[@trevorallen3212]

Just looking at the word has me going crazy pronouncing it.

[@livedandletdie]

I'm divided on the subject, but it's a factor, and that can I double down on.

[@aasyjepale5210]

Woah that escalated quickly... Or did it?

[@trevorallen3212]

@@aasyjepale5210 exponentially or linearly?

[@demetriuspsf]

YES! More dividorial vídeos. I still can't wrap my head around how awesome the notation for the dividorial is. Every time I think of the upside down exclamation I laugh.

[@paulthiede]

same

[@sergey1519]

¿Qué quieres decir con eso? ¡Eso no es divertido!

[@rupachakraborty3681]

Papa Flammy's Math video with Python programming.........it's worth more than gold. #pymath

[@integralboi2900]

Python? It thought this was a maths channel. Just kidding, these videos are interesting.

[@PapaFlammy69]

=D

[@fractal_lynn]

Hey! I'd recommend using the Decimal class for longer decimals! Also python can handle converting to floating point automatically. the range in the for loop can specify an iteration, so "for i in range(x,1,-2): temp *= i" would solve the double factorial You can also slip expressions directly into the return statement! Some of the code that I wrote; def factorial( x, degree = 1 ): temp = 1 for i in range( x, 1, -degree ): temp *= i return temp def dividorial( x ): return float( factorial( x, 2 ) / factorial( x - 1, 2 ) ) You might be interested in Spyder, software included in the Anaconda package that comes with interesting packages built in for python data science. I am very impressed by how this converges to pi/2, this is incredible! I really enjoy your videos and look forward to more videos like this! :D

[@PapaFlammy69]

Thank you for the python bits and the information!!! :)

[@fractal_lynn]

@@PapaFlammy69 Oh, and I apologize. I was mistaken for thinking pi/2 was Tau. Tau is actually 2pi, but I'm sure you already know that haha. Anyways, you might be interested to know Python has a built in complex number datatype. I used it to generate the fractal in my pfp lol.

[@paulthiede]

@@fractal_lynn i thought that was just a zoomed in screenshot of the mandelbrot set. but yeah, i dont think the mandelbrot set looks like that, at least not exactly.

[@fractal_lynn]

@@paulthiede Its a warped version of the mandelbrot set zoomed in to a specific spot; i.imgur.com/LmK5kow.png Also heres strictly the mandelbrot set with 1024 iterations per pixel; (7200x5400) i.imgur.com/l7JFxrH.png

[@SugarBeetMC]

@@fractal_lynn There is math.factorial().

[@jfbb3369]

I really liked the python part. I’m learning programming myself and this really helps understand math concepts in programming. Thanks Jens 😉

[@PapaFlammy69]

Nice :3

[@dozzco2827]

This video was incredibly interesting and incredibly fun to watch! Looking forward to more of this soon!

[@PapaFlammy69]

Thank you! =)

[@MYCHEM1CALR0M4NCE]

I've followed you for a few years now and this is one of my favourite videos you've made :))

[@eliyasne9695]

Seeing your video, one might mistakenly think that the limit notation was invented in the middle east.

[@PapaFlammy69]

:'D

[@shoopinc]

Al-Cauchy

[@TheMahu13]

Since an integer is either odd or even, you can simply right : if x%2 == 0 : ... else : ... Your video was really cool I like the addition of programming in it :) ! Keep it up

[@sergey1519]

you can also do just: if x%2: ... else: ... but much cleaner solution was suggested for double factorial, that avoids cases altogether: df = 1 for i in range(x, 1, -2): df *= i return df

[@angelmendez-rivera351]

Also, here are some relatively simple formulae that extend the dividorial to complex arguments via analytic continuation. z¡ = 1/sqrt(2)·(z/2)!/[(z - 1)/2!]·sqrt(π/2)^[cos(πz)/2] z¡ = 1/sqrt(2π)·2^z·sqrt(π/2)^[cos(πz)/4]·(z/2)!^2/z! z¡ = 1/sqrt(2π)·exp[ln(2)·z + ln(π/2)/4·cos(πz)]·Π(z/2)^2/Π(z) The second formula removes the quotient of difference from the first, just in case that this is preferrable, and expresses the dividorial more directly in terms of the factorial, and the third formula addresses notational concerns: namely, it uses the Π notation for complex numbers instead of using the factorial notation, which is, in most rigorous circles, reserved for natural numbers only, and it also addresses the concern that complex exponentiation is technically multivalued, and that the analytic exponential function should be used instead.

[@SugarBeetMC]

0:44 On the one hand, Pi. On the other hand, Py.

[@m1kr0kosmos]

thanks for the python stuff. Since you, Lex Fridman & Andrew Dotson & 3B1B have some shared interest, I am surely not the only one hoping you all bro it out

[@budtastic1224]

I've decided to call n!/(n-1)! the unfactorial

[@MathIguess]

Bruh my math degree is killing me and then I watch more math to unwind What am I even? Lol

[@PapaFlammy69]

bruhv

[@NathanThomasisepic]

You didn't have to split it up into odd and even cases. You can just do: def dfac(x): df = 1 for i in range(x, 1, -2): df *= i return df This loops through the range from x to 1 and subtracts 2 for each iteration.

[@PapaFlammy69]

That is nice, I didn't know about that, thank you! :)

[@VACatholic]

@@PapaFlammy69 No papa. Just do: def dfac(x): return prod(range(x,1,-2)) "prod" is the product of all elements in a list, no need to write loops :(

[@VACatholic]

no need to loop. Use "prod"

[@nitroh7745]

I LOVE this new format it is amazing!! Thanks papa

[@fedebonons8453]

The CIAO at the end is poetry for my italian ears

[@MathIguess]

I like the math + python format, it's nice Also, well made video, damn :D

[@AggieRinse]

Love the python/coding bits! More!

[@mastershooter64]

papa flammy are you gonna invent fractional rank tensors and fractional dimensional calculus in the next upload?

[@bmw123ck]

Love python being introduced to the channel!!

[@Kdd160]

Lol 4:05 what is wafu do u mean wife?? 🤣🤣🤣🤣😂😂 I think dividorial is pi's boyfriend and pi is dividorial's girlfriend what a beautiful relationship 😍😻💟 And u pronounce 5 like FAF just like y as whayy(inspiriert von Ozarks bester Netflix-Serie) And u pronounce pi over 2 like paaaye over two ("""")

[@emonaetothexlisa9194]

Paaayeee over 2 lmao

[@paulthiede]

lmao

[@davids841]

Using modulus two inside some of the loops can make for some much faster and nicer branchless coding. This would help if you were to try very large values of n.

[@matrefeytontias]

This `elif` with the entire test of whether x is odd after we already know it's not even hurts my little programmer heart. At least the cool math makes up for it ;_;

[@paulthiede]

yeah, he could just have made an else. I noticed that as well and I dont understand it too.

[@gabriel7233]

He didn't even have to consider 2 cases and just made loop with step 2

[@matrefeytontias]

Yeah, but for someone who's not a programmer and only follows the mathematical logic, it makes sense to have a disjunction of cases. Otherwise, just loop from n down to 1 exclusive in steps of 2 and multiply throughout.

[@Daniel-nl3ug]

This code makes sense because if it's not even, it could be a decimal, in which case it'd make sense to have an else statement and then raise an input error exception, papa just didn't include that part because it's not very relevant to the video.

[@paulthiede]

@@Daniel-nl3ug you're right! i didnt think of that

[@paulthiede]

def double_fac(x): product = 1 for i in range(x,1,-2): product = product * i return product This would also work as a double factorial

[@VACatholic]

so would: def double_fack(x): return prod(range(x,1,-2))

[@zh84]

Or just import scipy, which has a "factorial2" function that calculates the double factorial for you.

[@PapaFlammy69]

lel

[@dhoyt902]

Hi Papa FM I have a question. You know how we accelerate summation series with partial sums(like Euler transform(ie van vigngaarden) and Shanks? Are there methods to accelerate this product? Or, would you have to treat the product as a series of repeated sums and then apply acceleration?

[@angelmendez-rivera351]

Most likely the latter.

[@angelmendez-rivera351]

What I find neat about this is the fact that you unknowingly found a zeroth order asymptotic expansion for the dividorial. lim (2n)¡^2/(2n + 1) = π/2 implies lim (2n)¡/sqrt(2n + 1) = sqrt(π)/sqrt(2). sqrt(2n + 1) = sqrt(2)·sqrt(n + 1/2), so lim (2n)¡/sqrt(n + 1/2) = sqrt(π). Incidentally, this implies that lim (2n)¡/sqrt(n + 1/2) = (-1/2)!, so there is another cool connection between the dividorial and the factorial. Anyhow, the point I am making is that (2n)¡ ~ sqrt(π)·sqrt(n + 1/2), or equivalently, (2n)¡ = O(sqrt(n + 1/2)). This is useful because it gives you an understanding of growth that is now useful for calculating many other limits.

[@iridium8562]

Im seriously getting pissed at youtube, i got this recommended but I wasn’t notified..? What..?

[@PapaFlammy69]

I hate it so god damn much :(((

[@jack_papel]

12:00 I might be wrong but I'm pretty sure implications don't work like that

[@mkdk6831]

Is there a way to write the dividorial as an integral formula just like the gamma- or pi function for the factorial?

[@PapaFlammy69]

haven't found one yet sadly :3

[@aaronneugebauer4752]

There exists an integral representation of 1/dividorial through n!!/(n-1)!!= Integral from 0 to pi/2 of sin(x)^n dx (*2/pi if n is even). So you can conclude that 1/divdorial=Integral from 0 to pi/2 sin(x)^n/((3*pi/4)+pi/4*(-1)^(n+1)) dx. WolframAlpha then gives an even better formular for dividorial with some transformations: www.wolframalpha.com/input/?i=1%2F%28%283*pi%2F4%29%2Bpi%2F4*%28-1%29%5E%28n%2B1%29%29+sin%28x%29%5En+dx+from+0+to+pi%2F2

[@angelmendez-rivera351]

Arjun Vyavaharkar I am not a fan of the fact that you defined a shifted version instead of simply defining an extension of the original. It is needlessly cumbersome.

[@angelmendez-rivera351]

Arjun Vyavaharkar *Use it or lose it* Nobody is going to be using such an arbitrary formula of an arbitrarily named function that doesn't really show up in applications. Also, I guess you're too damn immature to accept constructive criticism. It's not my fault you're so hopelessly arrogant.

[@angelmendez-rivera351]

Arjun Vyavaharkar *I don't care about what you think* And I don't care that you don't care. Grow up, kid.

[@diegoalbertoduransalas4806]

that was so... clean!!

[@gabriel7233]

Why didn't you just made a loop from n to 1 with step -2 in dfac(x), it would be just df *= i and you wouldn't have to consider 2 cases

[@Arnodorian6125]

Pls do a live Q&A

[@PapaFlammy69]

On my list already! :)

[@PositronQ]

In the seco d Python program you can implement the code into a Github repository and I like This implementation

[@hugh1643]

this was stupid why not just print(math.pi/2) smh my head

[@PapaFlammy69]

;_;

[@paulthiede]

xD

[@tomkerruish2982]

"smh my head"? :D

[@ElZedLoL]

Shaking my head my head

[@ycombinator765]

@@tomkerruish2982 shaking (my head) ^2

[@hugh1643]

so 1/dividorial looks a bit like the integral of sin^n(x) or cos^n(x) from 0 to pi/2 just looking at the wikipedia page for double factorials.

[@maksimpobedinsky5458]

nice

[@victorrizkallah6014]

The python part is great

[@hhtt37]

For python, maybe jupyter could be better for video, mixing live execution and cached results. Math plots could be enlightening!

[@diegowielandt7089]

great format

[@comradeelmo5272]

I'm so used to visual studio it felt weird not seeing visual studio

[@PapaFlammy69]

xD

[@enzoqueijao]

He back ¡

[@elijahcamcam4611]

where do you buy a chalkboard

[@torment808]

But what are their applications?

[@nnniv]

Ayy python

[@juansalazar6849]

Math goes wild

[@OriginalSuschi]

Ok so if we want to write the double factorial in terms of factorial, we would get this, right? : n!! = (2*(n/2))(2*(n/2 -1))...*(2*2)(2*1), which then equals to (n/2)!*2^(n/2) For even n's, right?

[@paulthiede]

though im not sure I think you made an error with "2*(n/2-1)" I think it would need to be "2(n/2)-1" in order to be correct

[@OriginalSuschi]

Clashgamer04 if it was according to my assumption, for n=6 it would be: (2*(6/2))(2*(6/2 -1))(2*(6/2 -2))=(2*3)(2*2)(2*1), which is the thing PapaFlammy wrote as well. With your assumption it would be: (2*(6/2))(2*(6/2)-1)(2*(6/2)-2 = 2*3*2*5*2*4 PapaFlammy already wrote (2*(n/2 - 1)) out to (2*n/2 -2)

[@angelmendez-rivera351]

Yes, n!! = 2^(n/2)·(n/2)! for even n is correct indeed.

[@Milldyria]

Is (ai)^2+(bi)^2=(ci)^2 Or any sutch identity?

[@paulthiede]

i mean (ai)² + (bi)² = aiai + bibi = -a² - b² and (ci)² = cici = -c² so -a² - b² = -c² which would be a² + b² = c² | *(-1). So just the pythagorean formular multiplied by -1. But then again its obvious that that would be true if you do the following: a² + b² = c² | *i² a²i² + b²i² = c²i² = (ai)² + (bi)² = (ci)². But because i²= -1 we are looping around with the same expressions. To answer your question, idk if there would be any use to your identity, but I think it would be a funny way of calculating the length of a hypothenuse in the 9th grade. :-)

[@___c5645]

Clashgamer04 Dividorial notation, not i. a¡ is different from ai

[@Milldyria]

@@paulthiede I was thinking about the dividorial. It would be very interesting to see if there are any identities with it. Like I dunno, pascals triangle but with dividorials. Or golden ratio in dividorials terms.

[@tablonhobbs45]

Still love ya, papa

[@OriginalSuschi]

Wait but which dividorial definition do you mean here? 🚰 Asking before watching the video

[@angelmendez-rivera351]

Why would you ask before the watching the video? That is so pointless.

[@OriginalSuschi]

Angel Mendez-Rivera yeah I know, but I like pointless things. Life is pointless as well... And I am 14 and this is deep

[@diederickfloor4261]

Poopietorial

[@rokarus7658]

So since we have lim ((2n)¡)^2/(2n+1) -> pi/2, does that mean that we have n¡ ~ n^1/2 ?

[@angelmendez-rivera351]

I think it would be more accurate to say n¡ = O(n^(1/2))

[@Milldyria]

What is the plot of Sin(xi) or cos(xi)?

[@hughmungous1539]

Sinh(x) and cosh(x) respectively I'm pretty sure

[@PapaFlammy69]

kinda, yes

[@angelmendez-rivera351]

I think he meant the plot of sin(x¡) and cos(x¡). But honestly, there are no "good" continuous representations of the double factorial... let alone the dividorial. And no, I'm not counting that one strange continuous representation that Wikipedia uses for the double factorial. It hardly makes any sense to use, given the context.

[@Supremebubble]

I noticed

[@PapaFlammy69]

[@OmegaBlackjack]

Ok, why didn’t we just multiplied wa*2 to seek for π, not π/2 ? xD

[@matthewcapstick6242]

At 7:48 is that a sin or a k or a h?

[@angelmendez-rivera351]

Yes

[@phat5340]

Papa flammy do you have depression ?

[@PapaFlammy69]

wat? Hell no lol

[@phat5340]

You look like you're having a Vietnam flashback on almost every video

[@PapaFlammy69]

xD

[@klofat]

do not mix up П and п.

[@djvalentedochp]

Nice

[@Frankie13074]

33th!!!

[@PapaFlammy69]

nice xD

[@matron9936]

But what is !i¡!?

[@michaelempeigne3519]

pi^2 does not equal g

[@angelmendez-rivera351]

It does, though

[@michaelempeigne3519]

@@angelmendez-rivera351 only to one decimal place

[@angelmendez-rivera351]

Michael Empeigne Yeah. That's good enough.

[@spaghetti1383]

Nice Jojo reference

[@sovietcat919]

Papa answer this question I bet you can’t What’s 2+2

[@PapaFlammy69]

4 apples

[@sovietcat919]

@@PapaFlammy69 oh my I underestimated your intelligence

[@PapaFlammy69]

indeed my dear son.

[@sovietcat919]

@@PapaFlammy69 Dear Papa, I am sorry to inform you that the answer is actually 22 apples

[@ricardoparada5375]

I never liked python syntax (c++/java looks so much better lmao) but this is still some really cool stuff

[@tryphonunzouave8384]

4 views 4 thumbs up, talk about ratio

[@shlokmahajan8439]

Hello bro love from India,your vedio helps a lot

[@OonHan]

haha yes

[@footballknight7]

I am second NGL!

[@fujatv503]

*laughs in engineering* i havent understood a word you said, but i know that it is useless

[@paulthiede]

is it really that hard? i mean i always though like "why is he explaining that 6/2 =3?"

[@fujatv503]

@@paulthiede oh that is what he was trying to tell in this video, thanks for translating

[@mattchris3421]

GENIUS likes to use blackboard not White board....using white board u become dumb cuz keep asking ' Why Board? '

[@mudkip_btw]

I think a better name of the dividorial would be fictorial, short for fiction-orial, get it, fact-orial, fiction-orial cuz the dividorial is kind of made up ok imma leave

[@mudkip_btw]

Do like the number pahhh

[@mudkip_btw]

Please stop using the upside down ! :d just use a ?

[@RealLifeKyurem]

Randy ? is already used for the termial, the addition version of the factorial. The function is introduced by Donald Knuth in 1997. Example: 5? = 1 + 2 + 3 + 4 + 5 5? = 15 In general: n? = n(n+1)/2

[@mudkip_btw]

@@RealLifeKyurem ah I didn't know that, ty

[@integralboi2900]

n. is a good option, but so is n, and maybe n;.

[@twiddle7125]

May I make a suggestion to just slow down your speech just a little bit. Actually, you don't need to slow down, but you should take a breath between sentences, lol.

[@VACatholic]

Your code sucks papa. You should just: def dfac(x): return prod(range(x,1,-2)) I'm disappointed.