The Joint Mathematics Meetings are taking place this week (January 6 – 9) in Seattle, Washington, U.S.A. This will be the 99th annual winter meeting of the Mathematical Association of America (MAA) and the 122nd annual meeting of the American Mathematical Society (AMS).

Maplesoft will be exhibiting at booth #203 as well as in the networking area. Please stop by our booth or the networking area to chat with me and other members of the Maplesoft team, as well as to pick up some free Maplesoft swag or win some prizes.

Given the size of the Joint Math Meetings, it can be challenging to pick which events to attend. Hopefully we can help by suggesting a few Maple-related talks and events:

Maplesoft is hosting a catered reception and presentation ‘Challenges of Modern Education: Bringing Math Instruction Online’ on Thursday, January 7th at 18:00 in the Cedar Room at the Seattle Sheraton. You can find more details and registration information here: www.maplesoft.com/jmm

Another not to miss Maple event is “30 Years of Digitizing Mathematical Knowledge with Maple”, presented by Edgardo Cheb-Terrab, on Thursday, January 7 at 10:00 in Room 603 of the Convention Center.

Here’s a list of Maple-related events and talks:

Exploration of Mathematics Teaching and Assessment through Maple-Software Projects of Art Diagram Design as Undergraduate Student Research Projects

Wednesday, Jan 6, 10:20, Room 2B, Convention Center

Lina Wu

30 Years of Digitizing Mathematical Knowledge with Maple

Thursday, Jan 7, 10:00, Room 603, Convention Center

Edgardo Cheb-Terrab

MAA Poster Session – Collaborative Research: Maplets for Calculus

Thursday, Jan 7, 14:00, Hall 4F, 4th Floor, Convention Center

Challenges of Modern Education: Bringing Math Instruction Online

Thursday, Jan 7, 18:00, Cedar Room, 2^nd Floor, Sheraton Center

Using Maple to Promote Modelling in Differential Equations

Friday, Jan 8, 10:40, Room 617, Convention Center

Patrice G Tiffany; Rosemary C Farley

If you are presenting at Joint Math and would like to advertise your Maple-related talk, please feel free to comment below, or send me a message with your event and I’ll add it to the list above.

See you in Seattle!

Daniel

Maple Product Manager

http://www.mapleprimes.com/posts/38019-Calling-Out-To-C-From-Maple#

if i can use maple to call c# function such as AForge.QLearning

how to set some tasks for it to guess some system of polynomials to fit hibert series criteria?

how to set a game for it to run itself to discover itself?

Hello,

I have a worksheet that makes 38 Mo. I find it a bit huge because my code is only composed by 20 code lines. I guess that Maple has saved a important volume of data which makes the file very heavy.

I would like to reduce the size of my worksheet so as to be able to send it to colleagues.

How can I remove all the data save in a file?

I find how to display only the inputs but the data in term of size seems to be still present.

Thanks a lot for your help

what is the reason for negative value is displayed on the table? 

thank you for helping>>>>>mapel_prime.mw

Good evening,

I am trying to solve the following sum equation:

sum((factorial(n)/factorial(2*n))^n, n = 1 .. infinity)

but I require the step-by-step solution as to how it is done, but can't seem to find that option. 

If anone could help me out, that would be great.

Kind regards,

A

Dear, need help about Implicit finite difference method 

plz reply as soon as possible

Thanks

Hi,

Need help in plotting the function in the attached file. m, h and k are parameters.

Thanks

Bessel_Function.mw

> restart;
> u := -3*sqrt(mu)*tanh*(A+sqrt(-mu)*(x-(1/2)*mu*t))/sqrt(6*a)+3*sqrt(mu)*sech*(A+sqrt(-mu)*(x-(1/2)*mu*t))/sqrt(-6*a);

>
> mu := .5;
0.5
> lambda := 1;
1
> a := 0.5e-1;
0.05
> A := 1.5;
1.5
> plot3d([abs(u)], x = -3 .. 3, t = -3 .. 3);

Hello guys and gals!

I'm not strong enough with maple to get what the result I want.

It seems that it's because I'm asking for two lenths, and not a lenth and an angle, but I have no Idea how to tackle it diferently.

If you know a trick, please share it!

Here's an image:

http://imgur.com/xavAUoB

And here's the maple file attached (I think)

 complex_problem_from_the_internet.mw

Thanks,

Happy new year!

This post is my attempt to answer the question from here .  

The procedure  ContoursWithLabels  has 2 required parameters: Expr  is an expression in  x  and  y  variables,  Range1  and  Range2  are ranges for  x  and  y . In this case, the output is the list of floats for the contours and 8 black contours (with labels) (the axis of coordinates as a box). 

The optional parameters: Number is positive integer - the number of contours (by default Number=8),  S is a set of real numbers  C  for contours (for which Expr=C) (by default  S={}),  GraphicOptions  is a list of graphic options for plotting (by default  GraphicOptions=[color = black, axes = box]),  Coloring  is an equality  Coloring=list of color options for  plots[dencityplot]  command (by default Coloring=NULL). 

The code of the procedure:

restart;

ContoursWithLabels := proc (Expr, Range1::(range(realcons)), Range2::(range(realcons)), Number::posint := 8, S::(set(realcons)) := {}, GraphicOptions::list := [color = black, axes = box], Coloring::`=` := NULL)

local r1, r2, L, f, L1, h, S1, P, P1, r, M, C, T, p, p1, m, n, A, B, E;

uses plots, plottools;

f := unapply(Expr, x, y);

if S = {} then r1 := rand(convert(Range1, float)); r2 := rand(convert(Range2, float));

L := [seq([r1(), r2()], i = 1 .. 205)];

L1 := convert(sort(select(a->type(a, realcons), [seq(f(op(t)), t = L)]), (a, b) ->is(abs(a) < abs(b))), set);

h := (L1[-6]-L1[1])/Number;

S1 := [seq(L1[1]+(1/2)*h+h*(n-1), n = 1 .. Number)] else

S1 := convert(S, list)  fi;

print(Contours = evalf[2](S1));

r := k->rand(20 .. k-20); M := []; T := [];

for C in S1 do

P := implicitplot(Expr = C, x = Range1, y = Range2, op(GraphicOptions), gridrefine = 3);

P1 := [getdata(P)];

for p in P1 do

p1 := convert(p[3], listlist); n := nops(p1);

if n < 500 then m := `if`(40 < n, (r(n))(), round((1/2)*n)); M := `if`(40 < n, [op(M), p1[1 .. m-11], p1[m+11 .. n]], [op(M), p1]); T := [op(T), [op(p1[m]), evalf[2](C)]] else

if 500 <= n then h := floor((1/2)*n); m := (r(h))(); M := [op(M), p1[1 .. m-11], p1[m+11 .. m+h-11], p1[m+h+11 .. n]]; T := [op(T), [op(p1[m]), evalf[2](C)], [op(p1[m+h]), evalf[2](C)]]

fi; fi; od; od;

A := plot(M, op(GraphicOptions));

B := plots:-textplot(T);

if Coloring = NULL then E := NULL else E := ([plots:-densityplot])(Expr, x = Range1, y = Range2, op(rhs(Coloring)))  fi;

display(E, A, B);

end proc:

Examples of use:

ContoursWithLabels(x^2+y^2, -3 .. 3, -3 .. 3);

ContoursWithLabels(x^2-y^2, -5 .. 5, -5 .. 5, {-20, -15, -10, -5, 0, 5, 10, 15, 20}, [color = black, thickness = 2, axes = box], Coloring = [colorstyle = HUE, colorscheme = ["White", "Red"], style = surface]);

The next example, I took from here:

ContoursWithLabels(sin(1.3*x)*cos(.9*y)+cos(.8*x)*sin(1.9*y)+cos(.2*x*y), -5 .. 0, 2 .. 5, {seq(-2 .. 2, 0.5)}, [color = black, axes = box], Coloring = [colorstyle = HUE, colorscheme = ["Cyan", "Red"], style = surface]);

There are many more examples can be found in the attached file. 

ContoursWithLabels1.mw

Edit. The attached file has been corrected.

Hi,

  I think similar question has been asked by several people, but I did not find a suitable thread. My question is, suppose I have a probablity distirubtion function like

  p(x,y) = exp(-alpha (x+y) ) x^2 y^2 / |x-y|  , alpha>0

 x,y goes from - \infty to + \infty. This function is normalizable but unbounded, which makes the rejection algorithm a bit difficult(?).

  How to generate samping points from this type of probability distribution function?

Thank you very much!

hi. i am tottaly new to the maple and i have a problem.

consider function f with variables x & y which are not independent. x & y are functions of t and the relation is unknown.

for example i wanna to the below job:

define f as f=x^2+y^2

differentiate it with respect to t : diff(f,t) which should give me 2(dx/dt)x+2(dy/dt)y

i've googled it alot and i couldn't find anything usefull.

(the problem is how to set x as a function of t with unknown relation and use it in another function and then differentiate it with respect to t)

thanks alot :)

has anybody experienced issues with Maple connectiong to Solidworks?  There doen't seem to be much info ion trouble shooting a failed "OpenConnection" call.

thanks,

Bill

From a Maple Primes answer two years ago:

f(x,y) is the equation of a line through point [m,n]. The solve command finds values of a and b for which f(x,y) are lines through [m,n] and tangent to x^2 + y^2 = r^2.

f := proc (x, y) options operator, arrow; a*(x-m)+b*(y-n) end proc

solve([f(0, 0) = r, a^2+b^2 = 1], [a, b])

These commands are far from the conventional solution. Why do they provide the correct answers?

Hello,

I would like to plot gait diagrams (the lines you can see on the picture belowà from the solutions obtained with a NL oscillator (composed with 8 coupled odes). Here the result that I would like to obtain.

Initial plot:

Desired plot

I would like to obtain 4 lines corresponding to the 4 elliptic trajectories obtained with the NL oscillator. The four lines should be done like this. When the trajectory is above 0, the line should be colored in green. When the trajectory is below 0, the line should be colored in black. 

May you help me to define this kind of graph called gait diagrams from the solution of the NL oscillator ?

Here you can find my maple code:

K:=Matrix([<0, -1, 1, -1>,<-1, 0, -1, 1>,<-1, 1, 0,-1>,<1, -1, -1,0>]);

for i to 4
do
r[i]:=sqrt((u[i](t))^2+(v[i](t))^2):
omega[i]:=omega[sw]/(1+exp(b*v[i](t)))+omega[st]/(1+exp(-b*v[i](t))):
Equ[i]:=diff(u[i](t),t)=Au*(1-r[i]^2)*u[i](t)-omega[i]*v[i](t):
Eqv[i]:=diff(v[i](t),t)=Av*(1-r[i]^2)*v[i](t)+omega[i]*u[i](t)+MatrixVectorMultiply(K,<seq(v[i](t),i=1..4)>)[i]:
EqSys[i]:=[Equ[i],Eqv[i]]:
end do:

paramsCycle:=omega[st]=4*2*Pi,omega[sw]=2*Pi,Au=5,Av=50,b=100;
params:=paramsCycle;

Differential system 
sys:=map(op,eval([seq(EqSys[i],i=1..4)],[params]));
ic:=[u[1](0)=0, v[1](0)=0,u[2](0)=0, v[2](0)=-0.1,u[3](0)=0, v[3](0)=0.1,u[4](0)=0, v[4](0)=0.1];
Résolution1
res:=dsolve([sys[],ic[]],numeric):
Initial boundaries
tcalc:=4;
ic2:=[seq(u[i](0)=eval(u[i](t), res(tcalc)),i=1..4),seq(v[i](0)=eval(v[i](t), res(tcalc)),i=1..4)];
Résolution2
res:=dsolve([sys[],ic2[]],numeric):

tmax:= 40:
numpts:=100*tmax:
plots:-odeplot(res,[t,v[1](t)],0..tmax,thickness=2, view=[0..5, -1.5..1.5],numpoints = numpts);
plots:-odeplot(res,[t,v[2](t)],0..tmax,thickness=2, view=[0..5, -1.5..1.5],numpoints = numpts);
plots:-odeplot(res,[t,v[3](t)],0..tmax,thickness=2, view=[0..5, -1.5..1.5],numpoints = numpts);
plots:-odeplot(res,[t,v[4](t)],0..tmax,thickness=2, view=[0..5, -1.5..1.5],numpoints = numpts);
plots:-odeplot(res,[seq([t,v[i](t)+i*5],i=1..4)],0..tmax,thickness=2,view=[0..5,0..25], numpoints = numpts);

Thanks a lot for your help