MaplePrimes Questions

I have three columns of data (real numbers) in Excel that have about one thousand rows and I wish to plot their relationships in Maple. The source data simply comprises three columns of numbers and so, are not ordered in triples, [x,y,z]. Does anyone know of a simple routine available that allows me to import the raw data and output a 3-D wire / solid plot in Maple?

Thank you! 

a1:=odeplot(p,[t,x(t)],0..5,color=black,thickness=2):

a2:=plot((2/125)*t^3-(3/25)*t^2+t+1,t=0..5,color=red,style=point,thickness=4):
display(a1,a2,labels=[x1,x2]);
Error, (in plots:-display) expecting plot structure but received: a1
 

I have an explicit mathematical expresssion in the form

m = f[1](x,y)*A + f[2](x,y)*B + f[3](x,y)*C

I wonder if it is possible to use the collect function on the vector z = [A,B,C]^T and write the resulting expression as the dot product

m = f^T.z 

Basically, I am interested in obtaining the vector f = [f[1](x,y),f[2](x,y),f[3](x,y)]^T as a Maple variable in closed-form.

 

Thank you in advance for your help.

 

Marco

I'm referancing the maple help page: https://www.maplesoft.com/support/help/maple/view.aspx?path=Groebner in my thesis, but I am not sure what fields to include, and what to put in them- particularly for the year.

I'm a newbie.I am trying to solve a differential equation symbolically.But I'm having  a\some problem.

First 3 lines works fine....but 4th line gives no output...

Secondly How I'm supposed to simplify the output of 3rd line

 

1.>ode := P*(diff(y(x), x, x))+diff( K(x)*(diff(y(x), x, x)), x, x ) = 0

2.>ics := y(0) = 0, y(L) = 0, ((D@@2)(y))(0) = 0, ((D@@2)(y))(L) = 0

3.>dsolve(ode, y(x))

3rd lines output--y(x) = _C2*x+_C1+Int(DESol({diff(_Y(x), x, x)+2*(diff(K(x), x))*(diff(_Y(x), x))/K(x)+(P+diff(K(x), x, x))*_Y(x)/K(x)}, {_Y(x)}), [x, x]) 

 

4.>dsolve({ics,ode},y(x))

5.>

4th lines gives no output.After hitting "Enter",I get a new line in my Worksheet(line 5).Why..can anyone explain please?

 

I want to match   a*x^b only if "a" and "b" are both not the symbol "y"

Therefore 5*x^y should fail to match, but this below matches.

restart;
expr:=5*x^y;
if patmatch(expr,conditional(a::anything*x^(b::anything),(a<>y and b<>y)),'la') then
   print("it matched, why??");
   map(assign,la):
   print("a = ",a, " b=",b);
fi;

it gives

                      "it matched, why??"
                      "a = ", 5, " b=", y

What did I do wrong? is my conditional wrong above?

 

I was working through some example problems and came across this statement.  Actually comparing Sympy results to Maple. (I prefer Maple).

"Calculate the volume integral of f(...) over the sphere of radius r"

I can solve the problem, but, got hung up on the exact meaning of the problem statement.

Tom Dean

restart
with(VectorCalculus):
with(LinearAlgebra):

## http://www.acme.byu.edu/wp-content/uploads/2017/08/Sympy.pdf
## Problem 7

f := proc(x, y, z)
    (x^2 + y^2 + z^2)^2
end proc;

(M, d) := Jacobian([rho*sin(phi)*cos(theta),
                    rho*sin(phi)*sin(theta),
                    rho*cos(phi)],
                   [rho, phi, theta],
                   'determinant' );
abs(d);
simplify(%);

eqn := f(rho*sin(phi)*cos(theta), rho*sin(phi)*sin(theta), rho*cos(phi));
eqn := eqn * abs(simplify(d));

soln := int(eqn,[rho=0..r,theta=0..2*Pi,phi=0..Pi]);
subs(r=3, soln);

 

How I can convert Root Of to conventional form?

Thanks

root_of.mw
 

Q1 := x*(x-delta)*(1-x)*(x+y)-alpha*x*y; Q2 := beta*x*y-Zeta*y*(x+y); SOLL := solve({Q1, Q2}, {x, y})

{x = 0, y = 0}, {x = 1, y = 0}, {x = delta, y = 0}, {x = RootOf(beta*_Z^2+(-beta*delta-beta)*_Z-alpha*Zeta+alpha*beta+delta*beta), y = -RootOf(beta*_Z^2+(-beta*delta-beta)*_Z-alpha*Zeta+alpha*beta+delta*beta)*(Zeta-beta)/Zeta}

(1)

``


 

Download root_of.mw

 

update_standings_ELO_draw.mw


This latest question involves ELO, which is used to rate chess players (and other sports).

The difference between the previous (win-loss) system and ELO is the allowance for a third option, a draw. Carl coded the previous rating system, and I have tried to alter it.

See code for ELO procedure.

Adding a D to the third field, gets the result I want.

ELO(Record(rating= 2400),Record(rating= 2000),D)

Adding anything other than D to the third field means the first player defeats the second.

My preference would be to leave nothing in this third field to signify this (but I don't know if this is possible). 

The main problem is getting the bulk processing to work (for a series of matches...)

See Update procedure.

I have this flawed notation:

Games:= [[Tom_C, David_J, D],[David_J, Tom_C]]: #Tom_C drew David_J then David_J defeats Tom_C.

Hope someone can help.

I am a newbie. I m having some problems in solving generalised eigenvalue problem for symbolic matrices. Help me to discover the problems with this code.

Here is my code:-

>with(LinearAlgebra):

> M := Matrix(2, 2, [4*(I*E/L^3), 6*(I*E/L^3), 6*(I*E/L^3), 24*(I*E/L^3)])

> G :=Matrix(2,2,[2/(15)(-1/(L)),1/(10)(-1/(L)),1/(10)(-1/(L)),12/(5)(-1/(L))])

> Eigenvalues(M, G, output = 'Vector[row]')

 

Error:-

Error, (in LinearAlgebra:-Eigenvalues) invalid input: LinearAlgebra:-Eigenvalues expects its 2nd argument, C, to be of type Matrix(square) but received G
 

If decimal place is 1 such as 123.1 as input matrix and output matrix is eigenvector

is there a easy way to see changes of eigenvector in graph ?

because 3 x 3 matrix has 9 elements 

and eigenvector each vector may not fix in place,

how do mathematician observe the perturbation in this case?

would like to find which changes can cause different in eigenvectors

 

 Dear Sir,

I trying to merge two which are one is an actual plotted graph and another one zoomed from the same graph. I want to copy zoomed graph into an actual graph, can anyone help in this regard?

My actual graph is like as follows

and after zooming this graph I got as follows

Now I want to change (Insert the zoomed graph inside the actual graph where the text is inserted) by using these two graphs as below

Please help me to plot like this...

I am currently facing an issue as regards stability analysis and I will be glad if anyone can offer me some suggestions.  My aim is to have a Theorem in my research stating that the Disease free equilibrium is locally asymptotically stable if R0 (Basic reproduction number) < 1 and unstable if R0 > 1. The same applies to the endemic equilibrium.  I am trying to carry out a stability analysis using Linearization. However, due to the nature of my model, not all my eigenvalues are negative. Below is my Jacobian Matrix:

ROW 1                           ROW 2                       ROW 3                         ROW 4

Column 1:=              -psi-gamma-mu                       varphi                              omega                                   betaS

Column 2:=                 -psi                 - varphi - varpi - mu                    0                                  epsilon * beta  M

Column 3:=                   gamma                         varphi                             - omega - mu                                  xi * beta * V

Column 4:=                 0                                0                                0                       betaS + epsilon * beta  M + xi * beta * V -tau - kappa - mu

By using the Routhz Hurwitz criterion, I reduced the matrix block from 4 by 4 to 3 by 3 by using column 4 as co-factor, however, after this matrix reduction, I was still unable to get all terms relating to my basic reproduction number for me to compare and see if the conditions in Routh Criterion is satisfied using the condition R0 < 1 holds.

My Basic Reproduction Number is as follows:

(varphi+mu+varpi)*(mu+omega-mu*rho)-(mu+varphi+omega)*mu*sigma+ epsilon ((omega+mu+gamma)*sigma*mu+(mu+omega-mu*rho)*psi) + xi*((varphi+varpi+mu+psi)*mu*rho+(mu*sigma+gamma+psi)*varphi+(mu-mu*sigma+varpi)*gamma)/((varphi+mu)*(psi+gamma+omega+mu)+(omega+mu+gamma)*varpi+psi*omega)

Are there any other methods to  carry out local stability analysis? and Can local stability analysis be done using maple?

Thank you in advance

Hi everybody!

 

I want to make the substitution  Int(x^n u(x),x=-Infinity..Infinity) by m_n (m index n), the subs command doesn't seems to work...

Thank you!

I can calculate real number eigenvector with QR algorithm from rotation matrix

but this real number is just magnitude of complex number

what is the rotation matrix for calculating complex number?

what algorithm do maple use to calculate fraction eigenvector?

are complex decimal and complex fraction eigenvectors are two different algorithm?

First 529 530 531 532 533 534 535 Last Page 531 of 2154