Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Strange behavior. I don't know if this can be made to work.  I replaced an equation with an interpolation object, f(x).

tInt := int(f, 133  .. 134, numeric);
tInt := 324.95206944305  <= Computes a result. I

tInt := int(f, 133  .. 312, numeric);
tInt:= int(fproc, 133. .. 312.) <= Won't compute it. The interpolation object has 3051 points, so that's not the problem. The name "fproc" is not something I have in my  code, just the name, "f".

If I try to make the  integral just slightly more complex by  passing a variable in, it doesn't work at all. This was to see if I could integrate with a varying divisor: f(x)/(x+1)
int(f(x), x = 133 .. 134, numeric); 

Error, (in int) wrong number (or type) of arguments: invalid options or option values passed to definite integration. <=

A fallback  was to use summation, but that doesn't work at all, same typ of error as the Error for int. (Not shown.)

So my fallback is to use a for loop to sum the set. The for loops produce another set of values, and I think I have to make interpolation objects for that as well. And so on. I really  was expecting  Maple to handle this better.

conjcl := proc(k, p, n) local t, L, q; L := []; t := 0; q := p^n - 1; while k*p^t <> k mod q do t := t + 1; L := [op(L), t mod q]; end do; return L; end proc

I'm fairly new to coding with Maple / in general, but basically for my course I need to do the following:
For each t=0,1,... I need to first check if kp^t = k mod (p^n-1), and if it isn't, I compute t mod (p^n-1) then add it to my list. The problem I have is I just can't seem to get it to work. I incliuded both the screenshot of the code and typed it out. Thanks in advance

Good day sirs,

I am having trouble identifying the error in the system of first-order ODEs. It gives the error code of

" Error, (in dsolve/numeric/process_input) received more than one indication of the dependent vars [theta[1](t), theta[2](t), theta[3](t), theta[4](t), theta[5](t), theta[6](t), theta[7](t), theta[8](t), theta[9](t)].

The document is attached below.

Thank you
"  Systems_error.mw

The attached worksheet contains expressions that cause a warning message to appear in my worksheet.  In another post it was suggested that the variable i should be explicitly declared local in the definition of the arrow expression. However, as shown in the same worksheet, this gives rise to an error message. Subsequently, I found a statement in help that suggested this approach would not work in 2-d math notation.  Are there any other ways of eliminating these warnings in 2-d math notation?

warningmessage.mw

Dear all

I run my code, I think evrything is well coded, but I get the following error 

invalid left hand side in assignment

Please, I need a help to solve this problem,  I haven't any idea about the origin of this error 

IOMMcode.mw

Thank you 

In a Maple session I have several worksheets open. When I click on the close cross of a worksheet tab the following can happen:

Closing the worksheet next to it works as normal

 

This occurs irregularly. For far I could not make it reproducible to report it to Maplesoft support. Anyone has seen something like this or ideas what the cause could be? 

Is this new in Maple 2022.2 or was it there before? I have no way now to verify. Could someone please check? I am not using assuming anywhere, so this must be internal.
 

interface(version);

`Standard Worksheet Interface, Maple 2022.2, Windows 10, October 23 2022 Build ID 1657361`

restart;

ode := diff(f(Z), Z$2) = 3602879701896397*(3860741894751515/1125899906842624 + (5048392920194975*f(Z)^2)/1073741824 - (2484212754397225*f(Z)^4)/1073741824 - (321579057930643*f(Z)^6)/2147483648 - (4936153155163025*f(Z)^8)/2251799813685248)/(4611686018427387904*f(Z)^3*(4178268760908915/1073741824 + 315761000*f(Z)^2));
dsolve([ode,f(0)=1,D(f)(0)=0])

diff(diff(f(Z), Z), Z) = (3602879701896397/4611686018427387904)*(3860741894751515/1125899906842624+(5048392920194975/1073741824)*f(Z)^2-(2484212754397225/1073741824)*f(Z)^4-(321579057930643/2147483648)*f(Z)^6-(4936153155163025/2251799813685248)*f(Z)^8)/(f(Z)^3*(4178268760908915/1073741824+315761000*f(Z)^2))

Error, (in dsolve) when calling 'assume'. Received: 'the assumed property or properties cannot be satisfied'

 


Maple 2022.2 on windows 10

Download dsolve_nov_10_2022.mw

This integral produces very large result and division by zero when simplified. I do not have earlier version of Maple to check if this is a new problem or not

f:=1/(3*u^(2/3)+3*((-4*u^(1/3)+1)*u^(4/3))^(1/2)-12*u);
anti:=int(f,u);
MmaTranslator:-Mma:-LeafCount(anti);
simplify(anti)

For reference, this is Mathematica's results with leaf count only 35 instead of 41,385 and no problem when simplifying.

f=1/(3*u^(2/3)+3*((-4*u^(1/3)+1)*u^(4/3))^(1/2)-12*u)
anti=Integrate[f,u]
LeafCount[anti]
Simplify[anti]

Could someone please check if this fails same way on earlier Maple version?

Hi,

I am fairly new to Maple and trying to figure out how to output text to a file, but nothing seems to be working. I would like to generate code files for Julia using Maple expressions.

I am in worksheet mode, using a 'code edit' region. If I do:

jlfile := fopen("./juliafunc.jl",WRITE,TEXT);

then the file is created, as intended. However, neither of the following commands seem to write any text to the file:

writeline(jlfile,"This is a test.");
fprintf(jlfile,"This is a test.\n");

(these are similar to the examples in the dcumentation)

Is there something I am doing wrong? I am using Maple 2020 on Linux.

Maple help stopped working, when I call the help from the menu or the corresponding button in the toolbar, the help is not displayed. What I can do?

Hi everyone, I have a question regarding the time in a plot.

The idea behind is, that i dont have to adjust the time after every change of parameters.

I have a function that only is valid until a certain time. I evaluete the time until cut off for example like this:

tempLrest := Optimization:-Minimize(tempLrest(t), t = 0 .. 2);
              tempLrest := [2.47883089786094 10^-8  , [t = 0.414486068771508]]

cutoff := tempLrest[2];

Now i want to call "cutoff" in an odeplot like this: 

odeplot(sol1, [t, z(t)], 0 ..cutoff , size = [300, 300])

unfortunately this doesnt work. 

Does everyone know how to solve this problem?

Thank you in advance!

Hello everybody. My lovely Dutch math book is progressing. Last week ive been scanning a copy of Advanced Engineering Mathematics from the library, a grueling task which took me 24 full hours to complete. But ive got a book that should be able to teach it with the right mindset and social contract (this Dutch book is too distracting, i cant even think). I will finish it though.. The sturdy that stands i guess. 

Any way the problem here is that the side from which a line interfaces with a shape in the 3D plane only has an x component and is lacking a y and z component. So the calculations are impossilbe, because you cant devide through 0. I suspect its something simple, but like said the repertoire is so terrible it feels like im scrambled eggs if i even remotely think about solving the problem. The facade makes the man right? Feels more like one of these spagetti westerns with only a frontage without anything behind it. 

The Question is as follows (translated from Dutch): "Given is the half ball with the equation z=sqrt(4-x^2-y^2). Check that the given point lays on the half ball and determin the eqation of the interface with the half ball in the given point."

So basically have a vector from that given point that will hit the surface of the half ball.   
 

#Opdracht 2

#a.

restart

z := sqrt(-x^2-y^2+4)

(-x^2-y^2+4)^(1/2)

(1)

"f(x,y):=sqrt(4-x^2-y^2)"

proc (x, y) options operator, arrow, function_assign; sqrt(-x^2-y^2+4) end proc

(2)

a := f(1, 1, sqrt(2))

2^(1/2)

(3)

partafgx := D[1](f)

proc (x, y) options operator, arrow, function_assign; -x/(-x^2-y^2+4)^(1/2) end proc

(4)

b := partafgx(1, 1, sqrt(2))

-(1/2)*2^(1/2)

(5)

c := (D[2](f))(1, 1, sqrt(2))

-(1/2)*2^(1/2)

(6)

raakvlak := z-a = b*(x-1)+c*(y-1)

z-2^(1/2) = -(1/2)*2^(1/2)*(x-1)-(1/2)*2^(1/2)*(y-1)

(7)

evalc(%)

z-2^(1/2) = -(1/2)*2^(1/2)*x+2^(1/2)-(1/2)*2^(1/2)*y

(8)

z = solve(raakvlak, z)

z = -(1/2)*2^(1/2)*x+2*2^(1/2)-(1/2)*2^(1/2)*y

(9)

#b.

restart

"f(x,y):=sqrt(4-x^2-y^2)"

proc (x, y) options operator, arrow, function_assign; sqrt(4-x^2-y^2) end proc

(10)

a := f(2, 0, 0)

0

(11)

partafgx := D[1](f)

proc (x, y) options operator, arrow, function_assign; -x/(4-x^2-y^2)^(1/2) end proc

(12)

b := partafgx(2, 0, 0)

Error, (in partafgx) numeric exception: division by zero

 

c := (D[2](f))(2, 0, 0)

Error, (in unknown) numeric exception: division by zero

 

raakvlak := z-a = b*(x-2)+c*y

z = b*(x-2)+c*y

(13)

z = solve(raakvlak, z)

z = b*x+c*y-2*b

(14)

#ik krijg er niets uit wat ergens op lijkt

#c.

restart

"f(x,y):=sqrt(4-x^2-y^2)"

proc (x, y) options operator, arrow, function_assign; sqrt(4-x^2-y^2) end proc

(15)

a := f(sqrt(2), sqrt(2), 0)

0

(16)

partafgx := D[1](f)

proc (x, y) options operator, arrow, function_assign; -x/(4-x^2-y^2)^(1/2) end proc

(17)

b := partafgx(sqrt(2), sqrt(2), 0)

Error, (in partafgx) numeric exception: division by zero

 

c := (D[2](f))(sqrt(2), sqrt(2), 0)

Error, (in unknown) numeric exception: division by zero

 

raakvlak := z-a = b*(x-sqrt(2))+c*sqrt(2)

z = b*(x-2^(1/2))+c*2^(1/2)

(18)

z = solve(raakvlak, z)

z = -b*2^(1/2)+c*2^(1/2)+b*x

(19)

``

Thank you very much. The first answer "a" was right, but for the other ones, this method they gave does not work. 

Greetings,

The Function 

Download Mapleprimes_Question_Book_2_Paragraph_5.5_Question_2b_2c.mw

After this code:

g1 := PermutationGroup({[[1, 2]], [[1, 2, 3], [4, 5]]})

g1:=<(1, 2), (1, 2, 3)(4, 5)>

I can get a string "<(1, 2), (1, 2, 3)(4, 5)>". Could I use the string to restore the PermutationGroup with an elegant method?

Here is a worksheet with a system of differential equations

Maple_DE_Example.mw

The system is

y'(x)=z(x)*g(x)

z'(x)=y(x)*g(x)

where g(x) is a random number generator seeded by x. That is, for the same x, g gives the same result. 

Using the classical method as an option to dsolve, a solution is plotted. Using rkf45, however, there is pretty much always an error that says

Warning, cannot evaluate the solution further right of ...

Why does it work with classical but not rkf45?

Is there a command that will count the number of digits in a binary number or any number for that matter? For example: 10110  has 5  digits in binary representation.

restart

I choose to define a number in base 10.

number := 22

22

(1)

Then I choose to convert that number into base 2.

number_binary := convert(number, binary)

10110

(2)

On visual inspection it is determined that number_binary has 5 digits in binary representation.

Is there a common command that can count the number of digits in a binary number? Such as, 10110 has 5 digits.

Download Count_digits_in_a_binary_number.mw

Count_digits_in_a_binary_number.mw

First 141 142 143 144 145 146 147 Last Page 143 of 2097