Hi everyone,
I have been trying to use maxima not more than one month so I
definately need some help even with the simplest tasks
Anyway my problem is about matrices.
I had defined some matrices which has variables inside and I am trying
to do simple calculations with them.
I have couple of 6x6 matrices and as I said the calculations between
those are very simple(only addition,substraction,multiplication).The
problem is I can not get the resultant matrix alone.You can see it
much better with the example below;
S1:(IDENT(6)-(Contact1_Coupling_FloorWheel1 .
Contact1_Coupling_Wheel1Floor)) . Contact1_Force_Floor;
All of these matrices are defined (like
Contact1_Coupling_FloorWheel1:matrix([1,0,0,0,0,1],[cos(B),1,-1,r,y,y],[-sin(A1),cos(A1),cos(B),0,0,1],[0,0,0,0,0,0],[r,0,0,0,0,0],[1,0,0,0,0,1])
After these calculation S1 appears including the other matrixes but
not as a whole matrix itself.
I would be realy pleased if anyone can help me with this problem.
P.S:apperance of S1 after calculations;
MATRIX([1, 0, 0, 0, 0, 0],
[0, COS(A1), - SIN(A1), 0, 0, 0],
[0, SIN(A1), COS(A1), 0, 0, 0], [0, 0, 0, 1, 0, 0],
[0, 0, 0, 0, COS(A1), - SIN(A1)],
[0, 0, 0, 0, SIN(A1), COS(A1)])
. (MATRIX([- mW a], [- mW a
- SIN(A) g mW], [- mW a - COS(A) g mW],
[betaT1 - Iwx alfa], [- (Iwx - Iwy) Wheel1_w
Wheel1_w - Iwy alfa],
[- (Iwy - Iwx) Wheel1_w w - Iwy alfa])
+ MATRIX([1, 0, 0, 0, 0, 0], [0, COS(A), SIN(A), 0, 0,
0], [0, - SIN(A), COS(A), 0, 0, 0], [0, R, 0, 1, 0, 0],
[- COS(A) R, 0, 0, 0, COS(A), SIN(A)],
[SIN(A1) R, 0, 0, 0, - SIN(A), COS(A)])
. Contact1_Force_Floor)
On Tue, 1 Feb 2005 16:05:20 -0600, maxima-request@math.utexas.edu
wrote:
> Send Maxima mailing list submissions to
> maxima@www.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@www.math.utexas.edu
>
> You can reach the person managing the list at
> maxima-admin@www.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: Working with more than 16 digits? (Barton Willis)
> 2. Re: Working with more than 16 digits? (Albert Reiner)
> 3. Re: Working with more than 16 digits? (Miquel Cabanas)
> 4. Re: Maxima problem with ASSUME_POS:TRUE; (Robert Dodier)
> 5. Re: maxima translations - newbie question (Vadim V. Zhytnikov)
> 6. Re: Maxima problem with ASSUME_POS:TRUE; (Antonio Amorim)
> 7. fast_linsolve (Fabrizio Caruso)
> 8. Bug in today's CVS (David Ronis)
> 9. Re: Bug in today's CVS (Barton Willis)
> 10. Re: Bug in today's CVS (Richard Fateman)
> 11. Re: fast_linsolve (Wolfgang Jenkner)
> 12. macsyma on sale on ebay! (Richard Fateman)
> 13. Re: macsyma on sale on ebay! (Paulo Ney de Souza)
> 14. Re: macsyma on sale on ebay! (C Y)
>
> --__--__--
>
> Message: 1
> To: Albert Reiner
> Cc: Maxima mailing list
> Subject: Re: [Maxima] Working with more than 16 digits?
> From: Barton Willis
> Date: Mon, 31 Jan 2005 09:45:16 -0600
>
> Three things:
>
> (1) The value of fpprec mysteriously affects the number of displayed
> digits
> for double floats.
>
> (2) Doing ev(e, numer) evaluates 'e 'to a double (or least it does
> sometimes).
>
> (3) To evaluate to a big float, use the function 'bfloat'
>
> (%i1) bfloat(%pi), fpprec : 50;
> (%o1) 3.1415926535897932384626433832795028841971693993751B0
>
> Barton
>
> Albert Reiner
> Sent by: maxima-admin@math.utexas.edu
> 01/31/2005 09:25 AM
>
> To: Maxima mailing list
> cc:
> Subject: [Maxima] Working with more than 16 digits?
>
> Hi,
>
> the documentation of FPPREC:
>
> ,----
> | - Variable: FPPREC
> | default: [16] - Floating Point PRECision. Can be set to an
> | integer representing the desired precision.
> `----
>
> seems to suggest that setting it to something higher than 16 would
> result in a higher working precision than default double precision
> numbers. However, while
>
> ,----
> | block([fpprec: 2], print(ev(%pi, numer))) $
> |
> | 3.14
> `----
>
> ...
>
> ,----
> | block([fpprec: 14], print(ev(%pi, numer))) $
> |
> | 3.14159265358979
> `----
>
> ,----
> | block([fpprec: 15], print(ev(%pi, numer))) $
> |
> | 3.141592653589793
> `----
>
> all work more or less as expected, a higher fpprec does not print out
> more digits as, e.g., N[Pi, 500] does in Mathematica or something like
> N(Pi, 500); does in yacas.
>
> Is there a simple way to get, say, the first 25 digits of pi in
> Maxima?
>
> Thanks in advance,
>
> Albert.
>
> _______________________________________________
> Maxima mailing list
> Maxima@www.math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>
> --__--__--
>
> Message: 2
> To: Maxima mailing list
> Subject: Re: [Maxima] Working with more than 16 digits?
> Cc: Barton Willis
> From: Albert Reiner
> Date: 31 Jan 2005 16:57:13 +0100
>
> Thanks for the clarifications!
>
> [Barton Willis , Mon, 31 Jan 2005 09:45:16 -0600]:
> > Three things:
> >
> > (1) The value of fpprec mysteriously affects the number of displayed
> > digits
> > for double floats.
>
> That mystery seems to be documented:
>
> ,----
> | - Variable: FPPRINTPREC
> | default: [0] - The number of digits to print when printing a
> | bigfloat number, making it possible to compute with a large number
> | of digits of precision, but have the answer printed out with a
> | smaller number of digits. If FPPRINTPREC is 0 (the default), or >=
> | FPPREC, then the value of FPPREC controls the number of digits used
> | for printing. However, if FPPRINTPREC has a value between 2 and
> | FPPREC-1, then it controls the number of digits used. (The minimal
> | number of digits used is 2, one to the left of the point and one to
> | the right. The value 1 for FPPRINTPREC is illegal.)
> `----
>
> Thanks again,
>
> Albert.
>
> --__--__--
>
> Message: 3
> Date: Mon, 31 Jan 2005 17:17:50 +0100
> From: Miquel Cabanas
> Subject: Re: [Maxima] Working with more than 16 digits?
> To: Albert Reiner
> Cc: Maxima mailing list
> Reply-to: Miquel.Cabanas@uab.es
>
> hi,
>
> On Mon, Jan 31, 2005 at 04:25:25PM +0100, Albert Reiner wrote:
> >
> > ,----
> > | block([fpprec: 15], print(ev(%pi, numer))) $
> > |
> > | 3.141592653589793
> > `----
>
> use bfloat instead of numer, i.e.
>
> (C1) block([fpprec: 32], print(ev(%pi, bfloat))) $
>
> 3.1415926535897932384626433832795B0
>
> (C2) block([fpprec: 128], print(ev(%pi, bfloat))) $
>
> 3.1415926535897932384626433832795028841971693993751058209749445923078164062862#
>
> 089986280348253421170679821480865132823066470938446B0
>
> Miquel
>
> --
> Miquel E Cabanas ------------------------------------------------------
> SeRMN, Universitat Autonoma de Barcelona (Miquel.Cabanas@uab.es)
> ------------------------------------------o-oo--ooo---ooo--oo-o--------
>
> --__--__--
>
> Message: 4
> Date: Mon, 31 Jan 2005 08:50:08 -0800 (PST)
> From: Robert Dodier
> To: Antonio.Amorim@fisica.fc.ul.pt
> Cc: maxima@math.utexas.edu
> Subject: Re: Maxima problem with ASSUME_POS:TRUE;
>
> Hi Antonio,
>
> I saw your message about the Maxima plug-in for your wiki.
> That is a very interesting idea.
>
> Given that one cannot foresee all the conditions in which
> Maxima will want to ask a question, I don't see that it's
> practical to use "assume" to try to preclude all questions.
>
> How about this. Instead of directly entering Maxima commands
> in the wiki text, open a little window in which users can
> interactively develop their answer, then cut and paste the
> result into the wiki page. It always takes some working
> with different expressions to get stuff just the way you
> want it, anyway.
>
> Incidentally
>
> assume_pos: true$
> integrate (x*log(x), x, a, a+da);
>
> yields the expected result without asking a question.
> In this case the difference, (a+da) - a = da, is a symbol.
> "? assume_pos" seems to imply (it's not very clear) that only
> symbols and subscripted variables are covered, not all expressions.
>
> Hope this helps,
> Robert Dodier
>
> __________________________________
> Do you Yahoo!?
> Take Yahoo! Mail with you! Get it on your mobile phone.
> http://mobile.yahoo.com/maildemo
>
> --__--__--
>
> Message: 5
> Date: Mon, 31 Jan 2005 22:27:26 +0300
> From: "Vadim V. Zhytnikov"
> To: Marco Ciampa , Maxima List
> Subject: Re: [Maxima] maxima translations - newbie question
>
> Marco Ciampa writes:
>
> >On Sun, Jan 30, 2005 at 12:07:52AM +0300, Vadim V. Zhytnikov wrote:
> >
> >
> >>Tcl has it's own translation mechanism called msgcat
> >>(since version 8.2 AFAIK). In the Maxima CVS source tree
> >>you can see two files:
> >>
> >>/plotting/msgs/ru.msg
> >>/interfaces/xmaxima/msgs/ru.msg
> >>
> >>
> >I've seen that but I'm trying to push the use of standard mechanism and
> >tools for translating maxima.
> >Why change and use gettext?
> >For many reasons...
> >
> > - It's easyer to update the translations using tools like msgmerge or
> > msguniq
> > - You can use standard tools for translating like emacs (po-mode),
> > kbabel, poedit etc that really improve productivity of about
> > 300%!
> > - You can use different encodings like UTF-8, iso-8859-1, etc
> >
> >etc.
> >
> >
> >
> Ok, all these advanced tools are very nice. And if
> you find it helpful in making new translation then fine!
> As far as I understand all things listed in
> http://groups.yahoo.com/group/tcl_announce/message/1750
> are just to simplify creation of new .msg files.
> But all xmaxima/omplotdata needs is a set of final .msg
> files for different languages and I really don't see any
> pressing reason to include any extra machinery in Maxima
> distribution.
>
> BTW, translation string in .msg should be in UTF-8.
> This is how it works for Russian - one ru.msg file
> gives me correct translation in all encodings I tried,
> even on Windows.
>
> There is yet one more strategical reason why I don't like
> to fiddle with xmaxima too much. I don't expect it evolve.
> We are expecting to have new standard graphical frontend
> for Maxima in some near future. On the other hand xmaxima
> is still there for a couple of years, especially on Windows,
> and making translation (not big deal - just a 20Kb work)
> may be helpful.
>
> --
> Vadim V. Zhytnikov
>
>
>
>
> --__--__--
>
> Message: 6
> Date: Tue, 01 Feb 2005 21:40:09 +0000
> From: Antonio Amorim
> To: Robert Dodier
> CC: maxima@math.utexas.edu
> Subject: Re: Maxima problem with ASSUME_POS:TRUE;
>
> Robert Dodier escreveu:
>
> >Hi Antonio,
> >
> >I saw your message about the Maxima plug-in for your wiki.
> >That is a very interesting idea.
> >
> >Given that one cannot foresee all the conditions in which
> >Maxima will want to ask a question, I don't see that it's
> >practical to use "assume" to try to preclude all questions.
> >
> >How about this. Instead of directly entering Maxima commands
> >in the wiki text, open a little window in which users can
> >interactively develop their answer, then cut and paste the
> >result into the wiki page. It always takes some working
> >with different expressions to get stuff just the way you
> >want it, anyway.
> >
> >Incidentally
> >
> > assume_pos: true$
> > integrate (x*log(x), x, a, a+da);
> >
> >yields the expected result without asking a question.
> >In this case the difference, (a+da) - a = da, is a symbol.
> >"? assume_pos" seems to imply (it's not very clear) that only
> >symbols and subscripted variables are covered, not all expressions.
> >
> >Hope this helps,
> >Robert Dodier
> >
> >
> >
> >
> >__________________________________
> >Do you Yahoo!?
> >Take Yahoo! Mail with you! Get it on your mobile phone.
> >http://mobile.yahoo.com/maildemo
> >
> >
> Dear Robert,
> Thank you for your message. You actually gave me a great idea. I will
> try to investigate if I can recover the questions that Maxima is posing
> and display them in the page. In this way only when there are no
> questions would the user get the result formula.
> Regarding the maxima WIKI plugin, the nice thing about it is that the
> image png is actually stored with a name that is the md5 signature of
> the input expression. This makes it very easy to find if the expression
> needs to be evaluated or if a png result already exists.
>
> All the best,
> Antonio Amorim
>
> --__--__--
>
> Message: 7
> Date: Tue, 1 Feb 2005 16:31:54 +0100 (CET)
> From: Fabrizio Caruso
> To: maxima@math.utexas.edu
> Subject: fast_linsolve
>
> Dear Maxima Mailing List
>
> The in line manual shows an entry for
> "fast_linsolve" which is (or was) an optimized
> linear system solver for sparse matrices.
>
> For me it would be great to have it
> because my summation package would
> have a considerable speed-up.
> I am sure many other packages would
> take advantage of it.
>
> Unfortunately "fast_linsolve" seems
> to have disappeared from Maxima.
> Does someone know what happened to it?
> Is it still there under a different name?
> Anyone working on sparse linear systems
> (where entries can possibly be rational functions)?
>
> Regards
>
> Fabrizio
>
> --__--__--
>
> Message: 8
> From: David Ronis
> Date: Tue, 1 Feb 2005 12:52:09 -0500
> To: maxima@math.utexas.edu
> Subject: Bug in today's CVS
> Reply-To: david.ronis@mcgill.ca
>
> I'm trying to use maxima to do a saddle point integration with the following:
>
> f(x):=x-2*log((exp(x)+1)/2) + x**2/4;
> exp(-x**2/4)*taylor(exp(f(x))*(x*kb+mu)**gamma,x,0,4);
> integrate(%,x,%minf,%inf)/4;
>
> As far as I can tell, the taylor expansion is done correctly, however,
> the multiplication by the additional exponential in the second line
> seems to be ignored. The resulting integral doesn't converge and I'm
> left with explicit %inf/%minf's at the integration (which is wrong in
> any event).
>
> David
>
> --__--__--
>
> Message: 9
> To: david.ronis@mcgill.ca
> Cc: maxima@math.utexas.edu
> Subject: Re: [Maxima] Bug in today's CVS
> From: Barton Willis
> Date: Tue, 1 Feb 2005 12:14:55 -0600
>
> Taylor expansions are contagious. When maxima multiplies exp(-x^2/4)
> times a
> Taylor polynomial, the product gets expanded as a Taylor polynomial. To
> do the calculation, use 'ratsimp' to convert from a Taylor polynomial to
> a non-taylor.
>
> Also, I think you want 'minf and 'inf instead of %minf and %inf.
>
> (%i3) kill(all)$
> (%i1) f(x):=x-2*log((exp(x)+1)/2) + x**2/4$
> (%i2) exp(-x**2/4)*ratsimp(taylor(exp(f(x))*(x*kb+mu)**gamma,x,0,4))$
> (%i3) display2d : false$
> (%i4) integrate(%o2,x,'minf,'inf)/4;
> Is gamma+4 positive, negative, or zero?
>
> pos;Is mu positive or negative?
> pos;Is gamma an integer?
> no;Is gamma+3 positive, negative, or zero?
> pos;Is kb*gamma positive, negative, or zero?
> pos;Is gamma+2 positive, negative, or zero?
> pos;Is kb^2*gamma positive, negative, or zero?
> pos;Is gamma+1 positive, negative, or zero?
> pos;Is gamma positive, negative, or zero?
> pos;Is kb zero or nonzero?
> nonzero;Is kb positive or negative?
>
> pos;(%o4) (24*sqrt(%pi)
> *(8*kb^4*mu^gamma*gamma^4-48*kb^4*mu^gamma*gamma^3
> +88*kb^4*mu^gamma*gamma^2
> -48*kb^4*mu^gamma*gamma+2*mu^(gamma+4))
> +4*sqrt(%pi)*(96*kb^2*mu^(gamma+2)*gamma^2-96*kb^2*mu^(gamma+2)*gamma)
> +384*sqrt(%pi)*mu^(gamma+4))
> /(768*mu^4)
>
> Barton
>
> David Ronis
> Sent by: maxima-admin@math.utexas.edu
> 02/01/2005 11:52 AM
> Please respond to david.ronis
>
> To: maxima@math.utexas.edu
> cc:
> Subject: [Maxima] Bug in today's CVS
>
> I'm trying to use maxima to do a saddle point integration with the
> following:
>
> f(x):=x-2*log((exp(x)+1)/2) + x**2/4;
> exp(-x**2/4)*taylor(exp(f(x))*(x*kb+mu)**gamma,x,0,4);
> integrate(%,x,%minf,%inf)/4;
>
> As far as I can tell, the taylor expansion is done correctly, however,
> the multiplication by the additional exponential in the second line
> seems to be ignored. The resulting integral doesn't converge and I'm
> left with explicit %inf/%minf's at the integration (which is wrong in
> any event).
>
> David
>
> _______________________________________________
> Maxima mailing list
> Maxima@www.math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>
> --__--__--
>
> Message: 10
> Date: Tue, 01 Feb 2005 10:29:25 -0800
> From: Richard Fateman
> Organization: University of California, Berkeley
> To: david.ronis@mcgill.ca
> CC: maxima
> Subject: Re: [Maxima] Bug in today's CVS
>
> What is %minf? Do you mean minf? Same for inf.
>
> Usually people use x^y rather than x**y.
>
> When you take a taylor series Q(x) and multiply it
> by exp(x), the exp(x) is converted to a taylor series
> and multiplied.
> Assuming you do not want this, you might try
> something like
>
> expand (exp(-x^2/4)* ratdisrep(taylor ...)))
>
> integrate(%,x,minf,inf)
>
> David Ronis wrote:
>
> > I'm trying to use maxima to do a saddle point integration with the following:
> >
> > f(x):=x-2*log((exp(x)+1)/2) + x**2/4;
> > exp(-x**2/4)*taylor(exp(f(x))*(x*kb+mu)**gamma,x,0,4);
> > integrate(%,x,%minf,%inf)/4;
> >
> >
> > As far as I can tell, the taylor expansion is done correctly, however,
> > the multiplication by the additional exponential in the second line
> > seems to be ignored. The resulting integral doesn't converge and I'm
> > left with explicit %inf/%minf's at the integration (which is wrong in
> > any event).
> >
> > David
> >
> > _______________________________________________
> > Maxima mailing list
> > Maxima@www.math.utexas.edu
> > http://www.math.utexas.edu/mailman/listinfo/maxima
>
> --__--__--
>
> Message: 11
> To: Fabrizio Caruso
> Cc: maxima@math.utexas.edu
> Subject: Re: [Maxima] fast_linsolve
> From: Wolfgang Jenkner
> Date: Tue, 01 Feb 2005 18:47:30 +0100
>
> Fabrizio Caruso writes:
>
> > Dear Maxima Mailing List
> >
> > The in line manual shows an entry for "fast_linsolve" which is (or
> > was) an optimized linear system solver for sparse matrices.
>
> > Unfortunately "fast_linsolve" seems to have disappeared from Maxima.
> > Does someone know what happened to it? Is it still there under a
> > different name?
>
> It's still there in share/affine (written by Bill Schelter). The CVS
> version is currently broken (mainly because of the sloop -> loop
> change). The version contained in the last stable release is supposed
> to contain some unspecified but reasonable amount of working stuff.
>
> If you use the stable version try `load("affine.lisp");'. Probably,
> this won't work because of some brain-dead logical-pathname stuff of
> mine. Plan B: Batch a file containing the following
>
> load("affine/compat");
> load("affine/amacros");
> load("affine/dummy-gc");
> load("affine/polybas");
> load("affine/sparsemat");
> load("affine/aquotient");
> load("affine/new-rat");
> load("affine/polya");
> load("affine/ndotsimp");
> load("affine/polyb");
> load("affine/polysmp");
> load("affine/sub-proj");
> load("affine/polyc");
> load("affine/polyd");
> load("affine/sheafa");
> load("affine/sheafb");
> load("affine/sheafc");
> load("affine/dim-3");
> load("affine/modsimp");
>
> Wolfgang
>
> --__--__--
>
> Message: 12
> Date: Tue, 01 Feb 2005 13:24:41 -0800
> From: Richard Fateman
> Organization: University of California, Berkeley
> To: maxima
> Subject: macsyma on sale on ebay!
>
> http://cgi.ebay.com/ws/eBayISAPI.dll?ViewItem&category=3789&item=7130355446&rd=1&ssPageName=WDVW
>
> $500. I just came across this and thought you guys might be
> interested.
>
> Don't rush, there are 99 copies available.
>
> RJF
>
> --__--__--
>
> Message: 13
> Date: Tue, 1 Feb 2005 13:46:58 -0800 (PST)
> From: Paulo Ney de Souza
> To: fateman@cs.berkeley.edu, maxima@math.utexas.edu
> Subject: Re: [Maxima] macsyma on sale on ebay!
>
> What a bargain!
>
> Paulo Ney
>
> >From maxima-admin@math.utexas.edu Tue Feb 1 13:25:32 2005
> >From: Richard Fateman
> >To: maxima
> >Subject: macsyma on sale on ebay!
> >
> >http://cgi.ebay.com/ws/eBayISAPI.dll?ViewItem&category=3789&item=7130355446&rd=1&ssPageName=WDVW
> >
> >$500. I just came across this and thought you guys might be
> >interested.
> >
> >Don't rush, there are 99 copies available.
> >
> >RJF
> >
> >_______________________________________________
> >Maxima mailing list
> >Maxima@www.math.utexas.edu
> >http://www.math.utexas.edu/mailman/listinfo/maxima
> >
>
> --__--__--
>
> Message: 14
> Date: Tue, 1 Feb 2005 13:59:17 -0800 (PST)
> From: C Y
> Subject: Re: [Maxima] macsyma on sale on ebay!
> To: maxima
>
> Heh - I guess David Schmidt is 71271.1417 on ebay. Well, it's the
> first active marketing I've seen from them since 1999, so that's gotta
> be a good sign. 99 copies available after almost eight days... I
> guess there isn't a rush to buy $500 software with no hard copy
> documentation!
>
> Hmm, new spiffy website for http://www.symbolics.com too. Looks like
> anybody wanting to upgrade to 2.4 from an earlier version can do it for
> $250.
>
> I hope they start active development of Macsyma again - it would really
> be a shame to have it remain static indefinitely. Maybe they could at
> least started digging out the Linux version from whatever server it's
> sitting on and help some of those folks that can't use their older
> purchased copies.
>
> CY
>
> --- Richard Fateman wrote:
>
> >
> http://cgi.ebay.com/ws/eBayISAPI.dll?ViewItem&category=3789&item=7130355446&rd=1&ssPageName=WDVW
> >
> > $500. I just came across this and thought you guys might be
> > interested.
> >
> > Don't rush, there are 99 copies available.
> >
> > RJF
> >
> > _______________________________________________
> > Maxima mailing list
> > Maxima@www.math.utexas.edu
> > http://www.math.utexas.edu/mailman/listinfo/maxima
> >
>
> __________________________________
> Do you Yahoo!?
> All your favorites on one personal page â Try My Yahoo!
> http://my.yahoo.com
>
> --__--__--
>
> _______________________________________________
> Maxima mailing list
> Maxima@www.math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>
> End of Maxima Digest
>