MaplePrimes Questions

I have to generate a code for carrying out the matrix form of the revised simplex method. I have a code in place but am struggling to convert the constraints into canonical form and introduce the penalty function. If anyone has any ideas I'd be very grateful!

Best Regards

Hello everyone,

The problem resides in the argument "l" (small l) in the function "ratio", more precisely in the function "expi" then "riskexpjs" then in "ratio", so after try-outs I noticed that, then maple couldn't execute the function plot as decribed in the script.

Best regards,


I am trying to animate a ”shadow” on an object in Maple. I have parametrized the following surface


and the “sun”


Ant thus evaluating the “shadow” to be where the dotproduct


Is less the zero, i.e NP.Nsun<0. Where I get the following:


My question is now, am I able to animate the shadow in a likewise manner as shown below

Animate( plot3d, [P(u,v),v=a..b,u=c..d,color=[f(u,v,t)<0]],t=0..T)?

thus using an implicit function as a colorfunction?

I have to prove the following:

So I do not need the explicit derivative of the function Psi(r,t) . The metric is:


I am in the case of a collapsing star that emit radiation during the collapsing.  And I do not need to have a rotating black hole so that the reason I dont have dt*dr term in the metric, and I fix theta and phi.  So if you look in the Maple file attach to this post, I don't manage to obtain what I need to prove the equality between the two aspect of the same calculation.

Plese, take into account that I am sort of novice with the Physcis package and that the question is not part of an exam.

Thank you in advance for your help. 

Mario Lemelin



I am trying to apply window functions from the SignalProcessing package to arrays. However, some window functions (e.g. Hann) appear to operate on a [0..N-1] index basis, others (e.g. Welch) work on [1..N].

I don't know how to make the latter ones work correctly, since after applying the windowing function the first entry in the array is not zero, as it should be.

Best regards

So below is a calculation im trying to create, but it just tells me there is an unexpected IF error, and returns me to the third For line. I have tried it straight and as a procedure but i just cant seem to get it to work. Any ideas or tweaks to make this work would be much appreciated. All of the other variables I've named work in their respective lines of code, which I have left out as it is long as it is!

ComlexAlgorithm:= proc(L)
local N, j, kf, o, DeltaMag, DeltaBond;
global NConfig, Configz1, Col, Row, OldConfig, NewConfig, NewTable;

# Set Delta values to zero.

# Creates loops to loop through the correct changes in values i.e. Add -1, then 1 to first point, move across each point one by one, repeat from the beginning for each configuration then repeat whole process for the amount of rows.
for N from -1 to 1 by 2 do
   for j from 1 to Col do
      for w from 0 to NConfig-1 do   
         for o from 1 to Row do

# Calculate DeltaMag.
if N=-1 then DeltaMag:=DeltaMag+1 else DeltaMag:=DeltaMag-1
end if:

# Calculate all relevant bonds i.e. 3 values for DeltaMag, except at "corner" lattice points where there is only 2 bonds.
if N= Configz1[w](o, j) then DeltaBond:= DeltaBond-1 else DeltaBond:= DeltaBond+1
end if:

if j-1>1 then
   if N <> Configz1[w](o, j+1) then DeltaBond:= DeltaBond+1
      if N = Configz1[w](o, j+1) then DeltaBond:= DeltaBond-1
else DeltaBond:= DeltaBond
end if:
end if:
end if:

if j+1<Col then
   if N <> Configz1[w](o, j+1) then DeltaBond:= DeltaBond+1
      if N = Configz1[w](o, j+1) then DeltaBond:= DeltaBond-1
else DeltaBond:= DeltaBond
end if:
end if:
end if:

# Calculate if, and by what degree, conifguration number changes using 2^(j-1) which is the arithmetric series for binary. Changes only occur when "incoming" spin is different.
if N= -1 then
   if Configz1[w](o, j) <> N then NewConfig[w]:= OldConfig[w-(2^(j-1))] 
else NewConfig[w]:= OldConfig[w]
end if:
end if:

if N= 1 then
   if Configz1[kf](o, j) <> N then NewConfig[w]:= OldConfig[w+(2^(j-1))] 
else NewConfig[kf]:= OldConfig[w]
end if:
end if:

NewTable[kf]:= [NewConfig[w], BondEnergy(w)+DeltaBond, MagEnergy(w)+DeltaMag];

# Set Delta values to zero to finsh.

end do:
end do:
end do:
end do:
end proc:

From the help page it says

To create a custom start page with a similar layout as the default start page, create a table in the worksheet, add a shortcut component to each table cell and modify each shortcut component appropriately.

I'm unable to find how to insert a shortcut component.  Can someone help?


Lets say im given 2 surfaces like this:

S1: z=x^2+y^2 and S2: z=-2x

To find all the points on this surface i set z=z and obtain the equation, but lets say I was given an arbutary vectorfield and was required to find the outward flux I would need and parametric equation for the surface which is now a solid enclosed by 2 surfaces.


Of course I could apply the divergence theoreme and integrate the divergens over the entire volume using tripple integral. This I suppose this requires a bit more work compared to using the vector calculus package.


Would anyone give me an hint on ho to acheive this? Im also looking to plot the object crated by the bounderies of the divergenc theorem bounderies

Int(Int(Int(4*z, z = x^2+y^2 .. -2*x), y = -sqrt(-x^2-2*x) .. sqrt(-x^2-2*x)), x = -2 .. 0)

I just happened to be searching older posts and could not find this one I created

  Why we buy - The science of shopping
September 18, 2009

It must have be zapped by a spamming filter or something. 

**edit added**  heres another I can't seem to find

Numb3rs, Maple and Wolframs Mathematica
July 23, 2009

   Consider a surface of the unit sphere about the origin, intersecting the volume of the cube [-a..a,-a..a,-a..a] 0<a<1.  What is the best Maple method to plot this surface with a direct mapping, that is without use of implicit plotting over a volume or rejecting grid points?

so yep pretty self explainatory, i was just wondering why it works when i use the asymptomatic expansion but not the iterated reciprocal substitution code shown in the that help interface for asympt.






subs(y = 1/y, series(subs(y = 1/y, Zeta(1/2+I*y)), y = 0, 6))

Error, (in series/Zeta) unable to compute series


asympt(Zeta(1/2+I*y), y, 6)







Hope everything going fine with you. I am facing problem to fine the exact (numerical) solution of the attached system of linear PDEs associated with BSc and ICs. I tried to solve it without BCs and ICs, with BCs and with ICs also all the times I failed. Please solve it either general, with ICs or BCs. You can try to solve it numrically. In attached file H(t) represent the unit step function. I am waiting your positive response. 

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

hi i can rewrite answer in compressive format. that result seen readible and incomplex.for example how i can use change of variable, in order to final result has a short format...


Why would firewalls cause issues with connecting to the maple kernel if the CAS is entirely self contained (as ive been told before on this forum)

4 5 6 7 8 9 10 Last Page 6 of 1392