Menu
▾
▴
Re: [Maxima-discuss] continued fraction display improvement
|
From: Barton W. <wi...@un...> - 2022-08-27 20:31:13
|
>I have hundreds of cf programs for Macsyma that I could share Have you considered putting them on GitHub? --Barton Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows From: Stavros Macrakis<mailto:mac...@al...> Sent: Friday, August 26, 2022 5:03 PM To: Richard Fateman<mailto:fa...@gm...> Cc: <max...@li...><mailto:max...@li...> Subject: Re: [Maxima-discuss] continued fraction display improvement Non-NU Email In the current version of Maxima (5.46, April 2022), you don't need to write ?mplus -- just "+" works. By the way, using "."["."["."]] as last makes for a cute display, but when you reduce the cf using expand or ratsimp, it shows up as an ugly constant name.... On Fri, Aug 26, 2022 at 3:31 PM Richard Fateman <fa...@gm...<mailto:fa...@gm...>> wrote: It seems that a slight change to makecf helps, then load(simplifying); makecf(nums,dens,%,[last]):=block([],negsumdispflag:false, last:if last=[] then "."["."["."]] else first(last), for i:1 thru % do last:simpfunmake(?mplus,[apply(nums,[%-i+1])/last,apply(dens,[%-i])]), last); makeit(nn):= makecf(lambda([n],x-x0),lambda([n],(a[n])),nn); makeit(4); will produce this display a[0]+(-x0+x)/(a[1]+(-x0+x)/(a[2]+(-x0+x)/(a[3]+(-x0+x)/.[.["."]]))) or (if your mail reader shows it as an image...) [cid:ii_l7avcbpw1] On Fri, Aug 26, 2022 at 11:02 AM Stavros Macrakis <mac...@al...<mailto:mac...@al...>> wrote: You aren't really talking about continued fractions in particular, but about generating unsimplified expressions. You can do that using the simplifying package: load(simplifying)$ ex: simpfunmake("+",[2,3]); 3+2 expand(%). << resimplify 5 On Fri, Aug 26, 2022 at 1:56 PM Henry Baker <hb...@pi...<mailto:hb...@pi...>> wrote: Bill Gosper, who no doubt is familiar to most of you, also has put a lot of work over the years into continued fractions (CF's), and consulted for many years with the Macsyma company to improve CF handling & display. He seems to currently be working mostly in Mathematica, but may have some older Maxima/Macsyma code that he might be willing to share. The latest email I have for Bill is <bil...@gm...<mailto:bil...@gm...>>. -----Original Message----- From: Jerry <wj...@gm...<mailto:wj...@gm...>> Sent: Aug 26, 2022 1:14 AM To: <max...@li...<mailto:max...@li...>> Subject: [Maxima-discuss] continued fraction display improvement I have worked on continued fractions for over 40 years and have many unpublished works. It would be of great help to put an option in funmake so that it does not simplify the representation into a mess. The Macsyma company that is now defunct did this so that continued fractions could be easily displayed in the traditional form. [cid:image004.png@01D8BA29.38FB04F0] They produced funmake_no_simp which does not simplify so that one can take numerators and denominators to produce a continued fraction rather than a strictly Regular continued fraction. In this simple example for e^x, the numerators are the same, but suppose I have a symbolic list such as [b[0], b[1], b[2],...] that I want to put in the numerators, I cannot seem to display it. eg makecf(nums,dens,%,[last]):=block([],negsumdispflag:false, last:if last=[] then "."["."["."]] else first(last), for i:1 thru % do last:funmake_no_simp("+",[apply(nums,[%-i+1])/last,apply(dens,[%-i])]), last); then makeit(nn):= makecf(lambda([n],x-x0),lambda([n],(a[n])),nn) Where all the numerators are x-x0 for an expansion around x0 and a[n] are the denominators with a[0] the separate additive term. Maxima produces for e^x a rather non-satisfying representation: --------------------------------------- + 1 x ----------------------------------- + 1 x ------------------------------- - 2 x --------------------------- - 3 x ----------------------- + 2 x ------------------- + 5 x --------------- - 2 x ----------- - 7 x ------- + 2 x --- + 9 . . . It is a great deal of difficulty to take arrays or lists and put them into an equation editor than to just paste the result as previously. To pique your interest, I attach one of my whiz bang results (unpublished). I have hundreds of cf programs for Macsyma that I could share. The bottom left corner is what I need in Maxima. Thank you. WJ Lentz mar...@gm...<mailto:mar...@gm...> 831 601 2120 _______________________________________________ Maxima-discuss mailing list Max...@li...<mailto:Max...@li...> https://lists.sourceforge.net/lists/listinfo/maxima-discuss<https://urldefense.com/v3/__https:/lists.sourceforge.net/lists/listinfo/maxima-discuss__;!!PvXuogZ4sRB2p-tU!BKxW7xnYbROymuFK9MKGB20Bnt3Jba1UuTmV1x6I710dzDg2bAsnS74xK9Z9Q_xw5eBGKn_PMnI4-LdSXQm4$> _______________________________________________ Maxima-discuss mailing list Max...@li...<mailto:Max...@li...> https://lists.sourceforge.net/lists/listinfo/maxima-discuss<https://urldefense.com/v3/__https:/lists.sourceforge.net/lists/listinfo/maxima-discuss__;!!PvXuogZ4sRB2p-tU!BKxW7xnYbROymuFK9MKGB20Bnt3Jba1UuTmV1x6I710dzDg2bAsnS74xK9Z9Q_xw5eBGKn_PMnI4-LdSXQm4$> |