Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi everyone

I'd be pleased if you could give a hand with the exploration assistant.

 

1. I want the exploration assistant to appear on the same document I am working on, but everytime I use it (either by right-clicking or by the explore command) it automatically appears on a new document.

 

2. can I manipulate a piecewise function when using embedded components? i.e.: plot the function and varying the parameters using sliders.

 

thanks

cesar

I have made a maplet. The aim is to take two numbers in textboxes and find their Jacobi symbol. The result should appear in a third textbox. Now I have done all the steps and called a procedure on buttonclick. Following is the procedure I have made

jacob := proc (x, y)

local A, B;

A := x;

B := y;

with(numtheory);

jacobi(A, B)

end proc:

 

Problem: The problem is that when I make any small modification to the code it runs and shows the result in the 3rd textbox but as soon as I save it, restart maple and recompile the program, the textbox 3 will show jacobi(value of textbox1, value of textbox2) instead of the result. What should I do, actually I am new to maple. Please guide me.

k := a*x6^7 + b*x6^6 + c*x6^5 + d*x6^4 + e*x6^3 + f*x6^2 + g*x6;
discrim(k, x6);

after know the discriminant equation, how to do next?

solve(Max(Min(x, 1), Min(x, 2)));
solve(Min(Max(x, 1), Max(x, 2)));

i do not know how to use Rootof something

 

any other expression

Please, I solved a pde system of equation problem numerically, using maple 17.

But I dont know how to plot multiple solutions on one graph.

I want to vary one of the parameters....

e.g Pr=0.71, Pr=7, Pr=10 where other parameters are kept constant

 

My working is attachedtobi_msc_solution.mw

restart

M := 1:

pde1 := diff(u(y, t), t)+Typesetting:-delayDotProduct(S, diff(u(y, t), y))-2*k^2*u(y, t) = diff(u(y, t), y, y)+theta(y, t)+Typesetting:-delayDotProduct(N, C(y, t))+Typesetting:-delayDotProduct(M, u(y, t))+u(y, t)/K:

                pde2 := theta(y, t)+t*(diff(theta(y, t), t))+S*(diff(theta(y, t), y)) = (diff(theta(y, t), y, y))/Pr-Typesetting:-delayDotProduct(alpha, theta(y, t)):

pde3 := C(y, t)+t*(diff(C(y, t), t))+S*(diff(C(y, t), y)) = (diff(C(y, t), y, y))/Sh-Typesetting:-delayDotProduct(R, C(y, t)):

PDE := {pde1, pde2, pde3}:

IBC := {C(0, t) = 1, C(1, t) = 0, C(y, 0) = 0, u(0, t) = 0, u(1, t) = 0, u(y, 0) = 0, theta(0, t) = 1, theta(1, t) = 0, theta(y, 0) = 0}:

pds := pdsolve(PDE, IBC, numeric)

module () local INFO; export plot, plot3d, animate, value, settings; option `Copyright (c) 2001 by Waterloo Maple Inc. All rights reserved.`; end module

(1)

pds:-plot[display](u(y, t), t = .5, linestyle = "solid", colour = "blue", legend = "Pr=0.71", title = "Velocity Profile", labels = ["y", "theta"])

 

``


Download tobi_msc_solution.mw

 

Please, Any help will be gracefully appreciated

 

hey i am new here and i have a question,

i have an partial differential equation diff(u(t,x),t$2)=diff(u(t,x),x$2) with the intial value problem u(0,x)=f(x)=1/(1+x^2) and diff(u(0,x),t)=0

and now my question i have already programm it:


> with(inttrans);

with(DEtools);

with(plots);


> with(PDETools);

> k := diff(u(t, x), `$`(t, 2)) = diff(u(t, x), `$`(x, 2));
bc := u(0, x) = 1/(1+x^2);
v := diff(u(0, x), t) = 0;
d / d \ d / d \
--- |--- u(t, x)| = --- |--- u(t, x)|
dt \ dt / dx \ dx /
1
u(0, x) = ------
2
1 + x
0 = 0

> pdsolve(k, u(t, x));
print(`output redirected...`);
u(t, x) = _F1(x + t) + _F2(x - t)

> c := pdsolve({bc, k, v}, u(t, x));
print(`output redirected...`); # input placeholder

 and now question at the last there is nothing does it means that maple can´t solve it with the intial value problem and how can solve it with Fourier-Transformation to x???

can anyone help me please and sorry my englisch is not so good ;)

Hello, 

I have a trigonometric equation.

I would like to isolate gamma[1](t) and to determine gamma[1](t) in fonction of alpha(t), beta(t) and z(t). The others variables in the equations are fixed parameters.

I have tried to use isolate function. But it doesn't work.

Of course, my expressions should be complex but that is not a problem if i manage to expresse gamma[1](t) in fonction of alpha(t), beta(t) and z(t).

Here my program

constraints_2.mw

Thank you for you help

superposition said that a vector is a linear combination of other vectors

but even if i calculated the coefficient, i do not know which vector is which other vectors's linear combination

how to prove?

InputMatrix3 := Matrix([[close3(t), close3(t+1) , close3(t+2) , close3(t+3) , close3(t+4) , close3(t+5)],
[close3(t+1) , close3(t+2) , close3(t+3) , close3(t+4) , close3(t+5) , close3(t+6)],
[close3(t+2) , close3(t+3) , close3(t+4) , close3(t+5) , close3(t+6) , 0],
[close3(t+3) , close3(t+4) , close3(t+5) , close3(t+6) , 0 , 0],
[close3(t+4) , close3(t+5) , close3(t+6) , 0 , 0 , 0],
[close3(t+5) , close3(t+6) , 0 , 0 , 0, 0],
[close3(t+6) , 0 , 0 , 0, 0, 0]]):
EigenValue1 := Eigenvalues(MatrixMatrixMultiply(Transpose(InputMatrix3), InputMatrix3)):
Asso_eigenvector := Eigenvectors(MatrixMatrixMultiply(Transpose(InputMatrix3), InputMatrix3)):
AEigenVector[tt+1] := Asso_eigenvector;

Matrix(6, 6, {(1, 1) = .514973850028629+0.*I, (1, 2) = .510603608194333+0.*I, (1, 3) = .469094659512372+0.*I, (1, 4) = .389872713818831+0.*I, (1, 5) = .279479324327359+0.*I, (1, 6) = -.154682461176604+0.*I, (2, 1) = .493994413154560+0.*I, (2, 2) = .306651336822139+0.*I, (2, 3) = -0.583656699197969e-1+0.*I, (2, 4) = -.417550308930506+0.*I, (2, 5) = -.566122865008542+0.*I, (2, 6) = .404579494288380+0.*I, (3, 1) = .449581541124671+0.*I, (3, 2) = -0.266751368453398e-1+0.*I, (3, 3) = -.529663398913996+0.*I, (3, 4) = -.359719616523673+0.*I, (3, 5) = .313717798014566+0.*I, (3, 6) = -.537405340038665+0.*I, (4, 1) = .386952162293470+0.*I, (4, 2) = -.351332186748244+0.*I, (4, 3) = -.390816901794187+0.*I, (4, 4) = .470032416161955+0.*I, (4, 5) = .231969182174424+0.*I, (4, 6) = .547134073332474+0.*I, (5, 1) = .306149178348317+0.*I, (5, 2) = -.530611390076568+0.*I, (5, 3) = .192717713961280+0.*I, (5, 4) = .291213691618787+0.*I, (5, 5) = -.562991429686901+0.*I, (5, 6) = -.431067688369314+0.*I, (6, 1) = .212576094920847+0.*I, (6, 2) = -.489443150196337+0.*I, (6, 3) = .553283259136031+0.*I, (6, 4) = -.488381938231088+0.*I, (6, 5) = .363604594054259+0.*I, (6, 6) = .195982711855368+0.*I})

Matrix(6, 6, {(1, 1) = .515428842592397+0.*I, (1, 2) = .515531996615269+0.*I, (1, 3) = .468108280940919+0.*I, (1, 4) = -.392394120975052+0.*I, (1, 5) = -.280467124908196+0.*I, (1, 6) = -.129613084502380+0.*I, (2, 1) = .494563493180197+0.*I, (2, 2) = .301273494494509+0.*I, (2, 3) = -0.622136916501293e-1+0.*I, (2, 4) = .438383262732459+0.*I, (2, 5) = .571041594120088+0.*I, (2, 6) = .377494770878435+0.*I, (3, 1) = .450886315308369+0.*I, (3, 2) = -0.323387895921418e-1+0.*I, (3, 3) = -.527636820417566+0.*I, (3, 4) = .332744872607714+0.*I, (3, 5) = -.322934536375586+0.*I, (3, 6) = -.549772001891837+0.*I, (4, 1) = .385916641681991+0.*I, (4, 2) = -.352066020655722+0.*I, (4, 3) = -.389655495441319+0.*I, (4, 4) = -.450049711766943+0.*I, (4, 5) = -.221529986447276+0.*I, (4, 6) = .568916672007495+0.*I, (5, 1) = .305485655770791+0.*I, (5, 2) = -.528766119966973+0.*I, (5, 3) = .201065789602278+0.*I, (5, 4) = -.310329356773806+0.*I, (5, 5) = .555973984740943+0.*I, (5, 6) = -.425730045170186+0.*I, (6, 1) = .210210489500614+0.*I, (6, 2) = -.488744465076970+0.*I, (6, 3) = .553484076328700+0.*I, (6, 4) = .494245653290329+0.*I, (6, 5) = -.364390406353340+0.*I, (6, 6) = .183130120876843+0.*I})
mm1 := 1;
solve(
[AEigenVector[mm1][2][1][6] = m1*AEigenVector[mm1][2][1][1]+m2*AEigenVector[mm1][2][1][2]+m3*AEigenVector[mm1][2][1][3]+m4*AEigenVector[mm1][2][1][4]+m5*AEigenVector[mm1][2][1][5],
AEigenVector[mm1][2][2][6] = m1*AEigenVector[mm1][2][2][1]+m2*AEigenVector[mm1][2][2][2]+m3*AEigenVector[mm1][2][2][3]+m4*AEigenVector[mm1][2][2][4]+m5*AEigenVector[mm1][2][2][5],
AEigenVector[mm1][2][3][6] = m1*AEigenVector[mm1][2][3][1]+m2*AEigenVector[mm1][2][3][2]+m3*AEigenVector[mm1][2][3][3]+m4*AEigenVector[mm1][2][3][4]+m5*AEigenVector[mm1][2][3][5],
AEigenVector[mm1][2][4][6] = m1*AEigenVector[mm1][2][4][1]+m2*AEigenVector[mm1][2][4][2]+m3*AEigenVector[mm1][2][4][3]+m4*AEigenVector[mm1][2][4][4]+m5*AEigenVector[mm1][2][4][5],
m1^2 + m2^2 + m3^2 + m4^2 + m5^2 = 1], [m1, m2, m3, m4, m5]);

[m1 = .4027576723+.5022235499*I, m2 = -.5922841426-1.043213223*I, m3 = -.1130969773+.9150300317*I, m4 = .9867039883-.5082455178*I, m5 = -1.400123192+.1536850673*I], [m1 = .4027576723-.5022235499*I, m2 = -.5922841426+1.043213223*I, m3 = -.1130969773-.9150300317*I, m4 = .9867039883+.5082455178*I, m5 = -1.400123192-.1536850673*I]

mm1 := 2;
solve(
[AEigenVector[mm1][2][1][6] = m1*AEigenVector[mm1][2][1][1]+m2*AEigenVector[mm1][2][1][2]+m3*AEigenVector[mm1][2][1][3]+m4*AEigenVector[mm1][2][1][4]+m5*AEigenVector[mm1][2][1][5],
AEigenVector[mm1][2][2][6] = m1*AEigenVector[mm1][2][2][1]+m2*AEigenVector[mm1][2][2][2]+m3*AEigenVector[mm1][2][2][3]+m4*AEigenVector[mm1][2][2][4]+m5*AEigenVector[mm1][2][2][5],
AEigenVector[mm1][2][3][6] = m1*AEigenVector[mm1][2][3][1]+m2*AEigenVector[mm1][2][3][2]+m3*AEigenVector[mm1][2][3][3]+m4*AEigenVector[mm1][2][3][4]+m5*AEigenVector[mm1][2][3][5],
AEigenVector[mm1][2][4][6] = m1*AEigenVector[mm1][2][4][1]+m2*AEigenVector[mm1][2][4][2]+m3*AEigenVector[mm1][2][4][3]+m4*AEigenVector[mm1][2][4][4]+m5*AEigenVector[mm1][2][4][5],
m1^2 + m2^2 + m3^2 + m4^2 + m5^2 = 1], [m1, m2, m3, m4, m5]);

[m1 = .4262845394-.5114193433*I, m2 = -.6313720018+1.072185334*I, m3 = -0.7337582213e-1-.9580760394*I, m4 = -1.036525681-.5400714113*I, m5 = 1.412710014+.1874839516*I], [m1 = .4262845394+.5114193433*I, m2 = -.6313720018-1.072185334*I, m3 = -0.7337582213e-1+.9580760394*I, m4 = -1.036525681+.5400714113*I, m5 = 1.412710014-.1874839516*I]

If A is a matrix 2*2 then how can decompose A as spectral decomposition.

I want to write maple code of the following algorithm with

the following parameters and initial values please help me.

T0 = 5.5556 × 107 cells, I0 = 1.1111 × 107 cells, V0 = 6.3096 × 109 copies/ml,

A1=A2=1,

c = 0.67, h = 1, d = 3.7877 × 10−3, δ = 3.259d,

λ = 2/3× 108d, R0 = 1.33,

p = (cV0δR0)/λ(R0−1)

and β = dδcR0/λp .

 

Algorithm
step 1 :
T(0) = T0, I(0) = I0, V (0) = V0 λi(100 ) = 0 (i=1, ..., 3), u1(0) = 0 =
u2(0).

step 2 :
for i=1, ..., n-1, do :
Ti+1=(Ti + hλ)/(1 + h[d + (1 − u1i)βVi]),

Ii+1 =(Ii + h(1 − u1i)βViTi+1)/(1 + hδ),

Vi+1 =(Vi + h(1 − u2i)pIi+1)/(1 + hc),

λ1n−i−1 =(λ1n−i + h[1 + (1 − u1i)βVi+1])/(1 + h[d + (1 − u1i)βVi+1]),

λ2n−i−1 =(λ2n−i+ hλ3n−i (1 − u2i)p)/(1 + hδ),

λ3n−i−1 =(λ3n−i + h(λ2n−i−1− λ1n−i−1 )(1 − u1i)βTi+1)/(1 + hc),

R1i+1 =(1/A1)(λ1n−i−1−λ2n−i−1 )βVi+1Ti+1,

R2i+1 =−(1/A2)λ3n−i−1 pIi+1,

 

u1i+1 = min(1, max(R1i+1 , 0)),

u2i+1 = min(1, max(R2i+1 , 0)),

end for

 

step 3 :
for i=1, ..., n-1, write
T∗(ti) = Ti, I∗(ti) = Ii, V ∗(ti) = Vi,

u1∗(ti) = u1i, u2∗(ti) = u2i.

end for

hello. before I used Mapple 15. But then I`ve run Mapple 16 and now I`ve a problem. I can`t use this program. I open the program, everthing is in the rule, but if I want to write any mathemathical function, or a letter, such as- x or x+2, the program does`t give any reaction. program only gives reaction the numbers.

Please, help me. (my english isn`t very good, and I don`t know I`ve explained my opinion).

Greetings!

I only recently started to work with Maple17 and tried to test the Explore command. In Classic worksheet I tried a very simple function with a parameter, but after setting the inteval of the parameter in the Java pop-up window I have no result. In return I take: "Error, (in Explore) invalid input: rtable_dims uses a 1st argument, A (of type rtable), which is missing"

Has anyone encounter the same issue before?

Thanx

https://drive.google.com/file/d/0B2D69u2pweEvMV92SGhtRGZONFk/edit?usp=sharing

a error and code in this attachment mw

i can pdsolve it, but numeric pdsolve it get error

hey maple followers,

i need, please to find the method used in command "minimize".
i looked into help maple and i found "theorema mean values" as example
some help please
thanks!

P.s: minimize not Minimize

First 1420 1421 1422 1423 1424 1425 1426 Last Page 1422 of 2261