MaplePrimes Questions

Hi there,

this may be a common task but did not find any helpful hint nor an answer: is there any way to combine to arrays A and B into another one element-wise, i.e.:
C = [[A[1],B[1]], [A[2],B[2]],...,[A[N],B[N]]]
without doing a for loop?


Thanks,
jon

Recently, in one of the old book on programming came across the following problem: to place on the chessboard 5 queens so that each free field was attacked by at least one queen. This problem is called the problem of the dominant queens. I have not seen the implementation of this task in Maple. Naturally to solve this problem for an arbitrary board N by N. I have 2 variants to solve the problem, but I am not going to to present them yet, so that everyone can enjoy the independent decision. Especially the interesting case is the case of the board 6 by 6, when the solution is unique (certainly up to symmetry).

Of course, it is interesting for each board  N by N to find the minimum number of queens that satisfy the above condition. It seems that for arbitrary board  N by N the exact value of this number is not known. I do not know any other way of solving the problem as a brute force method.

Hye,

Have a good day Sir.

can u help me how to solve my problem using the Implicit Crank Nicolson Finite different Method. I have attached the problem and what I have done in Word file. I have to modify my Equations in order to get the graphs (that is wrong actually). Can you help how to get the graphs by not change the sign of my equations?:( 

 

Thank you in advance.

God Bless you.

 

 

1. if parametric equations are
x0 = a^4, x1=a^3*b, x2=a^2*b^2, x3=a*b^3, x4=b^4

how to find ideal below

e12 = x0*x2-x1^2
e13 = x0*x3-x1*x2
e14 = x0*x4-x1*x3
e23 = x1*x3-x2^2
e24 = x1*x4-x2*x3
e34 = x2*x4-x3^2

2.
if let y_i = x_i/x0 for i = 1..4
how to find above parametrically by below equations
y1 = b
y2 = b^2
y3 = b^3
y4 = b^4

Dear any body,

I have tried several models/documents (with diagrams) using maplesim v7 32/64bit, but the simulation can not finish and in  equation generation a message appear as follow "protype must be a sructured record type or a record, but receipt squeezed".Same mesaages also appear if I try all examples exist in maplesim. What a big problem this?

Please help me..thanks in advance

I am trying to learn the <> notation to enter matrices and vectors. But I find this page very confusing

http://www.maplesoft.com/support/help/maple/view.aspx?path=examples%2FLA_Syntax_Shortcuts

it says:

but we see clearly the vertical bars are used to separate columns.

Isn't a column the thing that goes from the top to bottom and not from left to right in Maple LinearAlgebra?

 

Hi

My equations are:

diff(Th(z, t), t) = 7.1428*(diff(Th(z, t), z))-1397941.885*(279-Tw(z, t))-0.2160487e-1*(diff(Th(z, t), z, z))

diff(Tc(z, t), t) = -7.1428*(diff(Tc(z, t), z))+1298990.852*(Tw(z, t)-291)+0.189366e-1*(diff(Tc(z, t), z, z))

diff(Tw(z, t), t) = 3.3024901*(Th(z, t)-2*Tw(z, t)+Tc(z, t))+8.0029*10^(-4)*(diff(Tw(z, t), z, z))

and boundry conditions are:

Th(z, 0) = 296, (D[1](Th))(0, t) = 0, Th(1, t) = 296

Tc(z, 0) = 275, (D[1](Tc))(1, t) = 0, Tc(0, t) = 275

Tw(z, 0) = 0, (D[1](Tw))(1, t) = 0,Tw(0, t) = 0

I shoud solve this equations numerically and plot Th vs.t(and Th vs.z),

please help me.

thanks

Dear All,

I am solving ODE with BC .. then im having this error.. (Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system)

Anyone could help me please?? :)

34.mw

``

restart; with(plots); l := .1; M := .5; B := .6; c1 := 1; R := .5; Pr := .72; N := .5; rho := 1; Bt := .6; Ec := .5; lambda := .1; S := .5; blt := 5

Eq1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta))-2*(diff(f(eta), eta))^2+2*l*B*(diff(F(eta), eta)-(diff(f(eta), eta)))-M*(diff(f(eta), eta))-A*(2*(diff(f(eta), eta))+eta*(diff(f(eta), eta))) = 0;

diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(f(eta), eta))-2*(diff(f(eta), eta))^2+.12*(diff(F(eta), eta))-.62*(diff(f(eta), eta))-A*(2*(diff(f(eta), eta))+eta*(diff(f(eta), eta))) = 0

(1)

Eq2 := F(eta)*(diff(F(eta), eta))-2*(diff(F(eta), eta))^2+2*B*(diff(f(eta), eta)-(diff(F(eta), eta))-A*(eta*(diff(F(eta), eta))-2*(diff(F(eta), eta)))) = 0

F(eta)*(diff(F(eta), eta))-2*(diff(F(eta), eta))^2+1.2*(diff(f(eta), eta))-1.2*(diff(F(eta), eta))-1.2*A*(eta*(diff(F(eta), eta))-2*(diff(F(eta), eta))) = 0

(2)

Eq3 := (1+(4/3)*R)*(diff(theta(eta), eta, eta))+Pr*((diff(theta(eta), eta))*f(eta)-c1*(diff(f(eta), eta))*theta(eta))+2*N*Pr*Bt*(theta1(eta)-theta(eta))/rho+2*N*B*Pr*Ec*(diff(F(eta), eta)-(diff(f(eta), eta)))^2/rho-A*Pr*(4*theta(eta)+eta*(diff(theta(eta), eta)))+2*lambda*Pr*theta(eta)+Pr*Ec*(diff(f(eta), eta))^2 = 0;

1.666666667*(diff(diff(theta(eta), eta), eta))+.72*(diff(theta(eta), eta))*f(eta)-.72*(diff(f(eta), eta))*theta(eta)+.4320*theta1(eta)-.2880*theta(eta)+.21600*(diff(F(eta), eta)-(diff(f(eta), eta)))^2-.72*A*(4*theta(eta)+eta*(diff(theta(eta), eta)))+.360*(diff(f(eta), eta))^2 = 0

(3)

Eq4 := c1*(diff(F(eta), eta))*theta1(eta)-F(eta)*(diff(theta1(eta), eta))+2*gamma*Bt*(theta1(eta)-theta(eta))+A*(4*theta1(eta)+eta*(diff(theta1(eta), eta))) = 0;

(diff(F(eta), eta))*theta1(eta)-F(eta)*(diff(theta1(eta), eta))+1.2*gamma*(theta1(eta)-theta(eta))+A*(4*theta1(eta)+eta*(diff(theta1(eta), eta))) = 0

(4)

bcs1 := f(0) = S, (D(f))(0) = 1, (D(f))(blt) = 0, (D(F))(blt) = 0, F(blt) = f(blt)+blt*(D(f))(blt)-blt*(D(F))(blt), theta(0) = 1, theta(blt) = 0, theta1(blt) = 0;

f(0) = .5, (D(f))(0) = 1, (D(f))(5) = 0, (D(F))(5) = 0, F(5) = f(5)+5*(D(f))(5)-5*(D(F))(5), theta(0) = 1, theta(5) = 0, theta1(5) = 0

(5)

L := [.1, .15, .2];

[.1, .15, .2]

(6)

for k to 3 do R := dsolve(eval({Eq1, Eq2, Eq3, Eq4, bcs1}, A = L[k]), [f(eta), F(eta), theta(eta), theta1(eta)], numeric, output = listprocedure); Y || k := rhs(R[3]); YP || k := rhs(R[4]); YR || k := rhs(R[5]); YQ || k := rhs(R[6]); YA || k := rhs(R[7]) end do

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

 

``

``

 

Download 34.mw

the functions about vector field plot is so terrible,

is there any method to get it beautiful, like mathematica or matlab?

more options is needed for arrow's position, length, angle and color.

is there a plan for it?

And, the scalar field plot in 2d is better than vector field plot, but in 3d, isosurface plot are still need to be enhanced.

If possible, the plots in Plotting Guide are all need improve.

Dear all:

I wisth to replace in an expression that contains several sums:

sum(a[i],i=1..n+1) and sum(w[i],i=1..n+1)

and replace them with a single command into:

a[n+1]+sum([i],i=1..n) and w[n+1]+sum(w[n],1..n)

everywhere in the expression

Thank you

Claudio

PD := [diff(Th(z, t), t) = 7.1428*(diff(Th(z, t), z))-1397941.885*(279-Tw(z, t))-0.2160487e-1*(diff(Th(z, t), z, z)), diff(Tc(z, t), t) = -7.1428*(diff(Tc(z, t), z))+1298990.852*(Tw(z, t)-291)+0.189366e-1*(diff(Tc(z, t), z, z)), diff(Tw(z, t), t) = 3.3024901*(Th(z, t)-2*Tw(z, t)+Tc(z, t))+8.0029*10^(-4)*(diff(Tw(z, t), z, z))]

BCI := {Tc(0, t) = 275, Tc(z, 0) = 275, Th(1, t) = 296, Th(z, 0) = 296, Tw(0, t) = 0, Tw(z, 0) = 0, (D[1](Tc))(1, t) = 0, (D[1](Th))(0, t) = 0, (D[1](Tw))(1, t) = 0}

pds := pdsolve(PD, BCI, numeric)

Hallo. There is a package "Standard Form" (for Maple 3 & 4): http://www.cecm.sfu.ca/~wittkopf/

I am interesting if it can be converted for modern Maple 18?

 

The Maple help definition for spherical coordinates uses the triple (r, φ, θ) (Note the ordering!!) with φ in the range 0..Π and θ in the range 0..2Π. This means that the second entry in the triple is the zenith angle (latitude) and the third entry in th triple is the azimuth angle (longitude). This is confirmed by the relation to cartesian coordinates stated on the definition page as
x= r sin(φ) cos(θ)

y= r sin(φ) sin(θ)

z= r cos(φ)

However the help page for coords has spherical polars defined by the triple (u, v, w), with the relation to cartesian coordinates given as

x= u cos(v) sin(w)

y= u sin(v) sin(w)

z= u cos(w)

which suggests that this time it is the third entry in the triple (ie w) which is the zenith angle (latitude), with the second entry being the azimuth (longitude).

My simple-minded attempt to check which of these interpretations is correct is shown in the attached worksheet. This seems to confirm that the MapleHelp definitions page is correct and the help/coords page is incorrect - or am I missing something??

Download posVecChk.mw

 

 

Aslam-u-Alikum. Hope you will be fine. I have some error in the last command of the attached file. I want this type of answer

 

StressesMultiplication := Matrix(2, 2, {(1, 1) = 7, (1, 2) = 10, (2, 1) = 15, (2, 2) = 22})

n.mw

PhD (Scholar)
Department of Mathematics

f := x^2+y-z=0

f2:= y^2 +z-x = 0

after shift , solution shift too, can it be said it is invariant in parameter shift?

if not, any example to show a function which is invariant in parameter shift?

> solve(f);
/ 2 \
{ x = x, y = -x + z, z = z }
\ /
> f2 := y^2+z-x;
2
y + z - x
> solve(f2);
/ 2 \
{ x = y + z, y = y, z = z }
\ /
> f;
2
x + y - z = 0

First 1336 1337 1338 1339 1340 1341 1342 Last Page 1338 of 2434