i wrote a procedure to  find some functions but at the end  of the procedure gives unexpected string error. i didn't find any way out.

Best regards. 

Download procedureA.mws

How would I programmatically open a maple worksheet?

I thought this would work but it brought up a -1

ssystem("c:\\Users\laptop22\Documents\Start2Test.mw")

I am trying to solve the wave equation in polar coordinates.  The initial condition on u is given by f(r,theta) and the initial condition on u_t is zero.  The weight function is w(r).  I am not sure why it will not evaluate this as I know the solution remains finite on the domain (the unit disk).  Here is the code: 
 

Download Section_6.3.mw

Any assistance would be greatly appreciated. 

Hallo

I'm trying to plot several periods of this 2π-periodic function f (x) = 3x, x ∈]-π; π] with the following code which I got from an older post.

However, when using uneven multipla of PI as x it just gives me a verticle line. Can anyone shed some light to this or perhaps give a better solution?

Thanks in advance

Best regards

Thomas

In paper of Nail H. Ibragimov there is operator given at equation 2.10 which he call by name Euler-Lagrange operator. How I can use this operator in Maple to derive adjoint equations ?
I have attached JPEG in which I want to apply this operator to Lagrange identity "I" for adjoint equations.

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,

Hi,

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

P:=(u,v)->[u,v,-(u^2+v^2)]

and the “sun”

Sun:=t->[5*cos(t),0,5*sin(t)]

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

NP.Nsun

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 (r,t) . The metric is:

ds^2=(1-rg/r)*dt^2-(1-rg/r)^(-1)*dr^2

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 and .  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

dAlembertian.mw

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.
DeltaMag:=0:
DeltaBond:=0:

# 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.
printlevel:=4:
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.
DeltaMag:=0:
DeltaBond:=0:

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?

Hi.

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

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

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

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

   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?