Re: [sage-support] Re: Point-set topology in sage?

2010-08-05 Thread Alec Battles
Update: I think the G. Carlsson paper, Topology and Data, is quite a find. Thanks for suggesting the link, Justin. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options

Re: [sage-support] Re: Point-set topology in sage?

2010-08-05 Thread Alec Battles
Thanks. It looks promising. My earlier search didn't yield anything. Alec On Thu, Aug 5, 2010 at 2:13 PM, Justin C. Walker wrote: > > On Jul 29, 2010, at 06:12 , Alec Battles wrote: > >> Not specifically, no, but orms.mfo.de (an excellent list with frequent >> updat

Re: [sage-support] Re: Point-set topology in sage?

2010-07-29 Thread Alec Battles
Not specifically, no, but orms.mfo.de (an excellent list with frequent updates) lists 9 projects under topology. I will do some digging later, but feel free to help me search. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email

Re: [sage-support] Re: Point-set topology in sage?

2010-07-28 Thread Alec Battles
't see why it shouldn't be possible to do point-set topology with a computer. Is there something I'm missing? Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For

Re: [sage-support] Re: Graphics

2010-07-28 Thread Alec Battles
That's great info. Thanks. A On Wed, Jul 28, 2010 at 5:58 PM, kcrisman wrote: > > > On Jul 28, 11:51 am, Alec Battles wrote: >> Sorry if there's some good reason why this would be redundant, but I >> have an idea. It's a bit scary  that you can't simply

Re: [sage-support] Re: Graphics

2010-07-28 Thread Alec Battles
pykig.py), and draws geometry very well using scripted commands. It would give the program another input-output option, as well as bolstering sage's capabilities. I'd certainly make great use of it. Does that make sense? Anyway, KIG may be more what you're looking for, Santanu. Al

[sage-support] Point-set topology in sage?

2010-07-28 Thread Alec Battles
I've tried to find info about this. Can't find anything... Any suggestions/help? Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this grou

Re: [sage-support] Re: Does sage have support for category theory?

2010-07-27 Thread Alec Battles
Simon Thanks for sending so much great info. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support

[sage-support] Does sage have support for category theory?

2010-07-25 Thread Alec Battles
I've tried to find info on this. Sorry if it's already been written down. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this grou

Re: [sage-support] Re: Turning axes off

2010-07-14 Thread Alec Battles
Must be better than nothing. I'll try it. Thank you. Alec On Thu, Jul 15, 2010 at 12:19 AM, kcrisman wrote: > Dear Alec, > > The following will (partly) answer your question if you use 'axes' > instead of 'frame'.  I'm not sure (vis-a-vis that thr

Re: [sage-support] Re: adding more points to an existing polyhedron

2010-07-14 Thread Alec Battles
Sounds like an interesting project. Is it for school or for fun? I'm using Sage to make art and sorta failing. If I knew more programming I'd give it a lot more time, but sadly my only talents are in the GIMP/Inkscape department. Alec -- To post to this group, send email to sa

[sage-support] Turning axes off

2010-07-12 Thread Alec Battles
I am using SAGE, not so much as a mathematical tool, but to make art. With that in mind, I want to disable axes pretty much permanently. I am currently plotting complex functions such as sage: f(z) = z^5 + z - 1 + 1/z sage: complex_plot(f, (-3, 3), (-3, 3)) Thanks, Alec -- To post to this

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
On Apr 16, 9:48 pm, David Joyner wrote: > Sorry Alec, maybe there is a language problem here. Sorry, David, but there may be also a cultural problem as well. You, obviously, think that, being a moderator of sage-edu, you could reject your colleagues posts there, as you rejected my post th

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
On Apr 16, 8:55 pm, Alec Mihailovs wrote: > Maybe, I still don't understand it right, but it seems as if you are > saying, that the problem is to do a similar thing in cases when the > solution is not unique, producing the list of all the solutions. That > can be done similarl

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
en multiplying each of them by G, as in 'correct', If I again understood it wrong, could you, please, give a simple example (with smaller sizes), with the correct answer? Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this grou

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
On Apr 16, 5:28 pm, Alec Mihailovs wrote: > I looked at the C.decode (for C in the original post), and it does > something different - it corrects the errors. Perhaps, the better name > for it would be C.correct ? If that is needed, that could be done as def correct(w,erasures=[],num

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
On Apr 16, 5:28 pm, Alec Mihailovs wrote: > > Also, perhaps, just simple linsolve can be used (in my decode) if it > exists in Sage and works for overdetermined systems of linear > equations. Yes, I found it - it is called solve_left, and it can be used as follows, def decode

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
exists in Sage and works for overdetermined systems of linear equations. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.

[sage-support] Re: Vector solutions, a basic question.

2010-04-16 Thread Alec Mihailovs
x form A*x=b in which case the solution can be obtained as x=A.inverse()*b That should work better than solve. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options,

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
On Apr 16, 4:16 pm, Alec Mihailovs wrote: > For example, > > m=V.random_element(); m > >                 (2, 1, 0, 0, 1, 2, 0, 1, 0, 0) It should be V=VectorSpace(GF(3),10) before that. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
(2, 1, 0, 0, 1, 2, 0, 1, 0, 0) decode(w,range(15)) Too many erasures Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
On Apr 16, 1:50 pm, Alec Mihailovs wrote: > > Or even more simple -if A is the matrix consisting of all non-erased > columns, use A.transpose()*(A*A.transpose()).inverse() for decoding, > applying it to the non-erased part. And if A*A.transpose() is singular, the same thing, probab

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
On Apr 16, 1:39 pm, Alec Mihailovs wrote: > simple - find 10 other columns with non-zero determinant of the matrix > constructed from them, and use the inverse of this matrix for > decoding. Or even more simple -if A is the matrix consisting of all non-erased columns, use A.trans

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
e, the decoding seems to be rather simple - find 10 other columns with non-zero determinant of the matrix constructed from them, and use the inverse of this matrix for decoding. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to

[sage-support] Re: Coding theory:erasures

2010-04-16 Thread Alec Mihailovs
rrectly. Consider, for example, the last row with 5 zeros replaced with ones - if 5 erasures were allowed, then erasing them, we would get the last row, but erasing 4 remaining zeros in it instead, we would get the first row - so it couldn't be uniquely decoded. Alec -- To post to this group

[sage-support] Re: Vector solutions, a basic question.

2010-04-15 Thread Alec Mihailovs
For example, let the plane's equation is 2*x+3*y+4*z=5, and the line is given as =. Then the point of their intersection can be found in Sage as follows, var('t'); x,y,z=t-1,2*t-3,2-t plane=2*x+3*y+4*z-5 a=solve(plane,t)[0] [i.subs(a) for i in x,y,z] [1, 1, 0] Al

[sage-support] Re: Yamanouchi words

2010-04-14 Thread Alec Mihailovs
ut those words? The file devel/sage/sage/combinat/yamanouchi.py contains only that documentation string - nothing else. No classes or functions. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@goog

[sage-support] Re: Degree 120 splitting field

2010-04-14 Thread Alec Mihailovs
n by (x-b), then by (x-c). Now, the last step takes a long time, but it is not actually necessary - the polynomial is x^2 + (c + b + a)*x + c^2 + (b + a)*c + b^2 + a*b + a^2, so dividing it by (x-d), we get (x+a+b+c +d). Alec -- To post to this group, send email to sage-support@googlegrou

[sage-support] Re: Problems to susbstitute a function

2010-04-12 Thread Alec Mihailovs
ticular, in your example with h2, sage: function('f') sage: h2 = lambda t: sin(f(t)) sage: h2(x)=h2(x) sage: h2.diff() x |--> cos(f(x))*D[0](f)(x) Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-sup

[sage-support] Re: problems with sage and brian simulator

2010-04-12 Thread Alec Mihailovs
On Apr 12, 4:03 pm, Alec Mihailovs wrote: > Also, the show command should be replaced with savefig - something > like in the following example, That could be also done without changing the examples on the brian's webpage, by redefining show() as def show(): from matplotlib.py

[sage-support] Re: problems with sage and brian simulator

2010-04-12 Thread Alec Mihailovs
On Apr 12, 4:32 pm, William Stein wrote: > > You can also just do > > Action -> Restart Worksheet Thank you! Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com F

[sage-support] Re: problems with sage and brian simulator

2010-04-12 Thread Alec Mihailovs
By the way, the easy way to install brian package in Sage from the notebook is to execute the following cell, %sh easy_install brian Then save & quit the worksheet and open it again. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this g

[sage-support] Re: problems with sage and brian simulator

2010-04-12 Thread Alec Mihailovs
On Apr 12, 3:30 pm, Alec Mihailovs wrote: > > import brian_no_units > from brian import * Also, the show command should be replaced with savefig - something like in the following example, from matplotlib.pyplot import savefig eqs = ''' dv/dt = (ge+gi-(v+49*mV))/(20*ms) :

[sage-support] Re: problems with sage and brian simulator

2010-04-12 Thread Alec Mihailovs
om brian import * Then 1*mV 0.001 Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://w

[sage-support] Re: multidimensional numerical integration

2010-04-12 Thread Alec Mihailovs
35001e-13) Check it, integral(integral(x+2*y,y,1-x,1+x),x,1,2) 32/3 Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visi

[sage-support] Re: Question about numerical linear algebra.

2010-04-11 Thread Alec Mihailovs
tation for every method is available by entering A.(name of the method)? , for example, A.SVD? and Tab, and the source is available as A.SVD?? and then Tab. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...

[sage-support] Re: (1 - i) == -(i + 1) ??

2010-04-11 Thread Alec Mihailovs
information.        | > -- > sage: (1-I)/(sqrt(2)-1) > -(I + 1)/(sqrt(2) - 1) In Sage 4.3.4 (in Windows) that I am using, (1-I)/(sqrt(2)-1) (-I + 1)/(sqrt(2) - 1) Alec -- To post to this group, send email to sage-suppor

[sage-support] Re: Question about numerical linear algebra.

2010-04-11 Thread Alec Mihailovs
to Vh, as it supposed to be according to the help pages - it should be U*Sig*V.conjugate().transpose() == A Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visi

[sage-support] Re: points on elliptic curves mod N

2010-04-10 Thread Alec Mihailovs
2 does not exist Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.o

[sage-support] Re: Problems to susbstitute a function

2010-04-10 Thread Alec Mihailovs
On Apr 10, 10:37 am, Eckhard Kosin wrote: > > sage: sin.diff() > --- > AttributeError ... There is an interesting way to fix that, sin(x)=sin(x) sin.diff() x |--> cos(x) Alec -- To pos

[sage-support] Re: Arrangements().list() nice format

2010-04-10 Thread Alec Mihailovs
27;, 'eva']) (['eva', 'adam', 'otto', 'emil'], ['eva', 'adam', 'emil', 'otto']) (['eva', 'otto', 'adam', 'emil'], ['eva', 'otto', 'emil', &#

[sage-support] Re: Gap packages

2010-04-10 Thread Alec Mihailovs
lex.SimplicialHomology() [ [ 0 ], [ 1 ] ] gap("SimplicialHomology(SCNot2ConnectedGraphs(7),9)") [ [ 120 ] ] Seems to be working OK. Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+

[sage-support] Re: Arrangements().list() nice format

2010-04-10 Thread Alec Mihailovs
'adam'] ['eva', 'emil', 'adam', 'otto'] ['eva', 'emil', 'otto', 'adam'] ['otto', 'adam', 'eva', 'emil'] ['otto', 'adam', 'emil', 'eva'

[sage-support] Re: Problems to susbstitute a function

2010-04-10 Thread Alec Mihailovs
; x |--> cos(x) > > but > > sage: sin.diff() > --- > AttributeError ... Because S is a symbolic function, and sin is not. type(sin) The word 'symbolic' is missing. That, perhaps, can be

[sage-support] Re: Invoking Lisp from within Sage

2010-04-09 Thread Alec Mihailovs
On Apr 10, 1:32 am, Adam Getchell wrote: > Hi all, > > I realize this maybe a bit of an insane question, but I'm looking for a way > to use ecl within sage besides: > > ./sage -ecl For example, lisp("(def x 1)(defun f (x)(+ x 1))(f 2)") 3 Alec Mihailovs --

[sage-support] Re: Eigenvectors of a matrix

2010-04-09 Thread Alec Mihailovs
)).eigenvectors_right() Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org

[sage-support] Re: Inverses of Large Sparse Matrices

2010-04-09 Thread Alec Mihailovs
On Apr 9, 8:59 pm, William Stein wrote: > A 5000x5000 matrix just isn't really that big, IMHO... Actially, thinking about that, who knows what size its elements could be, if they are rational... They may be really big. Alec -- To post to this group, send email to sage

[sage-support] Re: Inverses of Large Sparse Matrices

2010-04-09 Thread Alec Mihailovs
irtual memory located on disk is used - it is faster to write on disk directly. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.goo

[sage-support] Re: Inverses of Large Sparse Matrices

2010-04-09 Thread Alec Mihailovs
on the disk incrementally, appending it to a file - that may be slower, but removes the restriction on the RAM size (but one would need to have enough space on the hard drive instead) - I am not volunteering though. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscr

[sage-support] Re: Problems to susbstitute a function

2010-04-09 Thread Alec Mihailovs
sqrt(R_b(t)^2-R_u^2)) 2*sin(psi(t))*cos(psi(t))*D[0](R_b)(t) - sin(psi(t))*R_b(t)*D[0](R_b)(t)/sqrt(-R_u^2 + R_b(t)^2) - (2*sin(psi(t))^2*R_b(t) - R_b(t))*D[0](psi)(t) Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this gro

[sage-support] Re: standard deviations in sage

2010-04-07 Thread Alec Mihailovs
)%RDF(5) Traceback (click to the left of this block for traceback) ... TypeError: unsupported operand type(s) for %: 'sage.rings.real_double.RealDoubleElement' and 'sage.rings.real_double.RealDoubleElement' Alec -- To post to this group, send email to sage-support@googlegr

[sage-support] Re: Can't solve equation with square roots...

2010-04-06 Thread Alec Mihailovs
=20, that the 2 triangles are similar. The left triangle has height 5, and the right one has hight 10, twice greater, so their sides on the x-axis also have the same ratio, 1:2, i.e. the length of the horizontal side of the left triangle is 1/3 of the distance from 0 to 20, i.e. 20/3. Alec

[sage-support] Re: Difference between sage and pyhton calculations

2010-04-05 Thread Alec Mihailovs
s 2 and 2.%3. is -1.00, also 2%3. is 2 and 2.%3 is -1.00. Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups

[sage-support] Re: Difference between sage and pyhton calculations

2010-04-05 Thread Alec Mihailovs
Or use Python floats, sage: float(6e-6) % float(10e-6) 6.0002e-06 In both cases the result is a Python float. To make it a Sage's real number, one can use RR, sage: RR(_) 6.00e-6 Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To un

[sage-support] Re: Interface to openssl hash functions?

2010-04-05 Thread Alec Mihailovs
d.p1 sage -f python-2.6.4.p7 Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-suppo

[sage-support] Re: nested "for" statements?

2010-04-05 Thread Alec Mihailovs
On Apr 5, 8:26 pm, Dana Ernst wrote: > Wow, thanks for all the feedback.  I learned a lot today! > > On Apr 5, 2010, at 5:31 PM, Alec Mihailovs wrote: > > u=[n for n in U if gcd(n,40)==1] By the way, one has to be careful using some of the methods of U. In partic

[sage-support] Re: nested "for" statements?

2010-04-05 Thread Alec Mihailovs
On Apr 5, 5:11 pm, Alec Mihailovs wrote: > u=filter(lambda n: gcd(n,40)==1, U) Also, that can be done as u=[n for n in U if gcd(n,40)==1] Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsub

[sage-support] Re: nested "for" statements?

2010-04-05 Thread Alec Mihailovs
On Apr 5, 5:11 pm, Alec Mihailovs wrote: > html.table([r]+[map(lambda j: k^j, r) for k in u], header=True) Or even more simple, html.table([r]+[[k^j for j in r] for k in u],header=True) Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this gr

[sage-support] Re: nested "for" statements?

2010-04-05 Thread Alec Mihailovs
r=True) Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org To unsub

[sage-support] Re: derivative instance

2010-04-05 Thread Alec Mihailovs
A: Yes. > Q: Are you sure? >> A: Because it reverses the logical flow of conversation. >>> Q: Why is top posting annoying in email? -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For mo

[sage-support] Re: Print only outputs?

2010-04-03 Thread Alec Mihailovs
opening it in a secure mode, i.e. with switch -Z if I remember correctly, which not that many people use.) Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options

[sage-support] Re: derivative instance

2010-04-03 Thread Alec Mihailovs
ives false, assuming f is not a derivative itself. One can do the following, for example, def is_diff(expr): try: return isinstance(expr.operator(), sage.symbolic.operators.FDerivativeOperator) except (AttributeError): return False Alec Mihailovs -- To post to this group, se

[sage-support] Re: List of "standard functions"?

2010-04-03 Thread Alec Mihailovs
On Apr 3, 11:55 am, Simon King wrote: >   sage: F = [f for f in globals().values() if callable(f)] I would write it as filter(callable, globals().values()) Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-supp

[sage-support] Re: Print only outputs?

2010-04-02 Thread Alec Mihailovs
it works with the current notebook. >   It doesn't affect the Print link, though. If one executes the following in a cell: html('div.cell_input_print {display: none}') and then clicks Print (at the lhs of the Worksheet button at the top), all input cells will be hidden in the hew

[sage-support] Re: Is there an efficient method of producing indexed variables?

2010-04-02 Thread Alec Mihailovs
2]] Vector(3, symbol=b); [b[1]] [] [b[2]] [] [b[3]] Alec -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-su

[sage-support] Re: Seeking Cython-101 lesson

2010-03-22 Thread Alec Mihailovs
t paid at least $200/hour or more, and will try to avoid any "interesting conversation" meaning that we have to listen to that, at any cost. Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sa

[sage-support] Re: Seeking Cython-101 lesson

2010-03-21 Thread Alec Mihailovs
eshmen, and 101 is usually a general education course taken by non-majors. Say Biology 101 means a general education biology course for non-biologists, and Math 101 is a general education course for non-mathematicians, as well as for people not taking any other math courses. Alec Mihailovs -- To p

[sage-support] Re: sage on slitaz linux?

2010-03-20 Thread Alec Mihailovs
, I guess. It is rather small, contains Sage 4.3.1 on Puppy Linux, and is available from http://mira.sunsite.utk.edu/sagemath/livecd/index.html as well on other mirrors. Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email

[sage-support] Re: atan2 throws "divide by zero"

2010-03-19 Thread Alec Mihailovs
Also, simplification might be improved. In particular, sage: atan(sin(1)/cos(1)).full_simplify() arctan(sin(1)/cos(1)) while atan(tan(1)).simplify() 1 Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage

[sage-support] Re: help with chanin dirivatives

2010-03-18 Thread Alec Mihailovs
q=diff(eta(x,t),x)==k(bx,bt) f(eta(x,t),bx,bt).diff(x,2).expand().subs(eq).subs(diff(eq,x)) delta^2*D[1, 1](f)(eta(x, t), delta*x, delta^2*t) + 2*delta*k(delta*x, delta^2*t)*D[0, 1](f)(eta(x, t), delta*x, delta^2*t) + delta*D[0](f)(eta(x, t), delta*x, delta^2*t)*D[0](k)(delta*x, delta^2*

[sage-support] Re: atan2 throws "divide by zero"

2010-03-18 Thread Alec Mihailovs
On Mar 17, 8:38 pm, G B wrote: > atan2(-pi,0) -->  RuntimeError: power::eval(): division by zero > > Any ideas how to get around this? A simple workaround is sage: from sympy import atan2 sage: atan2(-pi,0) -pi/2 Alec Mihaiovs -- To post to this group, send email to

[sage-support] Re: atan2 throws "divide by zero"

2010-03-17 Thread Alec Mihailovs
le it looks correct, it is not the answer that n should give. Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com

[sage-support] Re: Entering matrices whose elements are functions of row/column indices?

2010-03-17 Thread Alec Mihailovs
On Mar 16, 10:47 pm, Alec Mihailovs wrote: > By the way, vector doesn't work with integer numpy arrays, > > from numpy import array > vector(array([1,2])) > > Traceback (most recent call last): That shouldn't be too hard to implement - in particular, the following wor

[sage-support] Re: Syntax highlight in SAGE Calculator?

2010-03-17 Thread Alec Mihailovs
vector? vector?? Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org

[sage-support] Re: Entering matrices whose elements are functions of row/column indices?

2010-03-16 Thread Alec Mihailovs
e "", line 1, in File "free_module_element.pyx", line 309, in sage.modules.free_module_element.vector (sage/modules/ free_module_element.c:3172) TypeError: unsupported operand type(s) for ** or pow(): 'NoneType' and 'int' Alec Mihailovs -- To post to this group, sen

[sage-support] Re: Entering matrices whose elements are functions of row/column indices?

2010-03-16 Thread Alec Mihailovs
On Mar 16, 3:57 pm, Mike Hansen wrote: > On Tue, Mar 16, 2010 at 12:38 PM, Alec Mihailovs > > wrote: > > Will it work as matrix(5,f) as well, or only as matrix(5,5,f) ? > > It now works with both. > > --Mike That's great! What about vector(3,f), f.vector(3),

[sage-support] Re: Entering matrices whose elements are functions of row/column indices?

2010-03-16 Thread Alec Mihailovs
Mike Will it work as matrix(5,f) as well, or only as matrix(5,5,f) ? Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://grou

[sage-support] Re: Entering matrices whose elements are functions of row/column indices?

2010-03-15 Thread Alec Mihailovs
On Mar 15, 3:48 am, Alec Mihailovs wrote: > That seems to be a bug introduced in Sage, not in numpy, > > sage: fromfunction(lambda i,j:i-j+3,(1,1),dtype=int) > > array([[3]]) > > sage: matrix(_) > > [0 0 0] > [0 0 0] > [0 0 0] Yes, and that happens with other nu

[sage-support] Re: BoxPlot in Sage: how do we plot?

2010-03-15 Thread Alec Mihailovs
seems to be built without the png support in Sage. Here is a simple working example of a boxplot using matplotlib: from pylab import boxplot,savefig b=boxplot(range(10)) savefig("sage1.png") The original data list can be substituted in place of range(10) above. Alec Mihailovs -- To po

[sage-support] Re: Could someone help me.....

2010-03-15 Thread Alec Mihailovs
py.random.random_integers.html ) from numpy.random import random_integers dsums = random_integers(1, 6, 1) + random_integers(1, 6, 1) from matplotlib import pyplot counts = pyplot.hist(dsums, 11)[0]; counts array([ 279, 569, 845, 1072, 1398, 1683, 1355, 1106, 842, 576, 275]) pyplot.s

[sage-support] Re: Entering matrices whose elements are functions of row/column indices?

2010-03-15 Thread Alec Mihailovs
On Mar 15, 3:39 am, Alec Mihailovs wrote: > The second problem is that something strange happens for 1x1 matrices, > with any dtype, > > matrix(fromfunction(lambda i,j:i-j,(1,1),dtype=int)) > > [] > > matrix(fromfunction(lambda i,j:i-j+3,(1,1),dtype=int)) > > [0

[sage-support] Re: Entering matrices whose elements are functions of row/column indices?

2010-03-15 Thread Alec Mihailovs
On Mar 14, 6:15 pm, Alec Mihailovs wrote: > That could be also done as > > from numpy import fromfunction > > matrix(fromfunction(lambda i,j:i-j, (6,6), dtype=int)) > > [ 0 -1 -2 -3 -4 -5] > [ 1  0 -1 -2 -3 -4] > [ 2  1  0 -1 -2 -3] > [ 3  2  1  0 -1 -2] > [ 4  3

[sage-support] Re: Entering matrices whose elements are functions of row/column indices?

2010-03-14 Thread Alec Mihailovs
Integer Ring Alec Mihailovs -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org

[sage-support] Re: shift-return, shift-enter, click-mouse-on-evaluate. on evaluating a cell in a notebook

2009-04-21 Thread Alec Mihailovs
In Windows, such things, as well as any other key (and mouse) remapping can be done using an autohotkey script, http://www.autohotkey.com/ Alec Mihailovs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe

[sage-support] Re: plotting miltigraphs with pos

2009-04-10 Thread Alec
Also, show3d doesn't work with that graph, G.show3d() ZeroDivisionError: float division Alec Mihailovs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-uns

[sage-support] Re: Jmol causing Firefox crash

2009-04-10 Thread Alec
-plugin package. After that, I removed Jaunty from the software sources. Hopefully, soon Jaunty will be available for upgrading, and that will be even easier to do. Alec Mihailovs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To

[sage-support] Re: plotting miltigraphs with pos

2009-04-10 Thread Alec
e noticed that the degrees are printed not in the order of vertices in this graph, G.degree() [6, 6, 6, 3, 5] [i for i in G.degree_iterator()] [6, 6, 6, 3, 5] while it should be as in [G.degree(i) for i in G.vertices()] [6, 6, 6, 5, 3]

[sage-support] plotting miltigraphs with pos

2009-04-08 Thread Alec
I just tried to plot a multigraph with setting positions of vertices, G=Graph({'a':['a','b','b','b','e'],'b':['c','d','e'],'c': ['c','d','d','d&

[sage-support] Re: Maple in Sage on Mac

2009-01-17 Thread Alec Mihailovs
lpful, http://www.mapleprimes.com/blog/paul/usingopenmaplec Alec Mihailovs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more op

[sage-support] Re: NVIDIA Tesla

2008-11-23 Thread Alec Mihailovs
Thank you! > http://groups.google.com/group/mpir-devel/t/df88735e6d4e678c I should search the group before posting. Alec --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email

[sage-support] NVIDIA Tesla

2008-11-23 Thread Alec Mihailovs
Is it possible to run SAGE on NVIDIA Tesla (with 4 teraflops)? http://www.nvidia.com/object/personal_supercomputing.html http://tech.slashdot.org/article.pl?sid=08/11/23/068234&from=rss Alec --~--~-~--~~~---~--~~ To post to this group, send email to

[sage-support] Re: square of an inequality

2008-08-25 Thread Alec Mihailovs
or example, In[1]:= Reduce[x + y <= 3 && x >= 0 && y >= 0, {x, y}, Integers] Out[1]= (x == 0 && y == 0) || (x == 0 && y == 1) || (x == 0 && y == 2) || (x == 0 && y == 3) || (x == 1 && y == 0) || (x == 1 && y ==

[sage-support] Re: Hypergeometric sum

2008-08-23 Thread Alec Mihailovs
Ondrej and William, Thank you, it's nice to know about that. Alec --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this gro

[sage-support] Re: Hypergeometric sum

2008-08-23 Thread Alec Mihailovs
> What is your point? Actually, I meant that more about Maplesoft. Sage support is usually good in this group, just wasn't very good in this thread. Alec --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To uns

[sage-support] Re: Hypergeometric sum

2008-08-23 Thread Alec Mihailovs
very nicely, the problem was addressed, and a workaround was suggested). And the winner is - Wolfram Research. Alec --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED

[sage-support] Re: square of an inequality

2008-08-23 Thread Alec Mihailovs
E do that? I meant to solve. The answer should be either a list of points, {x=0,y=0}, {x=0,y=1} etc. - or their convex hull. That gives feasible points for integer (linear or convex) programming problems. I posted about that on Mapleprimes (Kinde

[sage-support] Re: Hypergeometric sum

2008-08-21 Thread Alec Mihailovs
ica/browse_thread/thread/b84ac9c2a48f1bc1?hl=en# Alec --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.

[sage-support] Re: Hypergeometric sum

2008-08-21 Thread Alec Mihailovs
be equal to the expression given by Mathematica if n was not a positive integer. Another form of the correct answer is (2 x)^(n+1/2) E^x BesselK[n+1/2,x] n!/(2 n)!/Sqrt[Pi] Alec --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com

[sage-support] Hypergeometric sum

2008-08-21 Thread Alec Mihailovs
Both Maple and Mathematica give wrong answers to the following sum, Sum[Binomial[n, k]/Binomial[2 n, k]/k! (2 x)^k, {k, 0, n}] I tried to find a way to calculate it in SAGE, but couldn't find symbolic sums in the documentation. Is Maxima supposed to be used directly?

[sage-support] Re: square of an inequality

2008-08-20 Thread Alec Mihailovs
should be a part of the standard test, I think. Maple can't do that. I didn't try Mathematica, but I believe that it is capable of that. Can SAGE do that? Alec --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To

  1   2   >