I have been using 'MemoryInUse()' to try to track down a memory leak that is crashing my Maple-based simulation. The problematic proc seems to be ProbDStarLNu which calculates the probablity of decay as a function of the angles of the decay product. It gets called a lot in acceptance rejection decay generation.

'MemoryInUse()' slows it down a lot. It also seemed to give weird results to me. I calculate the mem used on entry and then after the 4 trig functions have been called. The amount of mem 'in use' goes down after the trig functions are called. I was suspecting they were the source of the leak, but they seem almost an anti-leak.

This weird behavior caused me to write my own mem counting proc ('memUsed') based on summing over the results of 'mem:= kernelopts(memusage)'.

It is much faster. It gives nearly but not quite the same numbers. It still goes "backward" when the trig functions are called. Most of the mem growth comes from the algebra that follows the trig.

Why would that be?

I have included code for 'ProbDStarLNu'. It's pretty simple.

I have been having problem in solving this question because of the 'Whittaker'.

Can anyone please asist me how to integrate below problem with the 'Whitaker?

f := exp(-1.5*t)*(GAMMA(2)*t^(2-alpha)/GAMMA(3-alpha)-t^(1-alpha)/GAMMA(2-alpha)+t^2-t);
               /   (2 - alpha)        (1 - alpha)           \
               |  t                  t                 2    |
   exp(-1.5 t) |---------------- - ---------------- + t  - t|
               \GAMMA(3 - alpha)   GAMMA(2 - alpha)         /
int(f, t);
 

2.^(3.-1.*alpha)*3.^(alpha-3.)*(t^(-.5000000000*alpha)*3.^(-.5000000000*alpha)*2.^(.5000000000*alpha)*(alpha-3.)*(-2.+alpha)*exp(-.7500000000*t)*WhittakerM(-.5000000000*alpha, -.5000000000*alpha+.5000000000, 1.500000000*t)/(3.-1.*alpha)+t^(-.5000000000*alpha)*3.^(-.5000000000*alpha)*2.^(.5000000000*alpha)*(1.500000000*t-1.*alpha+2.)*(alpha-3.)*exp(-.7500000000*t)*WhittakerM(-.5000000000*alpha+1., -.5000000000*alpha+.5000000000, 1.500000000*t)/(3.-1.*alpha))/GAMMA(3.-1.*alpha)-1.*2.^(2.-1.*alpha)*3.^(-2.+alpha)*(t^(-.5000000000*alpha)*3.^(-.5000000000*alpha)*2.^(.5000000000*alpha)*exp(-.7500000000*t)*WhittakerM(-.5000000000*alpha, -.5000000000*alpha+.5000000000, 1.500000000*t)+t^(-.5000000000*alpha)*3.^(-.5000000000*alpha)*2.^(.5000000000*alpha)*(-2.+alpha)*exp(-.7500000000*t)*WhittakerM(-.5000000000*alpha+1., -.5000000000*alpha+.5000000000, 1.500000000*t)/(2.-1.*alpha))/GAMMA(2.-1.*alpha)+.1481481481+(0.9876543210e-1*(-6.750000000*t^2-9.*t-6.))*exp(-1.500000000*t)-(.2222222222*(-2.-3.*t))*exp(-1.500000000*t)

I want infinite series expansion of incomplete upper gamma function.

I clicked on the author's page link and was redirected to front page instead. Does this link work? Why does it not work for me?

How to define constants - MaplePrimes

I'm trying to use CodeTools[Profile] to investigate my memory leak and the effect or lack there of of 'forget' to improve it.

I copied my code from an example:

with(CodeTools[Profiling]);
 [Allow, Build, Coverage, GetProfileTable, Ignore, IsProfiled, 

   LoadProfiles, Merge, PrintProfiles, Profile, Remove, 

   SaveProfiles, Select, SortBy, UnProfile]

Profile([PickAngles(1.0, 1.0, 1.0)]);
Error, (in CodeTools:-Profiling:-Profile) unexpected argument(s)

The 'unexpected argument' is the 'result' from 'PickAngles'. I've tried with and w/o the '[]' but it does the same thing failing to recognize that 'PickAngles' is a procedure.

The example I used is here:

https://www.maplesoft.com/support/help/maple/view.aspx?path=CodeTools%2FProfiling%2FProfile procedure.

I did what to make of the subscript Profiles in the example so I used '[Profiles]' and that worked or seemed to.

Dear Users!

I hope you are doing well. In the attached file I want to convert the system of ODEs (attained the system against the value of M) into matrix form and need the matrices A, B and vector b. Remember the order of A, B and b vary as M vary. I am waiting for your kind response. Please take care and thanks

System_of_ODEs.mw

I can import maple, but I get an error when importing namespace/symbols:

import maple.namespace as mpl

The below is not properly rendering in my browser. You probably have to open the worksheet
 

Edit: Asking differently: Which symbol from the palettes can be used for equivalence (in a logic context)?

Edit2: Title changed. Original title: Why does Maple replace the equivalence symbol by an implication symbol

Download De_Morgans_law.mw

hi

I have a linear system with varibles trying to plot 3d the solutions x, y, z

here is my code: linear_var.mw

please any comment might help.

*******************************

restart;

Omega:=10:N:=0.5:M:=sqrt(N(N+1)):

a11:=0.5*(1+2*N)+M*cos(phi):
a12:=-0.5*((1+theta)^3+(1-theta)^3):a13:=-0.5*(N+M*cos(phi))*((1+theta)^3-(1-theta)^3): a21:=M*sin(phi): a22:=(-(1+2*N)+0.5*M*cos(phi))*((1+theta)^3+(1-theta)^3): a23:=-(Omega+0.5*((1+theta)^3-(1-theta)^3)*M*sin(phi)): a31:=0.25*((1+theta)^3-(1-theta)^3): a32:=Omega: a33:=-0.5-(N+0.25)*((1+theta)^3+(1-theta)^3): b1:=-0.5*a31: b2:=0: b3:=0.25+((1+theta)^3+(1-theta)^3)/8:

slove([a11*x+a12*y+a13*z=b1,a21*x+a22*y+a23*z=b2,a31*x+a32*y+a33*z=b3[,[x,y,z]);
Error, unable to match delimiters
Typesetting:-mambiguous(Typesetting:-mambiguous(slovelparlsqba11

  sdotx + a12sdoty + a13sdotzequalsb1commaa21sdotx + a22sdoty + 

  a23sdotzequalsb2commaa31sdotx + a32sdoty + a33sdotzequalsb3lsqb

  comma(xyz)rparsemi, 

  Typesetting:-merror("unable to match delimiters")))

plot3d(x, theta = .1 .. 5, phi = 0 .. 2*Pi, axes = boxed);
plot3d(y, theta = .1 .. 5, phi = 0 .. 2*Pi, axes = boxed); plot3d(z, theta = .1 .. 5, phi = 0 .. 2*Pi, axes = boxed);

I would like solving systems of congruences like the following one:

154x+69y = 0 mod 7^3

13x+716y = 0 mod 13^3

23x+3059y = 0 mod 23^3

I need the minimal non-trivial solution. I know that the solution of the system above is x=1103, y=26390.

How could I find the solution with Maple?

Hi!

I'm having problems with my maple not saving. I get no error message and no windows pop up.
The error occurred after the summer holidays in Maple 2022. I use Windows 11, everything is up to date. Have no Anti-virus programs.

  I have tried the following after I discovered the error:
- To install Maple 2023
- Uninstall and delete all maple folders, then reinstall Maple 2023
- Pressed Ctrl + s
- Press "Save as..."
- Pressed on the floppy disk/save icon
- Restarted computer and updated windows
- Run Maple as administrator

The only thing I can be allowed to do is:
- Ctrl + p
- Print to PDF/printer

Really hope you can help!

Dear power users, I am making the switch from Mathcad towards Maple and would like to know what is the most efficient alternative in Maple for a solve block. I have attached a work document to illustrate better my question.SolveBlockQuestion.mw Any help is highly appreciated.

Say I have ,

sol1 := C*vout*vin/(Iout*L*k^2)

and I want to simplify this expression by replacing  with .

Then I need to simplify it further by noting that , and that , so that the answer is only in terms of , , and .

I tried eval(sol1, vin/Iout = omega*Lm) initially to step through the first substitution, but that just returned the same expression.

Thank you.

I don't know where my last exchange with @mz6687  has been moved (not to the initial question  https://www.mapleprimes.com/questions/237066-Determinant-The-System-Hangs for what I see).
Nevertheless here is the reply I was sending to @mz6687  which ended with the message "Page not found".
Are_you_ok_with_that.mw

Could the one who moved the question meanwhile be so kind as to attach this reply?
TIA

If one writes Pi or Catalan, the output is Pi or Catalan, and writing evalf(Pi) or evalf(Catalan) we get approximations.

I would like to define other constants behaving in the same way. How to do it?