Maxima Digest, Vol 32, Issue 51



On Fri, Mar 27, 2009 at 6:03 PM, <maxima-request at math.utexas.edu> wrote:

> Send Maxima mailing list submissions to
>        maxima at math.utexas.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
>        http://www.math.utexas.edu/mailman/listinfo/maxima
> or, via email, send a message with subject or body 'help' to
>        maxima-request at math.utexas.edu
>
> You can reach the person managing the list at
>        maxima-owner at math.utexas.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Maxima digest..."
>
>
> Today's Topics:
>
>   1. Re: Problems with complex numbers (Volker van Nek)
>   2. Clozure CL + Windows XP (Barton Willis)
>   3. Re: Problems with complex numbers (Michel Talon)
>   4. Re: Fwd:  3D Plot scale (Mario Rodriguez)
>   5. Compiling GCL from CVS (Dieter Kaiser)
>   6. Re: Compiling GCL from CVS (Raymond Toy)
>   7. Re: Compiling GCL from CVS (Billinghurst, David (RTATECH))
>   8. plotting results of rk (Semyon Glazyrin)
>   9. Re: plotting results of rk (Martin Sch?necker)
>  10. Re: plotting results of rk (Jaime Villate)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 26 Mar 2009 18:57:51 +0100
> From: Volker van Nek <volkervannek at googlemail.com>
> Subject: Re: [Maxima] Problems with complex numbers
> To: Michel Talon <talon at lpthe.jussieu.fr>, Maxima at math.utexas.edu
> Message-ID: <49CBC21F.40206 at googlemail.com>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>
> (%i2) gfactor(%i*k^4 -2*k^3 -2*%i*k^2+2*k+%i);
>
> (%o2) %i*(k-1)*(k+1)*(k+%i)^2
>
> See ?gfactor
>
> HTH
>
> Volker van Nek
>
>
> Michel Talon schrieb:
> > Hello,
> >
> > maxima has difficulties factoring polynomials on the complex field.
> > As an example, compare the following maple and maxima computation:
> >
> > niobe% maple
> >     |\^/|     Maple 9 (IBM INTEL LINUX)
> > ._|\|   |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc.
> > 2003
> >  \  MAPLE  /  All rights reserved. Maple is a trademark of
> >  <____ ____>  Waterloo Maple Inc.
> >       |       Type ? for help.
> >
> >> factor(I*k^4 -2*k^3 -2*I*k^2+2*k+I);
> >>
> >                                       2
> >                (k - 1) (k + 1) (k + I)  I
> >
> >
> >> quit
> >>
> > bytes used=978068, alloc=851812, time=0.05
> >
> >
> > niobe% maxima
> > (%i1) algebraic:true;
> > (%o1)                                true
> > (%i2)  factor(%i*k^4 -2*k^3 -2*%i*k^2+2*k+%i);
> >                                            2
> > (%o2)                 (k - 1) (k + 1) (%i k  - 2 k - %i)
> > (%i3)  factor(%i*k^4 -2*k^3 -2*%i*k^2+2*k+%i,b^2+1);
> >                                                     2
> > (%o3)                    %i (k - 1) (k + 1) (k + %i)
> >
> > Note that it is necessary to specify that the computation is in the field
> of
> > b^2+1 to get the answer. Moreover in many cases, b is returned in place
> > of %i and many simplifications are missed.
> >
> > Is there any way to get more efficient computations?
> >
> > Thanks you very much
> >
> >
>
>
>
> ------------------------------
>
> Message: 2
> Date: Thu, 26 Mar 2009 14:32:46 -0500
> From: Barton Willis <willisb at unk.edu>
> Subject: Clozure CL + Windows XP
> To: "maxima at math.utexas.edu" <maxima at math.utexas.edu>
> Message-ID:
>        <OF89F979C1.0A26D628-ON86257585.006877A2-86257585.006B5F80 at unk.edu>
> Content-Type: text/plain; charset="US-ASCII"
>
> I tried building Maxima using Clozure CL + Windows XP. The tests rtest14,
> rtest15, and rtest_gamma abort; the rest of the tests
> pass. Deleting these three tests, the run times are:
>
> SBCL 1.0.22   109.9 seconds,               6.105 seconds GC   9 850 167
> 496 bytes consed
> Clozure CL      162.0 seconds,            28.063 seconds GC   3 745 442
> 024 bytes allocated
>
> Using Clozure CL, share_testsuite reports some errors, but it finishes:
>
> Errors found in rtest_numericalio.mac, problems: (50 52 57)
> Errors found in testprintf.mac, problems: (27 54)
> Error found in test_abs_integrate.mac, problem: (75)
> Error found test_pochhammer.mac, problem: (29)
>
> I used a 2.13 GHz, 2 GB RAM, Core 2 CPU. One bug that I know is (that
> shows up in rest14 or rtest15).
>
> (%i1) bessel_j(3,2.0);
>
> (%o1)                                 0.0
> (%i2) build_info();
>
> Maxima version: 5.17post
> Maxima build date: 13:56 3/26/2009
> host type: @host@
> lisp-implementation-type: Clozure Common Lisp
> lisp-implementation-version: Version 1.3-dev-r11537M-trunk  (WindowsX8632)
>
> The GNU Autotools build doesn't work, I think (the path names are wrong);
> I used the Lisp only build. I had
> to replace defsystem.lisp with the file that comes with Clozure CL.  Also,
> I think  (maxima-dump) doesn't work.
>
> For information on Clozure CL, see:
>
> http://trac.clozure.com/openmcl
> http://article.gmane.org/gmane.lisp.openmcl.devel/2878
>
> The first 90% of most projects is easy, but it's the final 9.9% that
> really counts. For (open-source, free, ...) CL on Windows we
> have about five ninety percent solutions (GCL, ABCL, SBCL, ECL, CLISP).
>
> Barton
>
> My Maxima CVS isn't 100% clean--I shared some results with a few of you
> yesterday that were somewhat different. I did
> these results  on a different machine.
>
> Barton
>
>
> ------------------------------
>
> Message: 3
> Date: Thu, 26 Mar 2009 21:03:36 +0100
> From: Michel Talon <talon at lpthe.jussieu.fr>
> Subject: Re: [Maxima] Problems with complex numbers
> To: maxima at math.utexas.edu
> Message-ID: <gqgn2p$5qh$1 at ger.gmane.org>
> Content-Type: text/plain; charset=us-ascii
>
> Volker van Nek wrote:
>
> >
> > (%i2) gfactor(%i*k^4 -2*k^3 -2*%i*k^2+2*k+%i);
> >
> > (%o2) %i*(k-1)*(k+1)*(k+%i)^2
> >
> > See ?gfactor
> >
> > HTH
> >
> > Volker van Nek
>
> Thanks, i had missed this one. Now the following problem which was
> in my computation:
>
>
> niobe% maxima
> (%i1) gfactor(sqrt(%i*k^4 -2*k^3 -2*%i*k^2+2*k+%i));
>                                 2               2
> (%o1)                 sqrt((%i (k  - 1) - 2 k) (k  - 1))
> (%i2) radcan(%o1);
>                                                 2
> (%o2)           sqrt(k - 1) sqrt(k + 1) sqrt(%i k  - 2 k - %i)
>
> (%i6) map(gfactor, sqrt(%i*k^4 -2*k^3 -2*%i*k^2+2*k+%i));
>                                                      2
> (%o6)                 sqrt(%i (k - 1) (k + 1) (k + %i) )
> (%i7) radcan(%);
>                                                 2
> (%o7)           sqrt(k - 1) sqrt(k + 1) sqrt(%i k  - 2 k - %i)
>
>
> This misses the opportunity to get (k+%i) sqrt(k^2-1) and to simplify
> (k+%i)
> with other factors in the computation. In fact i had to help a lot of the
> computation by hand to get the simple result at the end.
>
> Of course this is non trivial, and i don't pretend that maple is good for
> manipulations of square roots!
>
>
>
>
> --
> Michel Talon
>
>
>
> ------------------------------
>
> Message: 4
> Date: Thu, 26 Mar 2009 22:28:23 +0100
> From: Mario Rodriguez <biomates at telefonica.net>
> Subject: Re: [Maxima] Fwd:  3D Plot scale
> To: hagen.sam at gmail.com
> Cc: "maxima at math.utexas.edu" <maxima at math.utexas.edu>
> Message-ID: <1238102903.5880.7.camel at pc>
> Content-Type: text/plain
>
> Hi,
>
> > 1)  How should one go about printing a graph in gnuplot?
>
> You can save your graph in a file (png, gif, postscript, ..) and then
> print it.
>
>
>
> > 2)  Is there a way to plot polar coordinates in Maxima?
>
> See these examples made with package draw:
>
> http://www.telefonica.net/web2/biomates/maxima/gpdraw/noncartes
>
>
> See also that some examples make use of the terminal option this way:
>
> terminal = eps
>
> in order to create a postscript file.
>
>
>
> --
> Mario
>
>
>
> ------------------------------
>
> Message: 5
> Date: Thu, 26 Mar 2009 22:53:43 +0100
> From: "Dieter Kaiser" <drdieterkaiser at web.de>
> Subject: Compiling GCL from CVS
> To: "Maxima" <maxima at math.utexas.edu>
> Message-ID: <1CEAF7655E0E41E8A63C40056D0A9273 at Kaiser>
> Content-Type: text/plain;       charset="us-ascii"
>
> I have tried to install all what is needed to compile GCL from CVS. I have
> finished the installation of mingw, minsys, tcltk-8.4.1-1.exe,
> msysDTK-1.0.1.exe, ...
>
> I have followed this advise
>
> "The easiest way is to
> download and install ActivePerl binaries provided by
> <http://www.activestate.com>.";
>
> and have installed ActivePerl. (Question: Do I have to install ActivePerl
> at a
> special place?).
>
> The command ./CONFIGURE seems to works. The command MAKE stops- These are
> the
> last messages:
>
> gcc -D_GNU_SOURCE -c -fsigned-char -pipe -Wall
>  -fno-zero-initialized-in-bss
> -mms-bitfields  -O3  -march=i686 -mfpmath=387 -I/home/Dieter/gcl/gcl/o
> -I../h
> -I../gcl-tk cfun.c
> cfun.c:27:19: dlfcn.h: No such file or directory
> cfun.c: In function `fSdlsym':
> cfun.c:187: warning: implicit declaration of function `dlerror'
> cfun.c:191: warning: implicit declaration of function `dlsym'
> cfun.c:191: warning: assignment makes pointer from integer without a cast
> cfun.c:197: error: `RTLD_DEFAULT' undeclared (first use in this function)
> cfun.c:197: error: (Each undeclared identifier is reported only once
> cfun.c:197: error: for each function it appears in.)
> cfun.c:197: warning: assignment makes pointer from integer without a cast
> cfun.c:198: warning: assignment makes pointer from integer without a cast
> cfun.c: In function `fSdladdr':
> cfun.c:207: error: `Dl_info' undeclared (first use in this function)
> cfun.c:207: error: syntax error before "info"
> cfun.c:212: warning: implicit declaration of function `dladdr'
> cfun.c:212: error: `info' undeclared (first use in this function)
> cfun.c: In function `fSdlopen':
> cfun.c:232: warning: implicit declaration of function `dlopen'
> cfun.c:232: error: `RTLD_LAZY' undeclared (first use in this function)
> cfun.c:232: error: `RTLD_GLOBAL' undeclared (first use in this function)
> cfun.c:232: warning: assignment makes pointer from integer without a cast
> make[1]: *** [cfun.o] Error 1
> rm list.c
> make[1]: Leaving directory `/home/Dieter/gcl/gcl/o'
> make: *** [unixport/saved_pre_gcl] Error 2
>
>
> It seems to me that a header file is missing:
>
>   cfun.c:27:19: dlfcn.h: No such file or directory
>
> Has anyone an idea what I am doing wrong?
>
> Dieter Kaiser
>
>
>
>
> ------------------------------
>
> Message: 6
> Date: Thu, 26 Mar 2009 19:35:57 -0400
> From: Raymond Toy <toy.raymond at gmail.com>
> Subject: Re: [Maxima] Compiling GCL from CVS
> To: Dieter Kaiser <drdieterkaiser at web.de>
> Cc: Maxima <maxima at math.utexas.edu>
> Message-ID: <49CC115D.9080804 at gmail.com>
> Content-Type: text/plain; charset=ISO-8859-1
>
> Dieter Kaiser wrote:
> > I have tried to install all what is needed to compile GCL from CVS. I
> have
> > finished the installation of mingw, minsys, tcltk-8.4.1-1.exe,
> > msysDTK-1.0.1.exe, ...
> >
> I don't have any particular advice on how to do this.  Perhaps the gcl
> mailing lists can help.  But even on Linux, I have a very hard time
> building gcl.  I've basically given up and just use the old version that
> I somehow compiled many years ago.
>
> Maybe ecl would be a good alternative?  I have no experience with that
> on windows.  Works ok on Linux and MacOS.
>
> It certainly would be nice if the free versions of LW or ACL could
> compile maxima.
>
> Ray
>
>
>
> ------------------------------
>
> Message: 7
> Date: Fri, 27 Mar 2009 10:39:16 +1100
> From: RTATECH"
>        <David.Billinghurst at riotinto.com>
> Subject: Re: [Maxima] Compiling GCL from CVS
> To: Dieter Kaiser <drdieterkaiser at web.de>, Maxima
>        <maxima at math.utexas.edu>
> Message-ID:
>        <
> 026DCC31AB859648A6F16C0E5CD2580D018CE8AF at calttsv025.cal.riotinto.org>
> Content-Type: text/plain; charset="utf-8"
>
> > From:  Dieter Kaiser
>
> > I have tried to install all what is needed to compile GCL
> > from CVS. I have finished the installation of mingw, minsys,
> > tcltk-8.4.1-1.exe, msysDTK-1.0.1.exe, ...
>
> Deiter,
>
> This is probably a question for gcl-devel at gnu.org
> http://lists.gnu.org/mailman/listinfo/gcl-devel
>
> I have built CVS gcl quite a few times using msys.  I need to set up
> the build system on a new PC shortly - target is end of Feb :-(
> As I do this I will check the documentation.
>
> The build is fragile.  You need specific versions of gcc and binutils.
>
> You need perl in your path.  I think ActiveState does this by default.
>
>        David
>
> This email is confidential and may also be privileged. If you are not the
> intended recipient, please notify us immediately and delete this message
> from your system without first printing or copying it. Any personal data in
> this email (including any attachments) must be handled in accordance with
> the Rio Tinto Group Data Protection Policy and all applicable data
> protection laws.
>
>
>
>
>
> ------------------------------
>
> Message: 8
> Date: Fri, 27 Mar 2009 08:54:39 +0300
> From: Semyon Glazyrin <siglazyrin at gmail.com>
> Subject: plotting results of rk
> To: maxima at math.utexas.edu
> Message-ID: <20090327055439.GA5492 at voodoo.zuzino.mipt.ru>
> Content-Type: text/plain; charset=us-ascii
>
> Hello!
>
> How can I plot the result of solution by rk of system of diff equations:
> For example
>
> $ load(dynamics);
> $ res:rk([x+y,x-y],[x,y],[1,0],[t,0,1,0.1]);
>
> This returns me a list of [t,x,y]. How can I plot x(t) and y(t)?
>
> thanks,
> Semyon.
>
>
> ------------------------------
>
> Message: 9
> Date: Fri, 27 Mar 2009 08:44:46 +0100
> From: Martin Sch?necker <ms_usenet at gmx.de>
> Subject: Re: [Maxima] plotting results of rk
> To: maxima at math.utexas.edu
> Message-ID: <49CC83EE.7090208 at gmx.de>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>
> > This returns me a list of [t,x,y]. How can I plot x(t) and y(t)?
>
> Separate your values into pairs [t,x] and [t,y] and use the discrete
> option of plot2d:
>
> xtvalues: map(lambda([e], part(e,[1,2])), res);
>
> ytvalues: map(lambda([e], part(e,[1,3])), res);
>
> plot2d([
>     [discrete, xtvalues],
>     [discrete, ytvalues]
> ]);
>
>
>
> ------------------------------
>
> Message: 10
> Date: Fri, 27 Mar 2009 08:29:24 +0000
> From: Jaime Villate <villate at fe.up.pt>
> Subject: Re: [Maxima] plotting results of rk
> To: Semyon Glazyrin <siglazyrin at gmail.com>
> Cc: maxima at math.utexas.edu
> Message-ID: <1238142564.6447.4.camel at bethe>
> Content-Type: text/plain
>
> On Sex, 2009-03-27 at 08:54 +0300, Semyon Glazyrin wrote:
> > $ load(dynamics);
> > $ res:rk([x+y,x-y],[x,y],[1,0],[t,0,1,0.1]);
> >
> > This returns me a list of [t,x,y]. How can I plot x(t) and y(t)?
>
> y(t):
>
> plot2d([discrete, makelist([res[i][1],res[i][3]], i, 1, length(res))])$
>
> x(t), with axes labels:
>
> plot2d([discrete, makelist([res[i][1],res[i][2]], i, 1, length(res))],
> [xlabel,"t"], [ylabel,"y"])$
>
> Regards,
> Jaime
>
>
>
> ------------------------------
>
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>
>
> End of Maxima Digest, Vol 32, Issue 51
> **************************************
>



-- 
http://www.watch-movies-online-hollywoodkiller.com