MaplePrimes Questions

Hi,

How to insert the notation of a vector, with an arrow? thanks

S5VecteurTextPlot.mw

evalf(Pi*int(2*299792458^2*662607015*10^(-8)*10^(-34)/((exp(299792458*662607015*10^(-8)*10^(-34)/(1380649*10^(-6)*10^(-23)*lambda*5772.0)) - 1)*lambda^5), lambda = 0 .. infinity));

returns a complex number for the Solar radiant exitance in watts per square meter with an effective blackbody temperature of "5772.0" kelvins.

Replacing "5772.0" with an exact "5772" gives the correct answer though.  The nominal Solar temperature is defined as exactly 5772 kelvins whereas the physical Solar effective temperature is 5772.0 K (+/-) 0.8 K.

Is there a way to get Maple 2023 to avoid outputting these types of errors when using nonexact quantities?

I tried replacing "5772.0" with "5772.00000000001" and "5771.9999999999" and the error is not duplicated for those kelvin temperatures.

Interestingly, this error does not occur on Maple for Excel, only within the Maple desktop client....

https://scienceworld.wolfram.com/physics/Stefan-BoltzmannLaw.html

https://iopscience.iop.org/article/10.3847/0004-6256/152/2/41/

I'm currently addressing a problem related to modified Bessel functions using an older version of Maple (the specific version escapes my memory). In an attempt to resolve issues, I've experimented with the trial version of Maple 2023, but I've encountered an unusual phenomenon. Expressions that were previously simplifiable in Maple now resist simplification. The specific expression provided below, which should equate to 1, fails to be recognized as such by Maple. This poses a concern as it could lead to overly complex expressions in subsequent steps, considering this expression is only an intermediate stage. Is there a recommended approach to overcome this challenge?

f := (BesselI(0, alpha)*alpha-2*BesselI(1, alpha))/(BesselK(0, alpha)*BesselI(1, alpha)*BesselI(0, alpha)*alpha^2+BesselK(1, alpha)*BesselI(0, alpha)^2*alpha^2-2*BesselI(1, alpha))

(BesselI(0, alpha)*alpha-2*BesselI(1, alpha))/(BesselK(0, alpha)*BesselI(1, alpha)*BesselI(0, alpha)*alpha^2+BesselK(1, alpha)*BesselI(0, alpha)^2*alpha^2-2*BesselI(1, alpha))

simplify(f)

(BesselI(0, alpha)*alpha-2*BesselI(1, alpha))/(BesselK(0, alpha)*BesselI(1, alpha)*BesselI(0, alpha)*alpha^2+BesselK(1, alpha)*BesselI(0, alpha)^2*alpha^2-2*BesselI(1, alpha))

eval(f, alpha = .25)

1.000000000

NULL

Download question.mw

1. Sorry the question is not a maple product but just wanted to know if anyone was using the sharp equation calculator with a pc link interface as it is not with a very acceptable web presence and the kit is unaffordable. If Maple has an equivalent I would love to know.

2. Does Maple have a facility for using virtual drives and emulators on a portable basis (non cloud)?

Hi Math People,

I can't find an answer if it is possible to insert slider or any other component into document, adjust its value by sliding it and execute all following calculations in the document below the slider, those which depends on the value of slider parameter?

I know this is possible in graphs, or explore command, but could it be simply setting some value let's say "a". Lets set its value by slider and when slider is set, code should automatically update document? This part of code of slider "Action When Value changes"

Or am I expecting too much?

Best regards

Marcin

how to write the elements of the triangle of Pascal so that all 1 are aligned on a same vertical. Thank you.

this is my model. Please give me how to find the DFE and basic reproduction number from maple.

I am a little overwhelmed by the possibilities of ColorTools.
Is there a graphical overview of implemented palettes side by side?

please help me to fixed and find the analytic solution...

restart

with(linalg)

f1 := mu[1]*N-delta*r-beta*N*s(e+i)

mu[1]*N-delta*r-beta*N*s(e+i)

(1)

f2 := beta*N*s(e+i)+omega*v-(mu[1]+pi+gamma)*e

beta*N*s(e+i)+omega*v-(mu[1]+pi+gamma)*e

(2)

``

f3 := pi*e-(mu[1]+mu[2]+sigma)*i

pi*e-(mu[1]+mu[2]+sigma)*i

(3)

f4 := gamma*e+sigma*i-(theta+mu[1]+mu[2])*q

gamma*e+sigma*i-(theta+mu[1]+mu[2])*q

(4)

f5 := theta*q+xi*v-(mu[1]+delta)*r

theta*q+xi*v-(mu[1]+delta)*r

(5)

f6 := eta*N*s-(mu[1]+omega+xi)*v

eta*N*s-(mu[1]+omega+xi)*v

(6)

T := solve({f1, f2, f3, f4, f5, f6}, [s, e, i, q, r, v])

Error, (in solve) cannot solve for an unknown function with other operations in its arguments

 

``

Download cobacoba.mw

f := x -> 2*x/(x + 3)

p1 := plot(f(x), x = 0 .. 100);

p2 := with(plots);
inequal(58 < x, x = 0 .. 100, y = 0 .. 2);

display([p1, p2]);

Error, (in plots:-display) expecting plot structures but received:

How do i combine to p1 and p2 in the same graph. They print well individually..  Its maple 2021

Hello

I could not come up with a better title, so I apologize for that.  Let me explain what I am trying to implement using an example. 

Consider the following list of monomials (a 3D model,  4 monomials altogether):

model := [x^2*y*alpha[1, 11], x*z^2*alpha[2, 15], y^2*z*alpha[3, 17] + y*z*alpha[3, 8]]

The idea is to generate the set of all models by adding a non-repetitive monomial from the following list to model.

[alpha[i, 0], alpha[i, 1]*x, alpha[i, 2]*y, alpha[i, 3]*z, alpha[i, 4]*x^2, alpha[i, 5]*y*x, alpha[i, 6]*z*x, alpha[i, 7]*y^2, alpha[i, 8]*z*y, alpha[i, 9]*z^2, alpha[i, 10]*x^3, alpha[i, 11]*y*x^2, alpha[i, 12]*z*x^2, alpha[i, 13]*y^2*x, alpha[i, 14]*z*y*x, alpha[i, 15]*z^2*x, alpha[i, 16]*y^3, alpha[i, 17]*z*y^2, alpha[i, 18]*z^2*y, alpha[i, 19]*z^3]

where i indicates the coordinate where the monomial will be included.  Example: if alpha[i, 2]*y is to be added to the second coordinate of model it goes as alpha[2,2]*y. Note that alpha[i, 15]*z^2*x cannot be added to coordinate 2 since it is already there.  

The result will be a list of 56 models with 5 monomials. 

How can I do that efficiently?  

Many thanks

Ed.

PS. I have implemented something similar to the problem above, but I have used too many 'for' loops. 

I have a result:

u1 = RootOf(4*_Z^2 + (4*RootOf(60*_Z^3 - 60*_Z^2 + 15*_Z - 1) - 4)*_Z + 4*RootOf(60*_Z^3 - 60*_Z^2 + 15*_Z - 1)^2 - 4*RootOf(60*_Z^3 - 60*_Z^2 + 15*_Z - 1) + 1)

I tried: convert(u1,'radical');

but it reported error Error, (in sprintf) too many levels of recursion

how can I convert it to this form?

Thanks!!

Hello! We are using Maple to evaluate numerically expressions with very high precision and we are finding very strange results. After some playing with simple models we have come to conclusion that there might be a bug in Maple 2023. Here is a clean example of what we are experiencing. Plotting the function x*exp(-x) in [0,20]. Clearly the function cannot be negative in that interval! However, if Digits is set to 40 the results may be whatever. Any ideas? Thank you very much. 

Trying to figure out how to save a slice of an array to a repository.

for example: savelib(vmcofs[1],lib) [where lib is my repository location]

I'm getting the error: Error, invalid input: savelib expects its 1st argument, to be of type symbol, but received vmcofs[1]

Quoting it doesn't seem to help.

I need to do this programmatically as I will be incremented the row number in the array and saving partial data at each step.

That is...savelib(vmcofs[1],lib) then savelib(vmcofs[2],lib), etc.

I can save the whole array. Is there a more efficient way of doing this?

Hello everyone!

I'm using Maple 18. The problem is that fsolve does not solve equations, or more precisely, it does not pass the values of sU - from its range, and sV - from its range, to shootNL. As a result, R value:=[]; . I would like to ask for help in solving this problem.

I have attached the Maple source file.

fsolveProblem.mw

First 137 138 139 140 141 142 143 Last Page 139 of 2427