Gunter Königsmann | 25 Apr 08:32 2014
Picon

Another patch for draw


This truly is the last patch for now, I promise:
The documentation did tell that a terminal isn't supported - even if
recent versions of draw actually support it.

Kind regards,

    Gunter.
------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Gunter Königsmann | 25 Apr 08:07 2014
Picon

More line types for draw


Don't know how to allow users to use numbers as line types, as well so
this is the last patch for now. But did at least try to make all line
types accessible I need.

Kind regards,

  Gunter.
Attachment (0003-More-line-types.patch): text/x-patch, 2278 bytes
------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Gunter Königsmann | 25 Apr 07:47 2014
Picon

patch: Documented more line styles


Another thing that was doable without deep lisp knowledge.
------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Gunter Königsmann | 25 Apr 07:40 2014
Picon

Patch: Make the line style "dashed" work on more terminals


This task seemed doable even by a non-lisp mage so I gave it a try.
For the wxt terminal maxima sometimes did explicitely deactivate
dashes via "solid" and sometimes didn't. Fixed that, as well.

Testcase:

draw2d(   terminal=pdfcairo,line_width=2.25,color=black,
                dimensions=[1200,800],
                file_name="Test.pdf",line_type=dashes,
          explicit(sin(x),x,1,3)
);

now creates a dashed line instead of a solid one.

Kind regards,

  Gunter.
------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
_______________________________________________
Maxima-discuss mailing list
(Continue reading)

Dimiter Prodanov | 25 Apr 01:15 2014
Picon

delta function in maxima

Dear all,

I have written a short package implementing Dirac's delta function.

https://github.com/dprodanov/maximadelta

The rules can be used during integration and differentiation of expressions containing delta functions.

best regards,

Dimiter
------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Pankaj Sejwal | 24 Apr 15:31 2014
Picon

iterate over two lists

May be you can riffle information like this,
--------------------------------------------------------
riffle(x):=block([mat:transpose(apply(matrix,x)),out:[]],for i:1 thru first(matrix_size(mat)) do  (push(args(mat[i]),out)),reverse(out))$

Usage :: riffle([[a,b,c],[d,e,f],[g,h,i]]);
OUT :: [[a,d,g],[b,e,h],[c,f,i]]
----------------------------------------------------------
Than you can map 'print' or any other function.
--
Regards,
Pankaj Sejwal
_______________________________________________
"The more I read, the more I acquire, the more certain I am that I know nothing.” -
 Voltaire
------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Dan Gildea | 24 Apr 13:18 2014
Picon
Picon

quad_qawf nondeterministic

I intermittently get the following error in the maxima 
test suite with ecl 13.5.1:

Running tests in rtest8: 
********************** Problem 104 ***************
Input:
ev(e5, au = 0, omega = 2)

Result:
            - u
quad_qawf(%e   , u, 0, 2, sin, epsabs = 10.e-11, limit = 32, maxp1 = 100, 
                                                                   limlst = 10)

This differed from the expected result:
[0.4000000000000001, 22.16570948815925e-12, 175, 0]

If I run the following command at the maxima command line several
times in a row, I usually get the first result below, but
sometimes the second:

(%i46) ev( quad_qawf (exp (-u), u, au, omega, sin, limit=32), au=0, omega=2);

(%o46) [0.4000000000000001,22.165709488159253e-12,175,0]
(%i47) ev( quad_qawf (exp (-u), u, au, omega, sin, limit=32), au=0, omega=2);

(%o47) quad_qawf(%e^-u,u,0,2,sin,epsabs = 10.e-11,limit = 32,maxp1 = 100,
                 limlst = 10)

------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
Pankaj Sejwal | 24 Apr 10:41 2014
Picon

Re: Maxima-discuss Digest, Vol 3, Issue 44

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
I would start by looking at, say, share/lapack/dgeqrf.lisp or
eigensys.lisp to see how to implement an interface to the underlying
lapack routines.  This requires some knowledge of lisp, some knowledge
of how Fortran works, some knowledge of lapack, and some knowledge of
how Maxima stores matrices.
I can help with that.
Actually, there's a huge amount of boilerplate stuff and perhaps a lot
of this could be hidden behind some macros to make it more uniform and
consistent. This will also help with writing new interfaces.  I can
try implementing that.
Ray
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
I will investigate it and learn things you suggested.
Looks like it is going to take some good time  
seeing how learned I am with all that is referred to.

I am too thankful for your grace to offer me help.
-- 
Regards,
Pankaj Sejwal
_______________________________________________
"The more I read, the more I acquire, the more certain I am that I know nothing.” - Voltaire
------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Helmut Jarausch | 24 Apr 10:05 2014
Picon
Picon

How to iterate over 2 list?

Hi,

given two list tlist and ylist of equal length, how can I iterate of  
both list like
  for (t,y) in (tlist,ylist) do ....

My current solution is not appealing:

tlist:[1,2,3]$ ylist:[1,2,9]$
yc:copylist(ylist)$
for t in tlist do block([y],
   y:first(yc), yc:rest(yc),
   print(t,y)
);

Many thanks for a hint,
Helmut
------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
Stavros Macrakis | 24 Apr 04:33 2014
Picon

Subscripted variables WAS: How to make buildq evaluate its first argument - please help

On Wed, Apr 23, 2014 at 9:40 PM, Robert Dodier <robert.dodier <at> gmail.com> wrote:
On 2014-04-23, Stavros Macrakis <macrakis <at> alum.mit.edu> wrote:

> * not all of Maxima allows x[1] where it allows a variable
> ** for x[1] thru 3 => syntax error
> ** function, block, and lambda formals (do we really want x[i] and x[2] as
> formals?)
> ** ... I'm sure there are more

declare doesn't allow subscripted variables.

​That's an interesting one. If I say declare(x[i],integer), does that mean x[i] for all i? or just the literal i? or just expressions which are equal to i? (e.g. x[(i^2-1)/(i-1)-1]) If the last, what is the criterion for equality? Does assume(equal(i,j)) mean that x[j] is now declared integer? What if I declare x[i] to be an integer and later assume(equal(i,j)) or locally assume(equal(i,j))? What if I forget() the assumption that equal(i,j)?
For the record, I don't see any reason to disallow subscripted variables
​ ​
as local variables.
​ 
I can't think of any language that does allow subscripted variables as local variables.

Consider (x:[a,b,c],block([x[2]],x:[])? What is the value of x after that?

What about (x:[a,b,c],y:x,block([x[2]],y[2]:5))? The obvious implementation would result in x=[a,b,c], y=[a,b,c]. Is that reasonable? Or should the result be x=[a,b,c], y=[a,5,c], at the cost of changing the sharing relationships?

Maxima has computing/memoized subscripts. What if x[0]:1, x[i]:=i*x[i-1]. Now what is the value of block([x[0]:2],x[3])? This will depend on whether x[1] has been memo'ized before or not....  Is that a good idea, especially because now the value of x[1] outside the block will depend on this.

In short, I don't think allowing subscripted variables as local will simplify anything; rather it will make the language harder to understand and more implementation-dependent.

                    -s


------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Edwin Woollett | 23 Apr 22:04 2014
Picon

Chapter 2 Computational Physics Using R and Maxima


   Chapter 2, Initial Value Problems (Computational Physics Using
            R and Maxima) has been posted to my web page.

   From that web page:
-----------------------------------
Chapter 2, Initial Value Problems

 A variety of initial value problems are solved using both
    R and Maxima. Homemade codes for the Euler method and 
    the classical fourth order Runge-Kutta method are 
    demonstrated with numerous examples. The standard
    Maxima methods rk and rkf45 are discussed, and the
    R deSolve package ode method is extensively discussed
    with many examples. Finally, methods for solving initial
    value problems in which the derivatives are discontinuous
    and in which the dependent variables are discontinous are discussed.

    --cp2.pdf : Ch. 2, April 23, 2014, Maxima 5.31.0, R 3.0.2, 71 pages
    --cp2.tex : Ch. 2 Latex code file
    --k2util.mac : Maxima utility functions for Ch. 2, April 23, 2014
    --myode.mac : Maxima code for Ch. 2 , April 23, 2014, Maxima 5.31
    --myode.R : R code for Ch. 2, April 23, 2014, R 3.0.2
    --myrkf45.mac : Maxima code for Ch. 2 , April 23, 2014, Maxima 5.31 

Chapter 2 Topics

 1   The Euler Method, Truncation Errors, Round-off Errors, Instability
 2   4th Order Runge-Kutta Code for R and Maxima
 3   The Maxima Methods rk and rkf45
 4   The R package deSolve and Use of ode
 5   Solving O.D.E.'s with External Forcing Data
 6   Integrating O.D.E.'s with Discontinuous Derivatives
 7   Integrating O.D.E.'s with Discontinuous Dependent Variables Using R
      1  Events Specified with a data.frame
      2  Intravenous Drug Injection Model
      3  Using an Event Function at Specified Times
      4  Using an Event Function when a Root Condition is Satisfied
      5  Use of a Switching Parameter as a Dependent Variable 
--------------------------------------------------------------------------

Ted Woollett

http://www.csulb.edu/~woollett/

------------------------------------------------------------------------------
Start Your Social Network Today - Download eXo Platform
Build your Enterprise Intranet with eXo Platform Software
Java Based Open Source Intranet - Social, Extensible, Cloud Ready
Get Started Now And Turn Your Intranet Into A Collaboration Platform
http://p.sf.net/sfu/ExoPlatform

Gmane