Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi,

Seem to be a bit stuck. Here's my code:
 

Thanks in advance :-) 

It is suggested  

hypergeom([1/3, 2/3], [3/2], (27/4)*z^2*(1-z)) = 1/z

if z > 1. Here is my try to prove that with Maple:


 

a := `assuming`([convert(hypergeom([1/3, 2/3], [3/2], (27/4)*z^2*(1-z)), elementary)], [z > 1])

-(1/((1/2)*(27*z^3-27*z^2+4)^(1/2)+(3/2)*z*(3*z-3)^(1/2))^(1/3)-1/((1/2)*(27*z^3-27*z^2+4)^(1/2)-(3/2)*z*(3*z-3)^(1/2))^(1/3))/(z*(3*z-3)^(1/2))

(1)

b := `assuming`([simplify(a, symbolic)], [z >= 1])

2*(-(12*(3*z+1)^(1/2)*z-12*z*(3*z-3)^(1/2)-8*(3*z+1)^(1/2))^(1/3)+(12*(3*z+1)^(1/2)*z+12*z*(3*z-3)^(1/2)-8*(3*z+1)^(1/2))^(1/3))/((3*z-3)^(1/2)*(12*(3*z+1)^(1/2)*z+12*z*(3*z-3)^(1/2)-8*(3*z+1)^(1/2))^(1/3)*(12*(3*z+1)^(1/2)*z-12*z*(3*z-3)^(1/2)-8*(3*z+1)^(1/2))^(1/3)*z)

(2)

plot(1/b, z = 1 .. 10)

 

simplify(diff(1/b, z), symbolic)

-48*(((3*z-2)*(3*z+1)^(1/2)+z*(3*z-3)^(1/2))*((12*z-8)*(3*z+1)^(1/2)-12*z*(3*z-3)^(1/2))^(1/3)+((12*z-8)*(3*z+1)^(1/2)+12*z*(3*z-3)^(1/2))^(1/3)*((-3*z+2)*(3*z+1)^(1/2)+z*(3*z-3)^(1/2)))/((3*z+1)^(1/2)*(3*z-3)^(1/2)*((12*z-8)*(3*z+1)^(1/2)+12*z*(3*z-3)^(1/2))^(2/3)*((12*z-8)*(3*z+1)^(1/2)-12*z*(3*z-3)^(1/2))^(2/3)*(((12*z-8)*(3*z+1)^(1/2)-12*z*(3*z-3)^(1/2))^(1/3)-((12*z-8)*(3*z+1)^(1/2)+12*z*(3*z-3)^(1/2))^(1/3))^2)

(3)

``


 

Download simplification.mw

Maple provides efficient vectorization and automatic parallelization for many common operators. For example

x -> 2*~x*~cos~(x*~x)

But in my application it is common to want to create rather long vectorized operators starting from some complicated symbolic computations. Doing conversions by hand from symbolic expressions to element-wise operations is laborious and error prone.

As a very simple example consider that it is possible to obtain (almost) the same result as above by writing the following as a vectorized operation

D(x->sin(x^2))~

But there are at least two problems with this. First of all it is not nearly as efficient as the first operator and second, perhaps not unrelated, is that the datatype returned when applying this operator to a Vector/rtable of hardware floats (e.g. datatype=float[8]) becomes something  more general.

My question is how can I convert a complicated symbolic expression into an efficient numeric element-wise vector operation?

I have tried several different approaches but so far without success. In the case above for example it seemed natural to expect that the following derivative

D(x->sin~(x^~2))

would produce a vectorized result, but this is not the case. In another attempt I was unable to see how to perform substitions into an expression, e.g. like this

unapply(subs(`*`=`*`~, cos=cos~, diff(sin(x),x)), x)

I would be glad to receive suggestions and/or references to relevant documentation. 

 

Everything is simple, until you go underwater – This is what the University of Waterloo Submarine Racing team, or in short ‘WatSub’ coined as their motto. Never mind learning to scuba dive, and dealing with such things as rust, this newly formed team would have to compete against university teams with a decade or more of experience.

But that did not deter the team, and they started work on Ontario’s first submarine racing project. The team approached Maplesoft to be a sponsor and we are proud to have supported this ingenious venture. The team has used Maplesoft technology in the design and testing of the submarine.

“Maple has been our go-to calculations and analysis tool throughout the development of Amy (2015-2016 season), and we will continue using it throughout the development of Bolt (2016-2017 season),” said Gonzalo Espinoza Graham, President of the WatSub Team. “Its familiar interface and computing environment allowed us to set design benchmark targets from early on the design process and follow through with them on the later stage.”

What started as an engineering project in December 2014, becoming officially the first submarine racing team in Ontario. The team soon grew to over 130 general members and a tight core-team, who were eager to tackle new challenges.  The team resides inside the Sedra Student Design Centre, University of Waterloo’s state of the art facility that houses over 25 student teams, the largest of its kind in North America.  

WatSub made its first appearance on the European International Submarine Races (eISR) back in July 2016, with its 1st submarine ‘Amy’, where a single scuba diver piloted the submarine and propelled it through an unforgiving winding course marked by obstacles and turns 10 meters underwater. The team has since then participated in other competitions and is constantly improving the design and performance of the submarine, learning from each competition they participate in.  Next year Amy will participate in the 14th edition of the eISR international competition. “I think the greatest thing we learned is never to give up,” said Ana Krstanovic, a third-year political science student who manages communications for the team. “We’re more motivated now than ever.”

 

Ojaswi Tagore, Gonzalo Espinoza Graham, and Janna Henzl represented WatSub at the European International Submarine Race in Gosport, UK.

 

Another example of an innovative project that Maplesoft supported in 2016 is Waterloop: The Canadian SpaceX Hyperloop Competition Team, Canada's only SpaceX Hyperloop Pod Competition team. This project, which could change the way we travel in the future, is driven by a group of dedicated University of Waterloo students who have taken on the challenge to design and build a functional prototype Hyperloop pod. They will test it on a one-mile test track in Hawthorne, California in January 2017, pitting it against 22 of the 1200+ teams who originally entered the competition.

The Hyperloop is a conceptual next generation high-speed transit system that will take commuters between cities at speeds over 1,000 km/h. The technology will differ from previous rail transit by having pods ride on a cushion of air in a reduced pressure tube in order to reach greater speeds with a smoother ride, and is powered entirely by renewable energy.

 The Hyperloop Pod Competition was launched by Elon Musk, the billionaire engineer and founder of SpaceX and Tesla Motors.  The competition is separated into 3 rounds. The first one was held in late December, where selected teams sent in their initial designs to be reviewed. From there, 180 teams were chosen to compete at Texas A&M University. Each team set up a booth and a panel of judges critiqued them and chose 31 teams to move onto the final, build and test stage.

Waterloop Goose I

Waterloop Goose X

The GOOSE I is Waterloop’s half-scale, functional prototype vehicle pod, which will be the one in the competition.  The GOOSE X pod is a conceptual full size Hyperloop vehicle inspired by the prototype they are building. The full size pod will have a capacity of 26 passengers per pod.

"Our prototype has been designed to be as simple and economical as possible, while still performing all necessary functions for the full size Hyperloop. If it is successful, it has the potential to revolutionize the transit industry in the same manner the train and airplane has before it," said Montgomery de Luna, architectural design lead for Waterloop. “We would like to thank Maplesoft for their generous support.  Without sponsors like Maplesoft supporting our vision and encouraging innovative student projects, we wouldn’t be able to achieve our goal.”

Revolutionizing the transportation industry isn’t easy and is at times frustrating and time consuming for these teams, but having the best tools and resources will ensure that the teams have a good chance at excelling in competitions and creating innovative models that could change our future.

I use the example procedure when search. Clock in help 

but elapsed function can only run one time

because it return clock is not running

need to run clock start again and calculate from beginning again

how elapsed function can run more times

Hi!

Everyone,

I want to draw  phase plane of system of three fractional order equations. 

 

Note that 

Also want the  phase portrait when the values of alpha are not same....

Also

Thanks

 

 

 

I want to get solutions of this system ,can anyone help me ?solutions.mw

The Joint Mathematics Meetings are taking place this week (January 4 – 7) in Atlanta, Georgia, U.S.A. This will be the 100th annual winter meeting of the Mathematical Association of America (MAA) and the 123nd annual meeting of the American Mathematical Society (AMS).

Maplesoft will be exhibiting at booth #118 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.

There are also several interesting Maple-related talks and events happening this week:

 

Teaching Cryptology to Increase Interest in Mathematics for Students Majoring in Non-Technical Disciplines and High School Students

Wednesday, January 4, 0820, L401 & L402, Lobby Level, Marriott Marquis

Neil Sigmon, Radford University

 

Enigma: A Combinatorial Analysis and Maple Simulator

Wednesday, January 4, 0900, L401 & L402, Lobby Level, Marriott Marquis

Rick Klima, Appalachian State University

 

MYMathApps Calculus - Building on Maplets for Calculus

Thursday, January 5, 0800, Courtland, Conference Level, Hyatt Regency

Philip B. Yasskin, Texas A&M University 
Douglas B. Meade, University of South Carolina 
Andrew Crenwelge, Texas A&M University

 

Maple Software Technology as a Stimulant Tool for Dynamic Interactive Calculus Teaching and Learning

Thursday, January 5, 1000, Courtland, Conference Level, Hyatt Regency

Lina Wu, Borough of Manhattan Community College-The City University of New York 

 

Collaborative Research: Maplets for Calculus

Thursday, January 5, 1400, Marquis Ballroom, Marquis Level, Marriott Marquis

Philip Yasskin, Texas A&M University 
Douglas Meade, U of South Carolina

 

Digital Graphic Calculus Art Design in Maple Software

Thursday, January 5, 1420, International 7, International Level, Marriott Marquis

Lina Wu, Borough of Manhattan Community College-The City University of New York 

 

Maplesoft will also be hosting a catered reception and brief presentation on Teaching STEM Online: Challenges and Solutions, Thursday January 5th, from 6:00pm – 7:30pm, at the Hyatt Regency, Hanover AB, on the exhibitor level. Please RSVP at www.maplesoft.com/jmm or at Maplesoft booth #118.

 

If you are attending the Joint Math meetings this week and plan on presenting anything on Maple, please feel free to let me know and I'll update this list accordingly.


See you in Atlanta!

Daniel

Maple Product Manager

(1)How can i draw a graph by considering eta on x-axis and f'(eta),theta(eta)on y-axis in a single graph with respect to variatiation  in the parameter beta=0.01,0.1,1.0.. 

(2). how can i get different values of f'(eta) by varying values of eta .

restart; with(plots); beta := 0.1e-1; Bi := 10; Pr := 3.0; L0 := 1; w := 0.2e-1

Eq1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))^2+beta*H(eta)*(F(eta)-(diff(f(eta), eta))) = 0

diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))-(diff(f(eta), eta))^2+0.1e-1*H(eta)*(F(eta)-(diff(f(eta), eta))) = 0

(1)

Eq2 := G(eta)*(diff(F(eta), eta))+F(eta)^2+beta*(F(eta)-(diff(f(eta), eta))) = 0

G(eta)*(diff(F(eta), eta))+F(eta)^2+0.1e-1*F(eta)-0.1e-1*(diff(f(eta), eta)) = 0

(2)

Eq3 := G(eta)*(diff(G(eta), eta))+beta*(f(eta)+G(eta)) = 0

G(eta)*(diff(G(eta), eta))+0.1e-1*f(eta)+0.1e-1*G(eta) = 0

(3)

Eq4 := H(eta)*F(eta)+H(eta)*(diff(G(eta), eta))+G(eta)*(diff(H(eta), eta)) = 0

H(eta)*F(eta)+H(eta)*(diff(G(eta), eta))+G(eta)*(diff(H(eta), eta)) = 0

(4)

Eq5 := (diff(theta(eta), eta, eta))/Pr+f(eta)*(diff(theta(eta), eta))+(2*H(eta)*beta*(1/3))*(thetap(eta)-theta(eta)) = 0

.3333333333*(diff(diff(theta(eta), eta), eta))+f(eta)*(diff(theta(eta), eta))+0.6666666667e-2*H(eta)*(thetap(eta)-theta(eta)) = 0

(5)

Eq6 := G(eta)*(diff(thetap(eta), eta))+L0*beta*(thetap(eta)-theta(eta)) = 0

G(eta)*(diff(thetap(eta), eta))+0.1e-1*thetap(eta)-0.1e-1*theta(eta) = 0

(6)

bcs1 := f(0) = 0, (D(f))(0) = 1, (D(theta))(0) = -Bi*(1-theta(0)), (D(f))(5) = 0, F(5) = 0, G(5) = -f(5), H(5) = w, theta(5) = 0, thetap(5) = 0;

f(0) = 0, (D(f))(0) = 1, (D(theta))(0) = -10+10*theta(0), (D(f))(5) = 0, F(5) = 0, G(5) = -f(5), H(5) = 0.2e-1, theta(5) = 0, thetap(5) = 0

(7)

p := dsolve({Eq1, Eq2, Eq3, Eq4, Eq5, Eq6, bcs1}, numeric);

proc (x_bvp) local res, data, solnproc, _ndsol, outpoint, i; option `Copyright (c) 2000 by Waterloo Maple Inc. All rights reserved.`; _EnvDSNumericSaveDigits := Digits; Digits := 14; if _EnvInFsolve = true then outpoint := evalf[_EnvDSNumericSaveDigits](x_bvp) else outpoint := evalf(x_bvp) end if; data := eval(`dsolve/numeric/data/modules`[1]); solnproc := data:-Get("soln_procedure"); if not type(outpoint, 'numeric') then if outpoint = "solnprocedure" then return eval(solnproc) elif member(outpoint, ["start", "left", "right", "errorproc", "rawdata", "order", "error"]) then return solnproc(x_bvp) elif outpoint = "sysvars" then return data:-Get("sysvars") elif procname <> unknown then return ('procname')(x_bvp) else _ndsol := pointto(data:-Get("soln_procedures")[0]); return ('_ndsol')(x_bvp) end if end if; try res := solnproc(outpoint); [eta = res[1], seq('[F(eta), G(eta), H(eta), f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), theta(eta), diff(theta(eta), eta), thetap(eta)]'[i] = res[i+1], i = 1 .. 9)] catch: error  end try end proc

(8)

odeplot(p, [eta, f(eta)], 0 .. 5);

 

odeplot(p, [eta, diff(f(eta), eta)], 0 .. 5);

 

NULL

odeplot(p, [eta, thetap(eta)], 0 .. 5);

 

odeplot(p, [[eta, F(eta)], [eta, thetap(eta)]], 0 .. 5);

 

``


 

Download from_net_(1).mw

 

I am trying to plot f=[r^2 *cos(theta)+r*sin(theta)],as a DensityPlot[] in polar coordinates.

here r is a function of theta, r=r1(θ)..r2(θ) , for this reason, I have a problem to plot f in density plot

f=[r^2 *cos(theta)+r*sin(theta)];

r1(theta)= 0.3+0.1*cos(theta);
r2(theta)= 0.5+0.1*cos(theta);
display(changecoords(densityplot(f, r =r1(theta)..r2(theta), theta = 0 .. 2*Pi, style = patchnogrid, colorstyle = HUE), polar), axes = box, orientation = [270, 0], labels = [x, y, ``]);

(Error, (in plots/densityplot) bad range arguments r = .3+.1*cos(theta) .. .5+.1*cos(theta), theta = 0 .. 2*Pi )

The ODE diff(y(x),x) = sec(x)^2*sec(y(x))^3  can be solved as separable. So the answer should be 

simple_answer:=sin(y(x))*(cos(y(x))^2+2)=C_1+3*tan(x);

as can be seen by direct integration of each side of the differential equation. I am trying to make Maple give the same answer, but not having any luck. 

restart;
ode:=diff(y(x),x) = sec(x)^2*sec(y(x))^3;
sol:=dsolve(ode, y(x),implicit);

I tried simplify(sol,trig) and tried simplify(sol,size).  Both Maple answer, and the simple answer solve the ODE.

Is there a way to make Maple dsolve give the simpler answer, or simplify/convert the answer it gives to the simpler one? I am newbie in Maple.

I am using Maple 2015 and when assigning a name to a function or expression and hit enter, it appears like in the following picture.

However, I want to make it appear like in the following picture:

How can I make that ?

Thank you !

Can someone help me with this:

reverse_eng.mw

happy new year!

I was trying to solve a system of ODE using Maple, but to my surprise, Maple recognizes diff((phi(t), t)) as a variable which is different than t. 

My code is as following:

dsys := {2*m1*(a+l*sin(phi(t)))^2*(diff(diff(theta(t), t), t))+4*m1*(a+l*sin(phi(t)))*l*cos(phi(t))*(diff(theta(t), t))*(diff(phi(t), t)) = M, 2*m1*l^2*(diff(diff(phi(t), t), t))+4*m2*l^2*sin(2*phi(t))*(diff(phi(t), t))*(diff(phi(t), t))+4*m2*l^2*sin(phi(t))^2*(diff(diff(phi(t), t), t))-2*m1*(a+l*sin(phi(t)))*l*cos(phi(t))*(diff(theta(t), t))*(diff(theta(t), t))-2*m2*l^2*(sin(2*phi(t)))(diff(phi(t), t))*(diff(phi(t), t)) = -(2*(m1+m2))*g*l*sin(phi(t))-2*k*l^2*sin(2*phi(t)), phi(0) = 0, theta(0) = 0, (D(phi))(0) = 0, (D(theta))(0) = 0}

subs({M = 10, a = .5, g = 9.81, k = .1, l = .5, m1 = 10, m2 = 1}, dsys);

dsn1 := dsolve(dsys, numeric)

The error I got was Error, (in dsolve/numeric/process_input) input system must be an ODE system, got independent variables {t, diff(phi(t), t)}

I don't get why this is happening. Could you show me what's going on?

   

 

The code for the animation:

L:=[[-0.12,2],[-0.14,0],[0.14,0],[0.12,2]]:
L1:=[[0.05,2],[4,1],[2,4],[3.5,3.5],[1,7],[2,6.5],[0,10]]:
A:=plot(L, color=brown, thickness=10):
B:=plot([op(L1),op(map(t->[-t[1],t[2]],ListTools:-Reverse(L1)))], color="Green", thickness=10):
C:=plottools:-polygon([op(L1),op(map(t->[-t[1],t[2]],ListTools:-Reverse(L1)))], color=green):
Tree:=plots:-display([A, B, C], scaling=constrained, axes=none):
T:=[[-3.2,-2, Happy, color=blue, font=[times,bold,30]], [0,-2,New, color=blue, font=[times,bold,30]], [2.5,-2,Year, color=blue, font=[times,bold,30]], [-5,-3.5, "&", color=yellow, font=[times,bold,30]],[-2.5,-3.5, Merry, color=red, font=[times,bold,30]], [2.3,-3.5, Christmas!, color=red, font=[times,bold,30]], [0,-5, "2017", color=cyan, font=[times,bold,36]]$5]:
F:=k->plottools:-homothety(Tree, k, [0,5]):
A:=plots:-animate(plots:-display, ['F'(k)], k=0..1, frames=60, paraminfo=false):
B:=plots:-animate(plots:-textplot,[T[1..round(i)]], i=0..nops(T), frames=60, paraminfo=false):
plots:-display(A, B, size=[500,550], scaling=constrained);


Christmas_Tree.mw

 Edit.

 

First 1000 1001 1002 1003 1004 1005 1006 Last Page 1002 of 2215