MaplePrimes Questions

 

Hi, i want to plot stream tubes for the following problem, can anybody help me? Thanks

 

restart;
with(linalg);
Digits := 15;

velocity=[u(r,z),v(r,z),w(r,z)]

u(r,z)=cr, v(r,z)=r*omega, w(r,z)=0 at z=0

eq1:= 2*f(eta)+diff(h(eta),eta)
eq2 := (1/4)*(diff(h(eta),eta))^2-g(eta)*g(eta)-(1/2)*h(eta)*diff(h(eta),eta,eta)+(1/2)*diff(h(eta),eta,eta,eta);

eq3 := -g(eta)*diff(h(eta),eta)+h(eta)*diff(g(eta),eta)-diff(g(eta),eta,eta);

where f(eta), g(eta) and h(eta) are the solution of above system 

bc:=f(0)=1,f (N)=0,h(0)=0,D(h)(N)=0,D(h)(0)=-2,g(0)=R,g(N)=0; where R is constant 

 

Dear All

Is it possible to change the position of labels so that they do not overlap with curves of plots? Please see the attached file in which I want to shift label for vertical axes to slightly left and label for horizontal axes to the slightly right.

Plot_Label.mw

Why does this code return NULL? Which, earlier than 2017 version of Maple differs? What other Maple function could be used in place of eliminate for the general case?

 

eliminate({x[5] = x[2] / x[1], x[6] = x[3] / x[2], x[7] = x[1] / x[4], x[8] = (2 * x[2] + x[4]) / (2 * x[1] + x[3] + x[4])}, {x[1], x[2], x[3], x[4], x[8]});

 

In piecewise notation (real math symbol), how I could add more cases?

The program below is of two circles in the design of a simple ear-ring: baasically a large circle with a circle of varying radius removed.  I have tried to display the animation with the position of the centreof mass, but the animation is too fast and "jittery".  I'd like to slow the animation down to better see if my calculations are correct.  Is there some way of doing this, other than the Animation/Slower/Faster keys?.

 

restart:

# # # # # # # # # # # # # # # # # # # # # # # # # # # #

# Earing Two circles

# A small circle (sc) of radius r, is drawn inside and tangential to

# a bigger circle,(bc), of radius R.  The origin (0,0) is at the bottom of

# the big circle. The small circle is removed to obatain a shape in the

# form of an earing.  x_bar is the distance from the origin of the

# centroid.

# When the centroid is on the circumference of the small circle then the

# ratio R/r = the golden ratio 1.645(?)

# # # # # # # # # # # # # # # # # # # # # # # # # # # #

with(plots):

with(plottools):

Digits:=4:

R:=20:

#r:=10:

#circles

#sc := disk([0,2*R-r], r, color=white):  #in for loop

bc := disk([0,R], R, color=red):

 

for r from 1 by 1 to R do

  x_bar:=eval((R^2+R*r-r^2)/(R+r)):

printf("x_bar =%4.2f   R/r =%5.3f\n", x_bar, R/r);

  sc[r] := disk([0,2*R-r], r, color=white):

  Pt[r]:=point([0,x_bar],color=black):

#plots[display](sc,bc, scaling =constrained);

#t1[r]:=textplot([4*R/3,4*R/3,cat("xbar",convert(x_bar, string))]):

Rat:=R/r:

t1[r]:=textplot([4*R/3,4*R/3,convert(x_bar, float)]):   t2[r]:=textplot([2*R,4*R/3,convert(Rat, float)] ):

 

  p[r]:= plots[display](Pt[r],sc[r],t1[r],t2[r],bc ):       

end do:
 

#for loop to give individual images

#for r from 1 to R do
 

plots[display]({seq(p[r], r=1..R)}, scaling=constrained,axes=none,insequence=true);


#animate({seq(p[r], r=1..R)},t=1..4, scaling=constrained,axes=none);  #,frames=50);

#plots[display]({p[r]},scaling=constrained,axes=none,insequence=true);

#end do;

The attached worksheet performs two functions:

(1) It lets me print 4 × 6 Index Cards for the short entries in each table.

(2) It allows for easy storage and retrieval of syntax (code).

The worksheet has many tables, each separated by a Page Break.

Questions:

(a) Is there a way to sort all the different tables so they will be arranged in alphebetical order?

(b) When I select one table to print and open the Print Dialog, the "Selection" option is grayed out. (see graphic below).  (1) Is there a way to enable the selection option?  (2) Is there a way to determine what page I am on so I can use the "Pages from...to" option?  If I need to number the pages, will the page numbers reset to parallel a new alphabetic sort order.

Many thanks in advance.  See WC29_4_BY_6_NOTE_CARDS_UNSORTED.mw attached. And see image of Print Dialog below.

Les

 

Hi Dears,

I installed Maple 15 and 18 on my PC. But, when I want to use the PolyhedralSets package the following error is appeared. Is this package only on Maple 2015?

 

With the bests.

 

> solve({a <> 0, b>0});
 solve called with arguments: {a <> 0, 0 < b}, AllSolutions = FAIL, ConditionalSolutions = FAIL, DropMultiplicity = FAIL, Explicit = FAIL, MaxSols = FAIL, SolveOverReals = FAIL, SymbolicSolutions = true, TryHard = FAIL, UseAssumptions = false, domain = default, parameters = {}, parametric = false, split = false
 #(solve,127): _MaxSols := oldmaxsols

Error, (in PiecewiseTools:-Convert) unable to convert
 

 

Dears 

Hope you all are fine. I want to solve the following nonlinear system ODEs numerically

Eq1 := diff(F(eta), eta, eta, eta)-phi1*((diff(F(eta), eta))^2-(diff(F(eta), eta)+diff(G(eta), eta))*(diff(F(eta), eta, eta)))+alpha*(1-phi)^2.5*((diff(F(eta), eta, eta))^2+2*(diff(F(eta), eta))*(diff(F(eta), eta, eta))-(diff(F(eta), eta)+diff(G(eta), eta))*(diff(F(eta), eta, eta, eta, eta))+A*(2*(diff(F(eta), eta, eta, eta))+(1/2)*eta*(diff(F(eta), eta, eta, eta, eta))))-A*phi1*(diff(F(eta), eta)+(1/2)*eta*(diff(F(eta), eta, eta))) = 0; Eq2 := diff(G(eta), eta, eta, eta)-phi1*((diff(G(eta), eta))^2-(diff(F(eta), eta)+diff(G(eta), eta))*(diff(G(eta), eta, eta)))+alpha*(1-phi)^2.5*((diff(G(eta), eta, eta))^2+2*(diff(G(eta), eta))*(diff(G(eta), eta, eta))-(diff(F(eta), eta)+diff(G(eta), eta))*(diff(G(eta), eta, eta, eta, eta))+A*(2*(diff(G(eta), eta, eta, eta))+(1/2)*eta*(diff(G(eta), eta, eta, eta, eta))))-A*phi1*(diff(F(eta), eta)+(1/2)*eta*(diff(F(eta), eta, eta))) = 0;

assoicated with the following BCs

F(0)=0,D(F)(0)=1,G(0)=0,D(G)(0)=p, D(F)(L)=0, D( G(L)=0;

for phi1 := 1.2; alpha := 2; A := 1.5;phi:=0.1;

Please help me to find the solution

Special request to @acer@Carl Love @Preben Alsholm

I need help me, i need find maximum value,this form is false LL := proc (x, y) option inline; abs(wr(x, y)-w3(x, y)) end proc; CodeTools:-Usage(evalhf(maximize(LL(DD*i, DD*j), i = 0 .. N, j = 0 .. N))); Error, (in CodeTools:-Usage) cannot handle unevaluated name `i` in evalhf

Hi there,

I was given the following tasks (Only the last part of Ex 3,but Ex 2 is relevant,further explanation in my maple code:)

questionLAB.mw

Thanks!

EDIT:Just in case if anyone is completely baffled by the lack of context,here is the original version of Lab 1.

Context.mw

Hello,

How can i solve this integro-PDE(partial diffrential equation)??

regards...

eq := (1+6*(l/h)^2/(1+nu))*(diff(u(xi, tau), xi, xi, xi, xi)+int(-B*lambda*exp(-lambda(tau-s))*(diff(u(xi, s), xi, xi, xi, xi)), s = 0 .. tau))+diff(u(xi, tau), tau, tau) = alpha*(int((diff(u(xi, tau), xi))^2, xi = 0 .. 1)+int(-B*lambda*exp(-lambda(tau-s))*(int((diff(u(xi, tau), xi))^2, xi = 0 .. 1)), s = 0 .. tau))*(diff(u(xi, tau), xi, xi))+V^2*(sum(j*u(xi, tau)^(j-1), j = 1 .. 8))

jing-Fu.mw

hai,

 i want to plot a graph here is my codes and i want the graph like this, can any body help me

 

restart:
with(DEtools):
with(plots):

h(z):=1-(delta2/2)*(1 + Cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h^4)-(sin(alpha)/F)-h^2+Nb*h^4):
lambda:=int(K1,z=0..1):
F:=0.3:
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
evalf(lambda);

The program below is a high school problem, related to the area a horse can graze, given it is tethered to a rectangular barn.  The level of difficulty is related to the length of rope.  

   I wanted to display some graphics of the field, barn and tethered horse - this latter being the most difficult.  I experimented with a .png picture of a plain silouhette of a horse, imported this into Photoshop, then saved it as a .pdf file, importing this into Maple,  I managed to import this into the worksheet, but I wanted a scaled down version of the horse in the program plots[display] section.  I was unsuccessful in this.  Undeterred, I decided to try and draw a version of a horse using the plots/plottools packages.  The resulting "horse" looks more like a cat, warthog,  mouse or chameleon! 

   I understand later versions of Maple are able to import graphic images.   I'd appreciate some feedback as to how easy this is, and the quality of the resulting images in Maple output.

Thanks,

    David  .  .    

 

restart:

# # # # # # # # # # # # # # # # # # # # # # # # # # # #

# Horse tethered to barn - what area of grass?

# # # # # # # # # # # # # # # # # # # # # # # # # # # #

with(plots):

with(plottools):

macro(palegreen=COLOR(RGB, .5607, .9372, .5607)):

col1:=`black`:

print(`A horse is tethered to the corner of a barn, 10 m wide and 20 m long.  Find the area`);

print(`the horse can graze, if the length of rope is:`);

print(` i.)   5 m   ii.)  25 m  amd iii.)  50 m`);

#Length of rope

L:=10:

#Dims of barn

len:=20:wid:=10:

#Position of bottom left corner of barn

x_barn:=50:y_barn:=50:

len:=20:wid:=10:

 

#dimensions of field

flen:=120:fwid:=110:

 

rect_barn:=rectangle([x_barn,wid+y_barn], [len+x_barn,y_barn], color=brown):

rect_field:=rectangle([0,fwid], [flen,0], color=palegreen):

 

#Position of horse

x0:=37:y0:=32:

a := 4: b := 2.5:

belli := ellipse([x0,y0], a, b, filled=true, color=col1):

legf:=line([x0-1,y0-2], [x0-2,y0-6], color=col1, linestyle=1, thickness=1):

legf2:=line([x0-2,y0-2], [x0-3,y0-6], color=col1, linestyle=1, thickness=2):

rleg:=line([x0+1,y0-2], [x0+2,y0-6], color=col1, linestyle=1, thickness=2):

 

rleg2:=line([x0+2,y0-2], [x0+3,y0-6], color=col1, linestyle=1, thickness=2):

head := polygon([[x0-6,y0+3],[x0-5,y0+4], [x0-2,y0+5], [x0-3,y0+2]], color=brown, linestyle=3, thickness=2):

tail:=line([x0+6,y0-4], [x0+4,y0+1], color=col1, linestyle=1,thickness=2):

 

a := arc([x0+13,y0+3], 15, Pi/2..Pi+.1, color=blue):

plots[display](a,tail,head,belli,legf,legf2,rleg,rleg2,rect_barn,rect_field, scaling=constrained, axes=none);

for example, we need define  inverse function of g(x) where 

g(x)=int(exp(z^2), z = 0 .. x).

 

 

1 2 3 4 5 6 7 Last Page 1 of 1554