MaplePrimes Questions

Hi

Sometimes when I'm solving an equation I get the result in the form of a constant multiplied by zero, instead of just the number 0.

Is it possible to remove that option so I won't get a confusing answer? I've attached a picture of an example here

Thanks!

2.mw
 

416*(1-tanh(rho*y)^2)^2*tanh(rho*y)^3-272*(1-tanh(rho*y)^2)^3*tanh(rho*y)-32*tanh(rho*y)^5*(1-tanh(rho*y)^2)

416*(1-tanh(rho*y)^2)^2*tanh(rho*y)^3-272*(1-tanh(rho*y)^2)^3*tanh(rho*y)-32*tanh(rho*y)^5*(1-tanh(rho*y)^2)

(1)

``


 

Download 2.mw
 

416*(1-tanh(rho*y)^2)^2*tanh(rho*y)^3-272*(1-tanh(rho*y)^2)^3*tanh(rho*y)-32*tanh(rho*y)^5*(1-tanh(rho*y)^2)

416*(1-tanh(rho*y)^2)^2*tanh(rho*y)^3-272*(1-tanh(rho*y)^2)^3*tanh(rho*y)-32*tanh(rho*y)^5*(1-tanh(rho*y)^2)

(1)

``


 

Download 2.mw
 

416*(1-tanh(rho*y)^2)^2*tanh(rho*y)^3-272*(1-tanh(rho*y)^2)^3*tanh(rho*y)-32*tanh(rho*y)^5*(1-tanh(rho*y)^2)

416*(1-tanh(rho*y)^2)^2*tanh(rho*y)^3-272*(1-tanh(rho*y)^2)^3*tanh(rho*y)-32*tanh(rho*y)^5*(1-tanh(rho*y)^2)

(1)

``

 How can I collect these  for  simplest form 

 

 

 

Hi,

I would like to ask a very basic question. The solution would indeed simplify my life!

I have a function called (eu) :

eu := -2*n^2*B^2+2*n*B^2*p-2*n*B^2*p^2+n*B^2-B^2*p+B^2*p^3+n*c*d

I would like to collect it with respect to two specific terms: (n-p+p^2) and (n-p+p^3), that would give me at the end:

eu := -2*n*B^2*(n-p+p^2)+B^2*(n-p+p^3)+n*c*d

Collect command gives me the following error meesage: "Error, (in collect) cannot collect n-p+p^2".

Would you have a solution for this problem?

Thanks in advance.

Basak

The program below is meant to use the  solve command to find the sides a, b, c of a triangle, given the lengths of the three altitudes, hA, hB and hC.  This program came from a higher version of Maple and I wondered what the word 'explicit' is meant to do in the  Sol->[solve({hA=(...statement.  The output just gives the name Sol - when I would have liked it to either give a solution for the sides in terms of the altitudes hA,hB,hC - or numeric float values.  For the attempted values of the altitudes all being equal to sqrt(30/2, the triangle is equilateral of side 1.

  Any suggestions to get this to work in Maple 7would be most appreciated.

Thanks,

   David

 

 

> restart:

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

# Altitude exploration 

# Trying to obtain expressions for the sides a,b,c of a triangle
# given the altitudes.

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

#with(plots):

#with(plottools):

#with(geometry):

Sol->[solve({hA=(a+b+c)/a, hB=(a+b+c)/b, hC=(a+b+c)/c},{a,b,c}  )];  #,explicit)];

#  ???Is explicit needed?? Is it a Maple 7 feature?

 

abc:=simplify(eval(Sol,[hA=sqrt(3)/2,hB=sqrt(3)/2, hC=sqrt(3)/2]));  #should give a=b=c=1

gf:=eval(Sol, a=1);

 

Warning, the name changecoords has been redefined

 

I've always had problems with math and numbers.  My goal is to learn enough by myself so that I can take a college algebra class.  Do you have any products that I could buy as an individual which would assess my math skills and design a program to bring me up to speed so I can start with college math classes?

Any other advice on learning math as an adult is welcome.

Thanks

Hi Dears,

I know "randpoly" command of Maple and I need to generate some random polynomials with parametric coefficients. For this I used this command for input set {a,b,c,d,x,y,z,w} which {a,b,c,d} are parameters and {x,y,z,w} are variables. But, the generated polynomials were not interestting. For example, some polynomials such as $(a-b)xy+(c^2-a)x*y*z-a*b+(d-1)z^2w$ are interested for me. In fact, difference parametric coefficients is important. 

Could you please help me to generate them?

Sincerely yours 

I am interested in the 5 circle theorem of Miquel.  Search on the internet 'Miquel five circle theorem' for more details. I would like to prove this theorem using Maple, and also see if there is a generalisation to this for more than 5 circles.  I wish to find the points of intersection of the circles and am using the fsolve command:

fsolve({(x-x[i])^2+(y-y[i])^2-r[i]^2, (x-x[i-1])^2+(y-y[i-1])^2-r[i-1]^2}, {x, y});

I am using the curly braces for sets - as I can't seem to get it to work for [] lists.   The output gives something like {x=12.0005, y=4.65}.  I want to use these values to obtain straight line equations and verify that the lines formed by successive circles form a pentagram, with all vertices on the five circles.  I just want to get to the floating point values, without the x= part.  The type of the returned expressions is '=' - whatever that means!

    I'm also wondering if using the plottools and plots packages is sufficient - as opposed to the geometry pakage.

I'm interested in how many people have heard of this theorem.  Does it have any generalisations to 6, 7, ... circles?

Any help or comments gratefully received.

    David.   

Hellow dears!!!

Hope everyone is fine with everything. I want the animation of the numerical solution of ODE i,e., f4 for delta=[0,1.5]. Please see the attachment and fix my problem. I Shall be very thankful to you.

Graph.mw

Special request to  acer 13834@Carl Love ,Preben Alsholm 10271

As part of the class I am teaching in biofluid mechanics, the students are learning about the Buckingham Pi theory for finding dimensionless groups (pi-groups).  The process involves the dimensions of parameters in terms of their basic units of measure: mass, M, length, L, and time, T.  Fore example, fluid density is M/L3 .  A typical step is to represent a product of parameters in terms of these measures raised to exponents.  Here is an example

(M/LT)a(M/L3)b(L2)c

where the exponents a,b,c are to be determined. What command(s) will massage this form to look like

Ma+b  L2c-3b-a T-a  ?

Cheers

 

Greetings to all.

A basic Euler-MacLaurin Summation recently appeared at this Math.Stackexchange Link . Several users computed an asymptotic expansion for the sum given by sum(1/k^sigma, k=1..n) with sigma a parameter between zero and one. E.g. for sigma=1/3 we find

       2/3
    3 n                     1           1           7            13            247
    ------ + Zeta(1/3) + -------- - --------- + ---------- - ----------- + ------------
      2                     (1/3)       (4/3)         10/3          16/3           22/3
                         2 n        36 n        4860 n       26244 n       590490 n

Comparing this with Euler-Maclaurin Maple produces the following:

> eulermac(1/n^(1/3), n, 8);
       2/3
    3 n         1           1           7            13            247             1
    ------ - -------- - --------- + ---------- - ----------- + ------------ + O(-------)
      2         (1/3)       (4/3)         10/3          16/3           22/3      (28/3)
             2 n        36 n        4860 n       26244 n       590490 n         n

The sign on the second term is incorrect. I would be curious to know what is happening here. BTW Mathematica also yields a plus sign on this term, same as what was obtained at the link.

Thank you for any comments you may have.

Marko Riedel

with(plots):
xy := MatrixMatrixMultiply(Matrix([[cos(t),sin(t)],[sin(t),cos(t)]]),Matrix([[x],[y]]));
animate3d([t, xy[1][1], xy[2][1]], x = 1 .. 10, y = 1 .. 10, t = -Pi .. Pi,coords=spherical);
animate3d([x, xy[1][1], xy[2][1]], x = 1 .. 10, y = 1 .. 10, t = -Pi .. Pi,coords=spherical);
animate3d([t, xy[1][1], xy[2][1]], x = -10 .. 10, y = -10 .. 10, t = -Pi .. Pi,coords=spherical);
 
would like to see how it rotate a ball
 

Hi 

I need yours help to solve the given integral 
 

restart

A = 3*(q-1)^(-1/(1+r))*GAMMA(q-3/(2+2*r))*GAMMA(3/(2+2*r))/(GAMMA(5/(2+2*r))*GAMMA(q-5/(2+2*r)));

g := int(`ε`*(1+(`ε`/(A*t))^(1+r)/(q-1))^(-q)*c*exp(-`ε`/a), `ε` = 0 .. infinity);

``

NULL

NULL

q,c,a and r are parameters.
Γ is a gamma function.

Download Karim.mw
 

 

 

I am trying to draw a tetrahedron and square based pyramid, having its triangular faces the same asthat of the tetrahedron.  While the geom3d package in Maple 7 has several polyhedra, a pyramid, or pentahedron(?) seem to be missing.  The closest is the octohedron - but this is two square pyramids "glued" together!   I'm totally bamboozled by the Polyhedron command with all its options.  Is it possible to somehow cut the octohedron in two, to get a pyramid?  Below is my Maple program, including at the end aplottools command to draw a separate ocotohedron.  

  Any help gratefully received.

     David

restart:
> #Puzzle Problem with tetrahedron and pyramid.  How many faces will
> #resultant polyhedron have when a triangular face of the pyramid is "glued" to  the tetrahedron.
> with(geom3d):

> print(`A pyramid and tetrahedron are shown: the triangular faces of each being the same dimensions. The`);
> print(`triangular face of the pyramid is attached to a triangular face of the tetrahedron.`);

> print(`How many faces does the resulting polyhedron have?`);
> print(`  a.) 5    b.) 6     c.) 7     d.) 8     e.) 9`);

> # octahedron, hexahedron, cube, icosahedron, dodecahedron - but pentahedron not supported in Maple 7 :-(
> RegularPolyhedron(d,[3,4],point(o,4,4,0),3):
> #dodecahedron(t,point(p,5,0,0),3):

> tetrahedron(t,point(p,5,0,0),3):
> #draw([d(color=red),t(color=green)],cutout=7/8,lightmodel=light4,
> draw([d(color=red),t(color=green)],title=`Tetrahedron & pyramid`,orientation=[45,45]);
> with(plottools):
> f := octahedron([0,0,0],1):  #, octahedron([1,1,1],0.5):
> plots[display](f,style=patch);

Dear Team, I have used RKF45 to solve my ODE with Maple. now I am required to solve same ODE using RK4 for comparison of solution. Ps help me with an example. Pls find my Parameters, intial values and ODE below:


 

``

ODE*equations

ODE*equations

(1)

diff(s(t), t) = (1-phi)*epsilon+(1-rho)*a+(1-f)*alpha*v(t)-(lambda+theta[1]+a+epsilon)*s(t)

diff(v(t), t) = phi*epsilon+rho*a+theta[1]*s(t)-((1-f)*alpha+f*theta[2]+a+epsilon)*v(t)

diff(e(t), t) = lambda*s(t)-(delta+a+epsilon)*e(t)

NULL

diff(r(t), t) = eta*i(t)+v(t)*f*theta[2]-(a+epsilon)*r(t)

``

``

My*parameters

My*parameters

(2)

v(0) := .4

.4

(3)

NULL

s(0) := 0.6e-1

0.6e-1

(4)

e[0] := .24

.24

(5)

i[0] := .17

.17

(6)

r[0] := .13

.13

(7)

c := 0.4e-1

0.4e-1

(8)

f := .4

.4

(9)

beta := .2

.2

(10)

epsilon := .8

.8

(11)

theta[1] := .1

.1

(12)

theta[2] := .3

.3

(13)

alpha := .9

.9

(14)

rho := .7

.7

(15)

eta := .99

.99

(16)

delta := .3

.3

(17)

a := 0.4e-1

0.4e-1

(18)

phi := 1

``

``

``


 

Download ODE_EQNS.mw

How can I solve this non-linear ODE exactly or approximately in some series solution?

Non_linear_ode.mw

First 886 887 888 889 890 891 892 Last Page 888 of 2433