MaplePrimes Questions

Sometimes we have questions that apply to every execution we make, regardless of what the actual nature of the commands included in that execution are.

 

For this reason, I propose that in the same way that we have the feature of "start up" code, which automatically executes upon opening the specified worksheet, we also implement an option of automatically executed code for every execution that is carried out within the specified worksheet. 

for example, code that serves purposes of output such as 

" This execution took xxxx minutes and xx.x seconds to complete, and required an average of x % of your computer's total processor utilization capacity across all cores "

I know how i will be providing this feature in my worksheets, but i just feel that it would be a beneficial feature for those that use maple for purposes that do not involve a knowledge of how maple is programmed, and would not be able to do so.

Any discussion about this would be welcomed, which is why i am posting about it of course.

I am trying to see if Maple can solve Laplace PDE inside the disk in polar coordinates. Standard textbook problem. Radius of disk is `a`. The boundary conditions on the disk is `f(theta)`. One of the conditions needed also is that the solution is finite in the center of the disk.

I do not know how to tell Maple that the solution should be finite in the center of the disk. If I do not give this conditions, Maple gives me strange looking solution, which does not look like anything close to the standard series solution one gets from hand solution. There is not even a series solution.

This is what I tried

restart;
pde:=diff(u(r,theta),r$2)+1/r*diff(u(r,theta),r)+1/r^2*diff(u(r,theta),theta$2)=0;
bc:=u(a,theta)=cos(theta);
sol:=pdsolve([pde,bc],u(r,theta)) assuming r<=a,r>0

Now, how to tell it that `u(0,theta)` is bounded? So that the `ln(r)` solution do not show up? Adding `u(0,theta)<infinity` to the boundary conditions, gives error

restart;
pde:=diff(u(r,theta),r$2)+1/r*diff(u(r,theta),r)+1/r^2*diff(u(r,theta),theta$2)=0;
bc:=u(a,theta)=cos(theta),u(0,theta)<infinity;
sol:=pdsolve([pde,bc],u(r,theta)) assuming r<=a,r>0

The standard solution to this PDE is

Where `c0` and `cn` and `kn` above are found from boundary conditions at $u(a,\theta)$.

How can one get Maple to give the above solution? How to tell it that $u$ is bounded at $r=0$?

 

 

we know that the exponotial function is represented as a taylor series, that means that the exponontial of any value is approximatvly evaluated ,the same thing of sinus function when i calcute sin(Pi/6) it give me the exacte result,
 

sum(0^i/i!,i=0..infinity);

0

(1)

exp(0);

1

(2)

sin(0.2);

.1986693308

(3)

sin(Pi/6);

1/2

(4)

sin(Pi/2);

1

(5)

 


 

Download trying.mws

,can anyone explain to me these results ?

Came across this issue again.  While working within a worksheet, everything seems fine, using commands like imagetools outputing bmp files etc.  I had to, in the the Options->Precision, uncheck limit expression length to 1000000 during my session.  I saved it multiple times, however and unfortunately did not save as seperate versions, so I was left with one file.

The file size ended up being 64Mb.  The problem is trying to load it into Maple, everytime I try to load the worksheet my computer seems to freeze (laptop, 4Gb RAM, windows 7 64 bit, Maple 2017) I check task manager and my javaw.exe file is consuming 3.7Gb of memory!  So that's why it's locking up. 

I've tried multiple times to open the file and utilmately end up pressing and holding the power button to restart my system. 

 

C: y^2=10161/256*t^2+8829189/25600*t+7266953961/1024000 

how to find the (t, y) points on the curve C ?

THANKS.

Attractor.mw
Hello, anybody could help with this error?

restart

interface(imaginaryunit = j);

I

(1)

with(DEtools):

k := .5; 1; tau := .95; 1; mu := 0.1e-1; 1; pi := 116.1; 1; vartheta := 0.8e-2; 1; phi := 0.25e-2; 1; epsilon := 0.2e-2; 1; rho := 0.5e-1; 1; beta := 0.115e-1; 1; chi := 0.598e-2; 1; q := .5; 1; eta := .2; 1; delta := .1; 1; alpha := 0.57e-1; 1; p := .2; 1; Upsilon := 1.2

.5

 

.95

 

0.1e-1

 

116.1

 

0.8e-2

 

0.25e-2

 

0.2e-2

 

0.5e-1

 

0.115e-1

 

0.598e-2

 

.5

 

.2

 

.1

 

0.57e-1

 

.2

 

1.2

(2)

lambda := k*tau*(C(t)*Upsilon+I(t))/(S(t)+V(t)+C(t)+I(t)+R(t)):

DEplot3d({d*C(t)/dt = rho*lambda*S(t)+rho*`&epsilon;`*lambda*V(t)+(1-q)*eta*I(t)-(mu+beta+chi)*C(t), d*I(t)/dt = (1-rho)*lambda*S(t)+(1-rho)*`&epsilon;`*lambda*V(t)+chi*C(t)-(mu+alpha+eta)*I(t), d*R(t)/dt = beta*C(t)+q*eta*I(t)-(mu+delta)*R(t), d*S(t)/dt = (1-p)*pi+phi*V(t)+delta*R(t)-(mu+lambda+`&vartheta;`)*S(t), d*V(t)/dt = p*pi+`&vartheta;`*S(t)-(lambda*`&epsilon;`+mu+phi)*V(t)}, {C(t), I(t), R(t), S(t), V(t)}, t = 0 .. 300, stepsize = .1, [[S(0) = 8200, V(0) = 2800, C(0) = 200, I(0) = 210, R(0) = 200]], linecolour = t, axes = BOXED, scene = [S, V, C, I, R])

Error, (in DEtools/DEplot/CheckInitial) the 'number' option must be specified before initial conditions

 

``


 

Download Attractor.mw

 

How to convert maple file into the pdf file?For example the attacahed file is maple file.I want to convert this to a pdf file.

conservation_of_wave_eq.mw

If i was a good student i wouldnt be using the is function for these at all, but im quite tired now so i thought that it would be fine this time. But i honest cannot see what additional declarations need to be made here, and why a false is returned for the first query, then FAIL for the negation of the first in the second query, the only thing i can think of is declaring that indeed n and m are elements of N, but if maple was allowing for Gaussian n,m surely both queries would return fail, and so obvious it is indeed acknowledging the arguments are on R, anyway hopefully someone gets the chance to make me feel silly here

 


 

assume(n < m, n > 1)

is((3*n-1)/(n^3+2*n+1)-(3*m-1)/(m^3+2*m+1) < 0)

false

(1)

is((3*n-1)/(n^3+2*n+1)-(3*m-1)/(m^3+2*m+1) >= 0)

FAIL

(2)

``


 

Download MISSING_ASSUMPTION.mw

I'm stuck with applying boundary conditions for a cantilever beam please help with format

What is the best method of transferring figures of maple 18 to latex ?

Hi all,

I don't know why this file doesn't work anymore! It was ok last week. I even reinstall the Maple on my laptop, but it went wrong again...

 

Please help me.

 

Plot_Diff_Temperature.mw

Can Maple solve the chief or nurse shift problems? 

I have a positive function and I'm taking an integral from it. I don't get why Maple computes the integral over a bigger region but not over a subset of it! Here is the integral. Does anyone know what is going wrong?

int(15.91549431*abs(6.100147120*10^6/k[6]-2.027358424*10^6*(1.241833982*k[6]/T[1]+2.)/k[6]+2.002942410*10^6/T[1])*exp(-50.00000000*(1.241833982*k[6]/T[1]-1.5)^2)/(k[6]^2*T[1]^2), T[1] = -.5 .. 0, k[6] = -.5 .. 0, numeric)
int(15.91549431*abs(6.100147120*10^6/k[6]-2.027358424*10^6*(1.241833982*k[6]/T[1]+2.)/k[6]+2.002942410*10^6/T[1])*exp(-50.00000000*(1.241833982*k[6]/T[1]-1.5)^2)/(k[6]^2*T[1]^2), T[1] = -.9 .. 0, k[6] = -.9 .. 0, numeric)

Hi all

I have a special 2*2 block martix like follwing:

 

 

how can i define this by maple?

thanks in advance.

I'm optimizing some code and I have noticed that running a procedure within the Threads:-Task module seems to slow the procedure down. I want to use Threads:-Task to run many procedures in parallel and so finish a task faster.
I'm attaching an example of a simple "add" procedure run outside the Task module and within the Task module. As far as I know, the procedure is thread-safe, and this is shown by the similar time it takes to run 1 or 10 simultaneously.
What is going wrong?

``

restart

addproc := proc () local i, A; A := 0; A := add(i, i = 1 .. 1000); return A end proc

ts := time[real](); for t to 1000 do addproc() end do; TA := time[real]()-ts

ts := time[real](); for t to 1000 do Threads:-Task:-Start(null, Task = addproc) end do; TS := time[real]()-ts

ts := time[real](); for t to 1000 do Threads:-Task:-Start(null, Tasks = [addproc, addproc, addproc, addproc, addproc, addproc, addproc, addproc, addproc, addproc]) end do; TM := time[real]()-ts

TA; TS; TM

0.33e-1

 

2.049

 

2.006

(1)

TM/TA

60.78787879

(2)

``


 

Download TaskExample.mw

First 853 854 855 856 857 858 859 Last Page 855 of 2427