MaplePrimes Questions

I search saveplot and sleep function here

i had kept close filename function and move plot setup default to first line

but thread no sleep in maple 12

then I remove sleep and can save image file

but after I call it in a function and run a for loop for this function

the file is not updated after sleep 3 seconds 

how to run a for loop call it and it can refresh 

Is anyone accessing the https://www.mapleprimes.com/questions redirected to this page https://www.mapleprimes.com/errors/500.aspx?aspxerrorpath=/questions/default.aspx ?

https://www.mapleprimes.com/errors/500.aspx?aspxerrorpath=/questions/default.aspx

I tried to control z-axes by writing

plot3d(sech((1/2)*sqrt(3/5)*(x-2*t)), x = -10 .. 10, t = -10 .. 10, scaling = constrained, style = patchnogrid)

plot3d(sech((1/2)*sqrt(3/5)*(x-2*t)), x = -10 .. 10, t = -10 .. 10, scaling = constrained, style = patchnogrid)

 

but I could not change the height of z axes. I want to change the dimensions like in the following figure.

example-graph.pdfexample-graph.pdf

Although the range of z-axix is small but it appears very clear. Can anyone help me, please.

 

The following message is displayed at the beginning of my worksheet when i reopen it, no data is lost because the nature of my code retrives it's optimal boundaries for the loops but i dont know what the hell this is telling me this for:

"Warning, .hdb help databases are deprecated, 'C:\Users\the_r\maple\toolbox\Syrup\lib\Syrup.hdb' will not be used, see ?HelpTools,Migrate help page for more information"

It appears in blue standard size courier font if that is help thank you

Is it asking me to migrate? My passport was stolen by clandestine pharmacists 

 

 Dear fellows
I have the problme with first seq command.
for i while i <= M1-2 do for j while j <= M1-2 do for k while k <= M1 do Eq[i, jk] := simplify(eval(R, [x = i/(M1-2.), y = j/(M1-2.), t = k/(M1-1.)])) = 0 end do end do end do:

Sol := fsolve({seq(BC1[m8, m9]$m8 = 1 .. M1, m9 = 1 .. M1), seq(BC2[m10, m11]$m10 = 1 .. M1, m11 = 1 .. M1), seq(BC3[m12, m13]$m12 = 1 .. M1, m13 = 1 .. M1), seq(BC4[m14, m15]$m14 = 1 .. M1, m15 = 1 .. M1), seq(IC1[m4, m5]$m4 = 0 .. M1, m5 = 1 .. M1), seq(IC2[m6, m7]$m6 = 0 .. M1, m7 = 1 .. M1), seq(Eq[m1, m2, m3]$m1 = 1 .. M1-2, m2 = 1 .. M1-2, m3 = 1 .. M1-1)});

Got the following error.

Error, invalid input: seq expects its 3rd argument, step, to be of type numeric, but received m3 = 1 .. 1


Please help me in this regard.

Error, (in Optimization:-NLPSolve) the objective gradients at the initial point are too small

So this error indicates that the initial point is a local minimum or very very close to it right?

I'm very new to Maple (and coding in general). I'd like to know a very large coefficient of a certain 𝑞-series involving the Dedekind 𝜂-function. For example

s:=series(subs(q=q^6, eta^4), q, N);

where eta := series(q^(1/24)*product(1 - q^n, n = 1 .. 100), q, N). I know the command coeftayl(s, q=0, N) but for certain 𝑞-series it doesn't recognize the numerical value of this coefficient when 𝑁 is too large (even after subtracting the principal part of 𝑠, if there is any) - instead it realizes the coefficient as a HUGE limit of something. Cheers.

I need the final line of the uploaded worksheet to be the "Select Action", preprogrammed prior to the execution of the insert xml command, how do i do this?

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/ASK.mw .
 

Download ASK.mw

I am trying to do matrix multiplication with Dirac gamma matrices, but am getting a recursion error. 

the final input line, which is causing the error, doesn't seem to be inserted, but it is:

Library:-PerformMatrixOperations((mt . (Dgamma[2])) . m)

can anybody tell me what i'm doing wrong? i don't understand why 

Dgamma[2] . m 

works but 

mw . Dgamma[2] . m 

gets the recursion error. 

then i saw something that made me think i needed to use 'matrix', instead of 'Matrix'. so i tried that(gamult2) - now the 2nd multiply works, and the first fails.

any help greatly appreciated!

thanks,

larry

restart

with(Physics)

Physics:-Version()

`The "Physics Updates" version in the MapleCloud is 452 and is the same as the version installed in this computer, created 2019, October 23, 13:26 hours, found in the directory /Users/lstead/maple/toolbox/2019/Physics Updates/lib/`

(1)

with(LinearAlgebra)

``

m := Matrix(4, 1, {(1, 1) = 1, (2, 1) = 2, (3, 1) = 3, (4, 1) = 4})

Matrix(%id = 18446744078528784254)

(2)

mt := Transpose(m)

Matrix(%id = 18446744078392267830)

(3)

mt.m

Matrix(%id = 18446744078392260238)

(4)

Library:-PerformMatrixOperations(Dgamma[2].m)

Matrix(%id = 18446744078528770526)

(5)

Library:-PerformMatrixOperations(mt.Dgamma[2].m)

Error, (in Physics:-GetDefineINFO) too many levels of recursion

 

``


 

Download gammult.mw

gammult2.mw

hello

why this pdsolve not be answer

pde := diff(T(x, y, t), x, x)+diff(T(x, y, t), y, y) = (1/0.44e-2)*(diff(T(x, y, t), t));
conds := T(0, y, t) = 1, T(1, y, t) = 1, T(x, 0, t) = 1, T(x, 1, t) = 1, T(x, y, 0) = 0;
`assuming`([pdsolve([pde, conds])], [0 <= x and x <= 1, 0 <= y and y <= 1, t >= 0]);
 

Dear all

I have a fourth-order PDEs, with some boundary condition, I would like to compute an explicit solution or plot the solution in the interval [0,2]

 

PDEsolve.mw

thinks

Hello everyone,

Say I have the differential equation:

x1'(t) = 3C1*cos(t) + 2C2*sin(t)

x2'(t) = 5C1*cos(t) + 1C2*sin(t)

How can one draw phase portraits in maple?

 

 

I have a piecevise defined function where i wanna reference the condition inorder to use the maximum and minimum values of the defined range in other functions. Is there a function for this that i just cant find or do i have to use a workaround?

These are the values im trying to access. and heres the function:

f := proc (t) options operator, arrow; piecewise(-1 < t and t <= 0, t+1, 0 < t and t <= 1, 1-t) end proc

Hi,

I have written a module M which contains some procedures P1, P2, ...
This module contains the statement option package.

In a separate worksheet I type the command L := with(M) to load this module and stack in the list L the names of all the procedures M contains.
I would like now to use showstat to display the content of all these procedures..

I thus typed:
for p in L do  showstat(p) end do;

I got an error saying that showstat expects a name as argument.
In the command above p is not a name but a symbol.

Is there a possibility to "showstat" all the procedures M contains ?

Motivation:I'd like to count the number of command lines in module M (showstat removes the lines of comments) ... maybe there exist an other way than using showstat ?


Thanks in advance

from help, it says about option of 'spacestep'=numeric to pdsolve numerical solver the following

Specifies the spacing of the spatial points on the discrete mesh on which the solution 
is computed, and defaults to 1/20th of the spatial range of the problem, 
r-L, where L is the left boundary, and r is the right. Note: The spacing must 
be small enough that a sufficient number of points are in the spatial domain 
for the given method, boundary conditions, and spatial interpolation (see below). 
If the given value of spacestep does not fit into the input domain an integral 
number of times, the closest smaller step size that does is used. For problems 
using error estimates, or an adaptive approach, the total number of spatial 
points must be odd, so again, if this is not the case for the specified 
spacestep, then it is reduced to satisfy this requirement.

In this problem, the domain is [-1,1]. So the default is 1/20 of this which is 1/10.  Why is it when using spacestep smaller than 1/11, the animator generated from the solution does not work right?  i.e. playing the animate does not produce correct result as the case when using spacestep=1/11 or spacestep=1/10?   

When using spacestep=1/16 or 1/15 or 1/14 or 1/13 or 1/12, it all produce bad animation.  I do not see why that is. 

Any ideas?  Am I doing something wrong?
 

pde := diff(u(x,t),t$2)=diff(u(x,t),x$2):
bc  := u(-1,t)=u(1,t),D[1](u)(-1,t)=D[1](u)(1,t):
f:=x->piecewise(-1/2<x and x<0,x+1/2,0<x and x<1/2,1/2-x,true,0):
plot(f(x),x=-1..1):
ic  := u(x,0)=f(x),D[2](u)(x,0)=0:
sol:=pdsolve(pde,{bc,ic},u(x,t),numeric,time=t,range=-1..1,  timestep=1/16,spacestep=1/16)

_m1793827573184

sol:-animate(t=1,frames=50,title="time = %f");

sol:=pdsolve(pde,{bc,ic},u(x,t),numeric,time=t,range=-1..1,  timestep=1/16,spacestep=1/12):
sol:-animate(t=1,frames=50,title="time = %f");

sol:=pdsolve(pde,{bc,ic},u(x,t),numeric,time=t,range=-1..1,  timestep=1/16,spacestep=1/11):
sol:-animate(t=1,frames=50,title="time = %f");

 

 

Download solving_wave_pde_2.mw

 

 

First 510 511 512 513 514 515 516 Last Page 512 of 2308