god.one's profile - activity

Hi,

you can use line breaks line breaks

phi(epsilon, Q, r) = 1/(4*pi* \
                    epsilon \
)*Q/r

show(phi(1,2,3))

Worked for me in the sage command line v8.2

Why are there elements in the top right and bottem left corner of the matrix?

Or you could try numpy which is roughly 50 times faster

import time
import numpy as np
a=QQ(-1)

elapsed_time=0
for i in range(100):
    start=time.time()
    test = a*np.ones([1000,1000])
    end=time.time()
    elapsed_time=elapsed_time+(end-start)

print(elapsed_time/100)

But keep in mind that the example

from itertools import repeat
matrix(QQ,1000,1000,repeat(a,1000000))

generates a different type than a numpy array does

type(matrix(QQ,1000,1000,repeat(a,1000000)))
<type 'sage.matrix.matrix_rational_dense.Matrix_rational_dense'>

type(a*np.ones([1000,1000]))
<type 'numpy.ndarray'>

Hi,

it is possible, just try

var('x')
((cos(x)-sin(x))^2).trig_simplify()

which will result in

-2*cos(x)*sin(x) + 1

since sin(x)^2+cos(x)^2=1 (I think you mean this formula with the trigonometric identity).

Hi, I want to plot two region_plots into one plot. The idea is if you use different colors for the inequalities you can see how the regions change. A minimal example:

var('x,y')
plot1=region_plot(x<y,(x,0,1),(y,0,1),incol='red')
plot2=region_plot(2*x<y,(x,0,1),(y,0,1),incol='blue')
show(plot2+plot1)

I do not know how to manage this, since zorder or opacity are not working for region_plot. Does anybody know how to make this work (by the way I use sagenb.org)? Thanks in advance

Hello everyone,

I wanted to fit some data and I found a strange behaviour of solve(). After evaluating cell 1,2 and 3 and further evaluations of cell 2,3,2,3,2,3,... the solution for kf changes and gets smaller and smaller. I use Sage 5.9 in a VM in Win7Prof

cell1

# data=[another list]
data=((0.004,34.915),(0.018,34.665),(0.040,34.254),(0.071,33.695),(0.159,32.146),(0.282,30.323))
p1=list_plot(data,plotjoined=True,xmin=0,xmax=1,ymin=20,ymax=35)

cell2

# modell for fitting
lg=20
lm=35
A=(lm-lg)/(1-exp(-k))
c=lm-A
var('kf')
testmodel(t)=c+A*exp(-kf*t)

best_fit_par=find_fit(data,testmodel(t),initial_guess=[0.5],parameter=[kf],variables=[t],solution_dict=True)
print 'best_fit_par'
print best_fit_par

cell3

k=best_fit_par[kf]
lg=20
lm=35
A=(lm-lg)/(1-exp(-k))
c=lm-A
p2=plot(A*exp(-k*t)+c,(t,0,1),xmin=0,xmax=1,ymin=20,ymax=35)
show(p1*p2)

The values for kf started at 0.49 and are wandering towards 0.35

Thanks in advance

Thanks, but in this case using different assume scenarios resulting in different results is strange (escpecially in the both is negative case)

You can try the following with substituting b-a with c

var('a,b,c,A,B,C')
eq1 = (sqrt((c)^2 + (B-A)^2) == A+B)
assume(A>0)
assume(B>0)
solve(eq1,c)

which results in

[c == -2*sqrt(A*B), c == 2*sqrt(A*B)]

But can anybody explain me the results of the computation without the assume statements? For A,B both beeing negative the solver does not do anything and for A being bigger and B being smaller then zero, sage throws an error that maxima needs further inputs (the same as when you're not using assume)

Hi, could you please tell us which version of sage you are using? Going to sagenb.org

var('x y')
sol=solve(x + y, [x, y], solution_dict=True)
sol

gives

[{x: -y}]

which still surprises me since the system of equation is more than incomplete.

Try imag, imag_part or imaginary, they should do the trick. And remember to use help or ? often to find out what a function does.

Hi,

sage is correct saying that a negative number cannot be raised to a fractional power IF you want to plot it. For example

f(-2)=-1/2*sqrt(-2)
f(2)=sqrt(2)

The first one is imaginary, the other one real, so you cannot plot both in the same 2D plot. You could use the abs() function to avoid this problem. Another thing is

sage: f(x).limit(x=0,dir='+')
0
sage: f(x).limit(x=0,dir='-')
+Infinity

so that your plot for either f or g will have a huge peak which will suppress the rest of the plot. Try the following:

sage: f1(x)=x^(1/x)
sage: f2(x)=(-x)^(1/x)
sage: p1=plot(f1(x),0,5)
sage: p2=plot(f2(x),-5,-0.5)
sage: show(p1+p2)

I'm sorry, this forum is not here for solving your homework. If you have any non-working code which you could provide, anybody here is happy to help you.

And since this is a special use-sage-task, propably your instructor is lurking around here, too :)

Hello everybody,

for future calculations I need a CAS. Up to now, I worked with a free mathematica student licence. But now I need a officially bought version (which I can't afford). I'm now trying to work with Sage but my head of department is no fan of open source projects and thinks, among other things, that it is probably buggy and not as powerful as said mathematica. Could you please tell me how "stable" and "correct" the sage calculations are and give me further information on this topic (I already searched the forum for some information) so that I can present him/them why we should use Sage and nothing like mathematica/mathlab/...

Please tell me if you need more information or clarification.

Thank you all in advance for your effort.

Daniel

Hi, here is my piece of code

var('x,m')
assume(m, 'integer')
y = function('y', x)
desolve(diff(y,x,2) + 2*x*diff(y,x) - 4*m*y, y,contrib_ode=true,ivar=x)

which yields

[y(x) == k1*kummer_m(-m, 1/2, -x^2) + k2*kummer_u(-m, 1/2, -x^2)]

Internet(Wikipedia) sais that the solutions to this differential equation are the Kummer functions or Confluent_hypergeometric_functions. I now want to know how I can use these functions in sage since they must be defined in maxima.

Thanks in advance

Don't know if this helps but

x=var('x')
def P(f,m,n):
    return derivative(f,x,m+n)
f=x^20
P(f,1,2)

works fine with sagenb.org

Hi,

I got another problem relating to this question http://ask.sagemath.org/question/8137... and many other of that kind (sorry for opening another question). If I use

sum(-1**(i+1)/(2*i-1)**3*(sin(pi*(2*i-1)/2*(r1+r2))*sin(pi*(2*i-1)/2*(r1-r2))),i,1,oo)

in the shell, I get

TypeError: unable to make sense of Maxima expression 'imagpart(f[4,3]([1/2,1/2,1/2,1],[3/2,3/2,3/2],e^-(2*I*pi*r1)))' in Sage

I read in http://trac.sagemath.org/ticket/9825 that the problem concerns the temporary variables. But is there a possibility to reformulate the problem to get to the solution? Btw, I'm using Sage v6.2 on a LinuxMint machine.

On sagenb.org using var('x,y'); solve([y==sin(x),y+4*sin(x)==5],x,y) I get [y == sin(x), y + 4*sin(x) == 5] as a solution

On sagenb.org using

var('x,y')

solve([y==sin(x),y+4*sin(x)==5],x,y)

Sehr schoen. I will try this. If there is no easy way to turn .sws files into .sage files, the old copy&paste will have to do it for now.

Hello everyone.

I used sagenb.org for a very long time and now I wanted to migrate all my files to my PC. Is there any possibilty to load the .sws files into the sage shell? Searching the web did not bring any answers.

Another question: what is the best substitution to a sws in the shell? Should I write my sage code in another editor and then run it in the sage shell (eg. write all def in a py file, load it anew in every session,...)?

Thanks in advance

Hi, I got a question unrelated to this one but to your code. If I stupidly copy and paste it in a sagenb.org notebook, it throws different errors like 'invalid syntax'. The version that works for me is timeit('[randint(0,10) for i in xrange(1000)]') timeit('[randrange(0,11) for i in xrange(1000)]') import numpy timeit('[numpy.random.random_integers(0,10) for i in xrange(1000)]') timeit('numpy.random.random_integers(0,10,1000)') Can someone explain this to me or give me a hint for which keywords I have to search? Thanks in advance.