eq1:=( d)/(dt)u+(d^(2))/(dy^(2))u + s*( d)/(dy)u + delta * theta = 0;

eq2:=( d)/(dt)theta + (d^(2))/(dy^(2))theta + s*Pr*( d)/(dy)theta +lambda* exp(theta/(1 +(epsilon*theta))) = 0; 

initial and boundary conditons   

t <=0; u = theta = 0, for 0 <= y  <= 1   

t> 0;  u =0, theta = 0   at  y = 0;  

t> 0;  u =1, theta = 0   at   y = 1  ;

where, s, epsilon, Pr, lambda, delta are arbitrary parameters

Assuming I have an expression

you can know that the derivative with respect to y will vanish when y is a constant y0.

How can we use the mathematical software Maple to simplify this expression?

There appears to be a bug in Maple 2023 where constants beginning with a Greek letter cannot be saved.  Assigning a constant to a Greek letter works in the Maple interface but does not work when reading an MPL file encoded as UTF-8.  Example:

read "D:/Nikki/docs/Maple/constants.mpl" :

MPL file contents:

alias(asec=arcsec) : # arcsecant

ζX1 := asec(1) : # first airmass zenith angle
ζX2 := asec(2) : # second airmass zenith angle

The constant "ζX1" can be assigned/saved in Maple but shows as unassigned when imported from an MPL file encoded as UTF-8.

Any assistance in how to use Greek letters in MPL files is much appreciated.

Hello everybody, 

So, i managed to use the examples of https://www.mapleprimes.com/questions/233343-Drawing-Of-Complex-Numbers-On-A-Complex-grid- to get to a figure that works for me. If there is a better way, i would love to hear about it. 

I believe the question does not need translation because it is very clear due to the other examples that have been done before example 3. 

The real question is: how to combine these parameters to form the right figure? The parameters are in the red box in the 3rd picture. 


 

Download Toegepaste_Wiskunde_Deel_2_Hoofdstuk_6.mw

Hi all,

I would like to use the command ConvertIn() from the Galois Field package to convert polynomials to field elements.

This is what I've done:

G:=GF(2,4)
                        G := &Fopf;[16]

b:=G:-random()
                          b := T mod 2

s:=G:-ConvertIn(T)
Error, (in ConvertIn) only integer polynomials in T can be converted

Why is this error showing up?? My argument of the command ConvertIn IS an integer polynomial in T.

How can I avoid this error?

Thanks for your help!

BR David

Hello,

I am generating a code to describe how drug dosages (starting from 1 until dosage 2) are metabolized by the body.

The plot looks weird becasue it does not start from C0 = 4. I don't know how to set t=L so that C0 starts at 4. I suspect that I am misssing something in the piecewise function?  

My code is the follwoing:

C0 := 4;
k := 0.2;
L := 1;
C := t -> C0*(1 + exp(-k*L))*exp(-k*(t - L));
L := piecewise(L <= t and t < 2*L, 0, C(t));

plot(C(t), t = 1 .. 10) 

Thanks,

Andre

I'm trying to do a scatter diagram together with a line graph in one. But the function I putting in is not taking it. Is there anyone that can help me with this problem. And thank you.

I know trig functions are hard to simplify in general but was wondering if there was any simplification command that comes in handy for trig functions of the following form.

sin(theta)^(A - 2)*cos(theta)^2 - sin(theta)^(A - 2) + sin(theta)^A

I am trying to make a very rough animation of a rotating pulsar magnetosphere that would very roughly resemble this gif.
I was able to make a sphere and the dipolar field lines. I was wondering if anyone could provide any possible suggestion regarding.

1. How do i animate the rotation such that the rotation axis is not z but slightly misaligned to z.

2. Is there a better way to make the field likes from the poles (jets).

Download Pulsar_temp.mw

3. Any other suggestions that could help with the aesthetics.

Maple 2023 still can't verify its own series solution. This happens with Frobenius method when the indicial root are repated. I think I mentioned or reported this long time ago but unable to find a link now.

My question is: Since this looks like it will not be fixed by Maplesoft anytime soon, any one could suggest an alternative method to verify this solution?

restart;

ode:=2*x^2*diff(y(x), x, x) + 2*x*diff(y(x), x) - x*y(x) =0;
sol:=dsolve(ode,y(x),'series');
odetest(sol,ode,'series')

Maple's solution is correct, I solved it by hand and get same solution.

Maple 2023 on windows 10.

I could not find a way to check whether or not a given value belongs to an interval expressed as RealRange (or a union of them e.g. RealRange(Open(-10),-5), RealRange(Open(-6),infinity).

Ideally, I would like a Boolean function that takes two arguments (i) a value v and (ii) a (union of) RealRange r, and returns true/false depending on whether or not value v is in/out of the range r.

hallo  every body 
please i have a simple question 
how do i calculate this sum in maple 18 
thank you 

Download no_isolation.mw

I have access to a powerful GPU, a NVIDIA Tesla A100.

How to leverage it to solve my system of 3 nonlinear equations in 3 variables? Script: 160523_stylized_problem_GPU.mw

I don't think that it is just by adding the following two lines right before my solve() block, right?
CUDA:-Enable(true)
CUDA:-IsEnabled()
(https://www.maplesoft.com/products/maple/features/cuda.aspx)

I also read https://mcsr.olemiss.edu/docs/gpus-on-maple/ and https://www.mapleprimes.com/posts/36412-GPU-Programming-CUDA-OpenCL-And-Maple but it's still not clear to me. Thanks!

When trying to load/import/read an MPL file with Maple 2023, I get an error message of "Error, recursive assignment" which only links to a page which states "There is no help page available for this error."

https://www.maplesoft.com/support/help/errors/view.aspx?path=Error,%20recursive%20assignment

Presumably there is an error in the Maple syntax that prevents execution, but I cannot find a way to debug.  Maple Code Editor Diagnostics says "No Errors" but Maple Code Editor Console says "Error: recursive assignment."  The MPL file is 485 lines but I am assuming the error could be anywhere in the file.  Is there any way to fix this?

Update: I found the error by copying just a handful of lines at a time into Maple until the error was reproduced.  It was a typo of just one character out of 16,960 characters.  But this is very frustrating.  Why can't Maple 2023 tell me what line the error is on?  Is there a better way to debug MPL files?