Hi everyone:

I want to obtain only the vertical asymptotes this function in the interval [a, b] and i used the asymptotes command but i didn't get true answer, can everybody help me? a=0.3 and b= 0.36

I attached my code. 

tnx...

Download Asymptotes.mw

Where is the mistake?
 

I'm working with tensors using the DifferentialGeometry package (not yet checked out the Physics package since I plan on working with Lie algebras), and in particular I'm computing the Riemann tensor of a metric.

Raising the first index of CurvatureTensor(g), the Riemann tensor has the symmetries of a symmetric product of two 2-forms, and this is how I would like to present it, identifying terms like dw1dw2dw1dw2 with terms like dw1dw2dw2dw1 with the appropriate sign. Ideally I would arrive at an expression that is written in terms of (dw1^dw2) (dw1^dw2) (etc.) basis elements, provided that my tensor has these symmetries.

Is there some readily available functionality to do this, or a simple, standard method in Maple?

Dear Users! 

Hope you would be fine with everything. I want to define a matrix F of M+1 by M+1 order having element of the following form:

I derived the F[r,s] but confuse who to generate matrix now.

restart; M := 5; printlevel := 3; for r from 2 while r <= M+1 do for s while s <= r-1 do if type(r+s, odd) then F[r, s] := 2^(k+1)*sqrt((2*r-1)*(2*s-1)) end if end do end do

I have got the following function:

R = ((n1-n2)/(n1+n2))^2

I would like to plot R as a function of (n1/n2). How can I do this in maple?

Hi everybody:

How can i plot this function in maple? for x=0.2..0.4 and I is complex number(I^2=-1).

tnx.

f:=(x)->-(2.50000000000000000000000000000*(-2.424176758*10^30*x^5+5.982448001*10^29*x^3+(5.336360731*10^33*I)*exp((4.0405213687327619451756171962*I)/x)+(3.144879905*10^98*I)*exp((48.6318841588794957405477169504*I)/x)+(8.489643736*10^35*I)*exp((4.04052136873276194517561719622*I)/x)))*x^6/((1.009157478*10^101*I)*exp((48.6318841588794957405477169504*I)/x)*x^6-2.760182453*10^32*x^11-(7.651751413*10^37*I)*exp((4.0405213687327619451756171962*I)/x)*x^4+1.887831524*10^32*x^9-(8.567811466*10^36*I)*exp((4.04052136873276194517561719622*I)/x)*x^4-(1.207009672*10^100*I)*exp((48.6318841588794957405477169504*I)/x)*x^4+(6.538572988*10^29*I)*x^6+(1.246166933*10^31*I)*x^10-1.837651455*10^31*x^7+(1.504191760*10^38*I)*exp((4.04052136873276194517561719622*I)/x)*x^6+1.574891043*10^34*exp((4.0405213687327619451756171962*I)/x)*x^5-(7.612867062*10^30*I)*x^8+4.556143298*10^99*exp((48.6318841588794957405477169504*I)/x)*x^5-(2.488678360*10^105*I)*exp((52.6724055276122576857233341466*I)/x)+6.791123646*10^36*exp((4.04052136873276194517561719622*I)/x)*x^5-3.454612059*10^36*exp((4.0405213687327619451756171962*I)/x)*x^3-4.710342238*10^98*exp((48.6318841588794957405477169504*I)/x)*x^3+(3.488285966*10^35*I)*exp((4.0405213687327619451756171962*I)/x)*x^6+5.512268188*10^106*exp((52.6724055276122576857233341466*I)/x)*x)

Im trying to create a basis as similar to this document. I have created the 72 equations under root of unity method (aswell have 44 equations on a seperate document for the sum shidoku product way) but cant seem to generate a basis for it such as they did in the document with these equations. If anyone can help as i may be entering the Basis() wrong. If you need the 72 equations i can post them as well. Any tips that you have would be appreciated!

http://reu.dimacs.rutgers.edu/~ericlee/Shidoku.pdf

Hello everyone,

I have a project that uses several machines to make a computation. I wrote a package (called here "package") with procedures, "proc".  I wrote several .mpl files structured as:

read "package.mpl":

W:= Read("file1.mpl"):    #I wrote a procedure to read and assign to a variable. This is not the problem.

Export("output.mpl", proc(W));

The situation is as follows: in all but one of my machines, this works fine, and proc(W) is exported correctly. However, in this ugly duckling (which is running the GUI version of Maple18), maple prints as an output, in blue:

Export("output.mpl",...)

and in place of ... it gives the answer expected. But it does not export automatically. I can, of course, manually correct this, but this is not optimal. Moreover, if the output is too big, there will be display issues.

Help? Any light is appreciated, thanks.

Marcelo

PS: I do not expect an answer for the PS; but I also do not know how to make my package "official", in the sense that I always read the .mpl file with all the procedures. Is there a way to call it like the LinearAlgebra and so on and use the package:-proc syntax? 

I have a problem with this sequence, as defined : f(x)=1-1/x+2*trunc(1/x)

The sequence starts with the number x=50100 and with 100000 steps the aim of the sequence is to see how many numbers are between 0 and 1( and 1 and 2, 2 and 3, and 3 and 4). My idea is:

 f:=x->1-2*(1/x-trunc(1/x))+1/x;
       x:=50100;
       k:=0 : a:=0 : b:=0 : c:=0 : d:=0 : while k<100000 do x:=f(x): k:=k+1 : if(x>0 and x<1)   then(a=a+1)           elif(x>1 and x<2) then(b=b+1) elif(x<2 and x>3) then(c=c+1) elif (x>3 and x<4) then (d=d+1) end             if : end do: printf("%d %d %d %d", a,b,c,d)                                                                                                                                             But i have problems :(.  Does anyone see the mistake?

I want to solve this:
ode := diff(f(x), x, x)+(diff(f(x), x))/x-(1/2)*(1/x^2+3.6*10^17/(6.6))*sin(2*f(x)) = 0;

bcs := f(rho) = 0, f((4*10^(-8)-rho)^.5)

and I want to plot it for a range of rho between 10^-9..7*10^-8.
what should I do?

Due to the mechanistic process of our students and little creativity in analysis in schools and universities to be professionally trained is that STEM education appears (science, technology, engineering and mathematics) is a new model that is being considered in other countries and with very slow step in our city. In this work the methods with STEM will be visualized but using computational tools provided by Maplesoft which is a company that leads online education for adolescents and adults in the current market. In Spanish.

ECI_UNT_2018.pdf

ECI_UNT_2018.mw

Lenin Araujo Castillo

Ambassador of Maple

Hello,

I've a big amount of equations, i tried to solve them and get the answer as an array but for some reason it's not working, is solve the right one for this or what?

thank you

the 2 files aren't related, but both contain a system of equation i tried to solve but without avail

Advheat-ex53.mw

Advheat-ex54.mw
 

Download Advheat-ex54.mw
 

Download Advheat-ex53.mw

Hi,

I'm trying to help the daughter of one of my friends with some flux related exercises.

We have calculated this one by hand and then I wanted to repeat it in Maple to help her see how to use Maple.

My problem is that I cannot figure out how to represent the functions in Maple, so that the calculations fall out similar to what we have been doing by hand. Yes, we could just use the built in Flux operation, but that defeats the purpose here.

Please bear with me - my calculus is rusty and I only have access to Maple on my son's machine. Plus I have never used Maple before... yes, I'm that old ;-)

We want to calculate the flux of the vector field V through the surface given by r:

Then we can differentiate r wrt both u and v and calculate the normal vector for the surface:

The next step is then to do the double integral of:

This is - bar the evalVF call - almost identical to what you'd be doing by hand, so this is good.

Next is the flux out of the 3d object constrained by the square [0,1] and [-Pi/2, Pi/2] and the surface r.

The paramterisation of that is:

with w in [0,1].

At this point one applies Gauss' Divergence Theorem which means that in order to find the flux of V out of Omega we "just" have to calculate the spatial integral of the divergence of V over the body described by Omega.

But that can be calculated by integrating V(Omega(u,v,w)) times the Jacobi function.

The Jacobi function is the determinant of the output of Maple's Jacobian procedure:

In this case it is 1 - u^3.

Now comes the big problem: how do I get Div(V)(Omega(u,v,w)), ie, the divergence of V taken in the point Omega(u,v,w)?

If I do evalVF(V, Omega(u, v, w)) I get a column vector out that I cannot compute the divergence of using Divergence(%) since that only works on VectorFields.

And Divergence(V) is not a function so I cannot do

And if I change Omega to return a value in R^3 instead of a vector I cannot use Jacobian, so that stop even faster.

Defining two versions of Omega just to make Maple happy seems wrong.

I just want to write the tripple integral of Div(V)(Omega(u,v,w))*jacobi or something close to that.

Please let me know if I should provide the workbook or in some other way make all my calculations available.

Thanks in advance and sorry if my description is lacking,

Torben

Dear Users!

Hoped everyone fine with everything. I the following maple expression, I need a matrix A for each n. Like if I take k =1 I want A[1]; if I take k=2, I want A[1], A[2]; for k=3 I want A[1], A[2], A[3] and so on. A[i]'s is square matrix having order M-1 by M-1.

Further I want to generate a block matrix for k. Like for k=1 I want a block matrix as Vector(1, {(1) = A[1]}), for k=2 I want a block matrix as Matrix(2, 2, {(1, 1) = A[1], (1, 2) = 0, (2, 1) = 0, (2, 2) = A[2]}), for k =3 I want a block matrix as Matrix(3, 3, {(1, 1) = A[1], (1, 2) = 0, (1, 3) = 0, (2, 1) = 0, (2, 2) = A[2], (2, 3) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = A[3]}) and so on.

restart; alpha := 1;
k := 2; M := 3;
printlevel := 3;

for n while n <= 2^(k-1) do

for m from 0 while m <= M-1 do

for j from 0 while j <= M-1 do

Omega[m, j] := 2^((1/2)*k)*sqrt(GAMMA(j+1)*(j+alpha)*GAMMA(alpha)^2/(Pi*2^(1-2*alpha)*GAMMA(j+2*alpha)))*(sum((-1)^i*GAMMA(j-i+alpha)*2^(j-2*i)*(sum((1/2)*binomial(m, l)*(2*n-1)^(m-l)*(1+(-1)^(j-2*i+l))*GAMMA((1/2)*j-i+(1/2)*l+1/2)*GAMMA(alpha+1/2)/GAMMA(alpha+1+(1/2)*j-i+(1/2)*l), l = 0 .. m))/(GAMMA(alpha)*factorial(i)*factorial(j-2*i)), i = 0 .. floor((1/2)*j)))/2^(k*(m+1))

end do

end do;

A[n]:=???

end do;

I am waiting for your positive response.

Thanks
 

Hello i keep getting an error saying that maple expects an integer as it's argument but i has an integer so it is sort of the same thing like when a cat has a full bowl of biscuits and it is going mental at you to feed it, in every other circumstance i would'nt get this error and i even did an assume(k::integer) for the index of the series omega appears, still no. 

I am refering to the number theory package bigomega function btw