Hello,

I have a question: Hw can I add a straight line in the following plot?

2018-01-20_Biegeversuch_Ablauf.mw

Thank you very much!

Martin

Dear sir,

 i want to plot a graph with  tick marks on all sides , i am attaching codes and  a sample graph, see the attachmnent,

The graph range should be displace only on x and y axis and the vaues should not be on opposite to x and y axis

please see the sample graph

restart:
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4):
lambda:=Int(K1,z=0..1):
F:=0.3:
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
plot( [seq(eval(lambda, Nb=j), j in [0.1,0.2,0.3])], delta2=0.02..0.1);

hi everyone
 

I am trying to call the maple function fsolve in the matlab. I have used maple('fsolve(x^2-4)'). It worked.Now in the matlab code I have defined another variable. let us say a = 10.  Now I am trying to call the fsolve as follows.

a=10;

maple('fsolve(a*x^2-1)')

I ended up with the following error :

Error, (in fsolve) number of equations, 1, does not match number of variables, 

could any body help? tnx

How do you way text below

Maple 2015 

if i use function cat() is not 

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?

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

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