SageMath Puns

Let's look back at some memorable moments and interesting insights from last year.

Your top 10 posts:

• "Intro to Graph Theory with SageMath" by u/MCCPyTutorials
• "Intro to SageMath for Non-Programmers" by u/MCCPyTutorials
• "Mathematics behind cryptography l homomorphism and modulo group applicat..." by u/Siva_915544
• "Happy Cakeday, r/sagemath! Today you're 9" by u/AutoModerator
• "Sorry to bother you all with this question, but I broke typeset output by setting the renderer to SVG and now I can't figure out how to revert it." by u/SirMiba
• "macOS Big Sur Compatability and Download" by u/facearo
• "Writing a proof in Sage Math" by u/jssmith42
• "OEIS Sequence with SageMath (YouTube Short)" by u/MCCPyTutorials
• "How to take a generator of a polynomial ring in Sage? and specify an arbitrary ring polynomial?" by u/GPineda17
• "How do I find the zeros of a function in Sage Math?" by u/hwoodice

Hello All, I am quite new to sage math and working on verifying an answer of mine using the Chinese Remainder theorem.

https://preview.redd.it/1zeblazoume71.png?width=242&format=png&auto=webp&s=167d08b769ef5c2fba1cb8eabc624c44400edae2

Could anyone help me structure the equation for it to work in SageMath?

If I have a list [1,5,4,5,22,0,0,1]

And I need to generate a 7th degree polynomial from this, how should I proceed?

I tried using the methods given online, but it either crashes or gives me a clearly wrong answer ( I'm trying to find all the real roots by using real_roots() or roots() )

Any help would be appreciated.

i want to install sagemath in jupyterlab, I am running Pop_OS 20.10

I was playing with curve25519 in sage earlier today, and noticed that, the order of the elliptic curve group for 25519, is actually larger than the prime field over which it's defined:

ec = EllipticCurve(GF(2**255-19), [0,486662,0,1,0])
ec.order()     
57896044618658097711785492504343953926856930875039260848015607506283634007912


((2^255)-19)

57896044618658097711785492504343953926634992332820282019728792003956564819949

So the ec group contains about 221938542218978828286815502327069187963 more elements.

I understand that the ec group and the field are two different things:

ec group: pairs of points that satisfy the curve equation.

field: Galois field over which we define the curve mod p for certain operations.

I thought the field would put bounds on the curve order too, did I just make a silly assumption?

If I define the curve over a small field, like 17, the number of points is actually only 16.

I need to use Sage for the number theory course that I'm taking this semester, and I've been finding it difficult to use. I'm in my third week this semester and I've been unable to figure out how to use it for anything other than typing the commands from the textbook directly into sagecell.sagemath.org. Is there an easy guide to Sage that would help me use this software? I want to be able to a) figure out what the commands that I'm typing mean, b) figure out how to do things that aren't explicitly laid out for me by my textbook and professor (for my more advanced homework assignments) and c) have a solid foundation in Sage not just for this course but for if I need to use it in subsequent courses or when I go to grad school in a few years.

Hi,

I want to use the SageMath library on repl. (online IDE). Specifically I want to use the factor() function.

Anyone know how to import this successfully?

In my university we are learning to use mathematica for calculations such as flow, curvilinear integrals and limits with multiple variables. Personally I prefer to use open source software, so I was trying sagemath. Now for example, an integral of a function extended to a region of space, on mathematica I can perform it with 2 simple operations. Is this possible with sage? I'm reading the documentation, but I found only integral to one dimension. Here is an example:

https://preview.redd.it/r22epqdfq3m51.png?width=577&format=png&auto=webp&s=09030ae48d3f6cf88eafd68d19f6eb31a63d5114

Hello everyone!

During the course of my physics degree, I've seen a number of my friends use mathematica to easily solve equations. However, I can't afford a license for it and frankly, I wish to support open-source software and hence I'm planning on using sagemath.

I've been using python mostly (scipy, numpy, sympy, matplotlib) to solve equations and especially sympy for CAS. I was wondering if there is any difference in using sagemath as opposed to something like sympy. Also, is it as easy to solve equations in sagemath as it is in mathematica.

Thanks in advance! :)

What do you guys think will be the future of CAS software? Open source like SageMath or closed like Mathematica? I am asking because I am still undecided on which one I should learn to benefit me 10 - 20 years from now. If you could share with me the reasons you choose Sage over Mathematica I would appreciate it. Thanks.

I'm trying to solve a relatively straightforward problem. I would like to figure out when p(b, r) = (b / (r + b)) * ((b - 1) / (r + b - 1)) = 1/2, given the constraints that r and b are integers, and r + b is greater than some minimum value, which we can call m.

I started by solving for b in terms of r with solve([p(b, r) == 1/2], b). That gives me b = r + 1/2*sqrt(8*r**2 + 1) + 1/2, so I make a function blue(r) = r + 1/2*sqrt(8*r**2 + 1) + 1/2. So far so good. The troubles start now when I try to solve for when r + blue(r) >= m.

If I try solve([r + blue(r) >= 10**12], r), it gives me

[[4*r + sqrt(8*r^2 + 1) - 1999999999999 == 0],
 [4*r + sqrt(8*r^2 + 1) - 1999999999999 > 0]]

I can try to solve for r in the first one, but Sage seems to not want to do it. Given solve([4*r + sqrt(8*r^2 + 1) - 1999999999999 == 0], r), it just returns [r == -1/4*sqrt(8*r^2 + 1) + 1999999999999/4]. It won't isolate the r. That's the first problem.

If I take it to Wolfram|Alpha, it can solve it, and I get a value of 1/2 * (1999999999999 - sqrt(1999999999998000000000001)). I can bring that back to Sage. For my purposes, I did:

lowerBound = floor(1/2 * (1999999999999 - sqrt(1999999999998000000000001))) + 1

For convenience's sake, lowerBound is 292893218814. I would now like to solve for the first value of r that is greater than or equal to 292893218814 where p(b, r) = 1/2 and b and r are integers.

I have tried all sorts of things like assume(r > 292893218814), assume(r, integer), assume(b, integer), etc, but it will always either not work or tell me my assumptions are redundant, when they can't be redundant or it would give me the answer I'm looking for without them.

I think I'm just misunderstanding some fundamental things about how SageMath works. If anyone has any tips, I'd greatly appreciate it.

P.S. I'd rather not have any help with the problem itself. I just want to be better at using SageMath effectively.

Thanks guys :)

Hi everyone,

I noticed that SageMath 9.2 has been released for MacOS, but Windows is still on SageMath 9.1. Does anyone know when SageMath 9.2 will be released for Windows? Or has it been released and I'm just missing something?

Thanks!

I've been trying to compile the sagemath package. It failed on gcc 6, so I downloaded gcc-5 from AUR, and tried to compile with that. Another fail. Has anybody successfully managed to compile sagemath on Arch?

Does Sagemath work with OSX High Sierra 10.13.2, which version of Sagemath?

I'm using openSUSE Tumbleweed and I hope this Nix question is appropriate for here because I have SageMath 8.4 installed with Nix and I've been informed that the old SageMath notebook launched by sage -notebook (or just notebook() within a SageMath shell) is going to be removed so I've been trying to run SageMath in a Jupyter notebook. I thought, based on my past hiccups with Matplotlib that I should install Jupyter like this:

nix-shell -p gobjectIntrospection gtk3 'python27.withPackages(ps : with ps; [ matplotlib sage jupyter ])'

(in the hope that doing so would cause Matplotlib to also be available) that didn't work, with this error: https://gist.github.com/fusion809/503c3ad194b164b552056833af4d835d (can't post it here, the ` in it confuses the formatting, making what comes after it not formatted as code). So, therefore, I decided to install it with nix-env -i python2.7-jupyter instead, which ran correctly at the command-line, but when I launched Jupyter through SageMath, with the command sage -n jupyter SageMath launched Jupyter, which was opened in my web browser, when I selected a new SageMath notebook I saw a kernel error message: https://gist.github.com/fusion809/3885dd41b16c5fe8500a2d3ad80fcfc7, which seems oddly reminiscent of the above error message I received from the nix-shell command I tried earlier. I have been told that installing Python modules this way isn't ideal, but since I knew no better way than the ways I have tried I was kind of stuck.

Then, I read the Python segment of the Nixpkgs manual, https://nixos.org/nixpkgs/manual/#python, and inspired by that I created a file, build.nix, with the contents:

with import <nixpkgs> {};

python27.withPackages (ps: with ps; [ sage jupyter ])

(you'll note I omitted the matplotlib module, because at this stage I'll just be happy if I do not have a kernel error message in it) and tried to build it with nix-env -if build.nix but I just received the same error message as above, https://gist.github.com/fusion809/503c3ad194b164b552056833af4d835d.

Any ideas what I should do? Thanks in advance for using some of your precious spare time to help me. :)

Let's look back at some memorable moments and interesting insights from last year.

Your top 10 posts:

• "Whoever wrote the tutorials for sagemath.org is a cool dude" by u/BoomerFelonOwl
• "Video I Made on SageMath Help Features and Other Resources - May Be of Interest Here" by u/MCCPyTutorials
• "My SageMath setup, using vim-Jupyter plugin for interaction." by u/noscore6
• "SageMath vs Mathematica" by u/alexice89
• "Making examples of practical undergrad uses of sagemath?" by u/BoomerFelonOwl
• "Beginner to sagemath" by u/kn0xchad
• "Help with a probability problem SageMath seems to not want to simplify" by u/PM_ME_YOUR_FUN_MATH
• "truncating tensor/scalar perturbations at particular orders...(general relativity)" by u/deepdooper
• "Mathematical software: sagemath and mathematica" by u/TeoCol777
• "Change colors on docker"

In my university we are learning to use mathematica for calculations such as flow, curvilinear integrals and limits with multiple variables. Personally I prefer to use open source software, so I was trying sagemath. Now for example, an integral of a function extended to a region of space, on mathematica I can perform it with 2 simple operations. Is this possible with sage? I'm reading the documentation, but I found only integral to one dimension. Here is an example:

https://preview.redd.it/r22epqdfq3m51.png?width=577&format=png&auto=webp&s=09030ae48d3f6cf88eafd68d19f6eb31a63d5114