How I can plot torus structure in the following code instead of cylindrical.

Thanks.

 

>  >    >    >  (1) >  >  (2) >

 

Download toro.mw

 

 

\Hello,

How I can solve this algebraic to find unknowns ABCD?

I want to gain ABCD automatically without input the coefficients in rule by hand.

Because I should run the code for many input data

Thanks

 

>  restart; >  l:=0.5;a:=0.1; rho:=2700;h:=.0005; E:=72.4*10^9;v:= 0.3; n:=6; m:=1; AD:=10; mu:=(2*a*2.35)/l; nu:=sin(mu*l/(2*a))/sinh(mu*l/(2*a)); omega[m,n]:= 3067.173621; (1) >    >  E:=1:k[1,1]:=-5.660173062*10^10:k[1,2]:=-2.8552873062*10^10:k[1,3]:=-8.68528173062*10^10:k[1,4]:=-7.6788528173062*10^10:k[1,5]:=-1.52568528173062*10^10:k[2,1]:=-15.660173062*10^10:k[2,2]:=-21.8552873062*10^10:k[2,3]:=-18.68528173062*10^10:k[2,4]:=-71.6788528173062*10^10:k[2,5]:=-10.52568528173062*10^10: k[3,1]:=-5.65257260173062*10^10:k[3,2]:=-27.8552552873062*10^10:k[3,3]:=-81.6854428173062*10^10:k[3,4]:=-9.67858528173062*10^10:k[3,5]:=-3.52568528173062*10^10: k[4,1]:=-51.111660173062*10^10:k[4,2]:=-21.811552873062*10^10:k[4,3]:=-18.68528173062*10^10:k[4,4]:=-17.6788528173062*10^10:k[4,5]:=-11.52568528173062*10^10: k[5,1]:=-6.660173062*10^10:k[5,2]:=-61.852873062*10^10:k[5,3]:=-82.68528173062*10^10:k[5,4]:=-72.6788528173062*10^10:k[5,5]:=-21.52568528173062*10^10 (2)   >    >  S:=(Matrix([[rho*h*omega[m,n]^2+k[1, 1],k[1,2],k[1,3],k[1,4]],[k[2,1],rho*h*omega[m,n]^2+k[2,2],k[2,3],k[2,4]],[k[3,1],k[3,2],k[3,3]+rho*h*omega[m,n]^2,k[3,4]],[k[4,1],k[4,2],k[4,3],k[4,4]+rho*h*omega[m,n]^2]])).(Vector(1..4,[[A],[B],[C],[D]]))=-E*(Vector(1..4,[k[1,5],k[2, 5],k[3,5],k[4,5]])); (3) >

 

Download solve.mw

 

 

 

 

Setup   >  >  >  >  >  (1.1) >  (1.2) >  >  (1.3) >  (1.4) >  >  >  >  >  >  >  Boyle's law   >  (2.1) >  (2.2) >  (2.3) >  (2.4) >  (2.5) Conservation of Energy   db is much smaller then da, herfore va is neglected. >  (3.1) >  (3.2) >  (3.3) >  Conservation of Mass   >  (4.1) >  >  Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)   (4.2) >  (4.3) >  >  >  >

Hi i am trying to get a graph that shows height against time but it does not plot and I do not understand why. I copied the code of a friend and it does work on his computer. 

Download laatste_versie_berekeningen.mw

The worksheet below animates the flattening of a tetrahedron by expanding one of its faces, namely its triangular base.

I would like to animate the flattening of an octahedron so that it assumes the 2D figure resembling the Morley triangle which is included in the worksheet.

Are there documents on the web explaining the technique for doing so? Is there a Maple worksheet available on the web demonstrating the desired animation?

Flatten_a_tetrahedron.mw

n is a Carmichael number iff for every prime factor p of n, p-1/n-1.

Question: How to find odd squarefree composite numbers n having k distinct prime divisors, and the property that exactly k-1 prime divisors satisfy the Carmichael requirement, p-1/n-1 ?

Examples: 231,1045,1635. In these cases k=3 and the prime divisors satisfying the criteria are the greatest and smallest. I have a code for this but would like to compute the general case, where the criteria is satisfied for precisely any k-1 divisors.

Any assistance greatly appreciated.

David. 

 

Hi,

I have some trouble using the GlobalOprimization Toolbox (GoT)
When used with some options GoT returns an error message indicating the option is not valid.
For instance options (just a few among some others) "numexperiments", "nugget", "optsearch", "goal", ... are not recognized.

I'm using Maple 2018.0 (I don't know how to recover the number of the GlobalOptimization Toolbox).

Is it an installation problem (GotT has been bought legally) or a problem with this Toolbox ?

TIA

Is there any maple function to test whether infinite sums and integrals can be swapped say for example whether

this is true? for a given f(x):

Hi, my codes ran smoothly well until I changed some assignment value. It reported invalid object error.

Error, Got internal error in Typesetting:-Parse:-Postprocess : "internal error: invalid object "

Part of the codes are as follows:

""sigma:=0.00002:    m:=5:   n:=4:
Error, Got internal error in Typesetting:-Parse:-Postprocess : "internal error: invalid object "

Could anyone give some hints on what goes wrong here? It seems to be a configuration issue but I completely have no clue..

Thanks a million in advance,

Best,

Jie

 

 

In some previous version of Maple there was a function to convert an approximate numerical value, say 1.77245,

to the exact value (if it exists and is short), say sqrt(Pi). I cannot remember the name of this function. Is it still in Maple?

Thanks for any help

Pavel

restart; interface(rtablesize = 10): _EnvHorizontalName := 'x': _EnvVerticalName := 'y': eqPA := (y-b0)/(x-a0) = k: solPA := y=solve(eqPA, y): #k coefficient directeur de PA eqPB := (y-b0)/(x-a0) = -1/k: solPB :=y= solve(eqPB, y):#PB perpendicalaire à PA xA := solve(subs(y = 0, eqPA), x): yB := solve(subs(x = 0, eqPB), y): eqAB := x/xA+y/yB = 1; x k y k eqAB := --------- + --------- = 1 a0 k - b0 b0 k + a0 t := solve(a*(xM+(1/2)*t*a)+b*(yM+(1/2)*t*b)+c = 0, t); 2 (a xM + b yM + c) t := - ------------------- 2 2 a + b #Recherche des coordonnées de la projection d'un point sur une droite D #M(x,y)un point quelconque du plan, M'(x',y') son symé trique dans la symétrie orthogonale d'axe D #le vecteur MM' est colinéaire du vecteur normal n de D; vec(MM')=t.vec(n), n=

Hi! I can´t see what i´m doing wrong! :( Thank you for your help!!!

Can someone help me project ODE solution curves of an ODE onto velocity manifolds (u(x,t), v(x,t))?

Please see MAPLE worksheet below.

Thanks for any help..

Melvin
 

Help needed to graphically map (x(t), t) trajectories onto velocity manifolds (u, v) MRB: 16/6/19 The velocity manifolds u(x,t), v(x,t) are depicted in the figure below: >  The above velocity plot above shows classical pairs of negative (red) low and high velocity, and pairs of positive (blue) low and high velocity.  At any given instant there are fast and slow propagating solutions, from which classical trajectories on (x,t) are derived below for each of the red and blue velocity pulse solutions.  The cut away regions of the velocity manifolds represent complex  and We now solve the velocity ODE to plot a (x,t) trajectory on the red velocity manifold: >  >  (1) If the set of dependent variables is specified as a list, those variables appear in the same order in the output.  Now extract the solution points at given times >  (2) Here is a solution at a specified time: >  (3) Here is the the solution for the (x,t) trajectories: >  (4) We now plot trajectories on (x,t) space: >  >  We do not generate trajectories in the region  for which regions of the manifolds are no longer real-valued.  The above trajectories are solutions of the above coupled Burger's equations for (u(x,t), v(x,t)). We wish plot these trajectories on their  [t,x(t),D(x(t)] manifolds. >  The following is an example of how to project curves onto a surface, assuming an analytic form for the curves.  However, we wish to project the above two ODE solution trajectory curves onto the surfaces u(x,t) and v(x,t) which are solutions of two PDEs.   Below is an example of projecting a specified curve onto a surface  Our requirement is to project a curve which is an integral solution of the ODE (4) above onto one of the velocity manifolds u(x,y), v(x,y).  Below is an example of an algebraically specified trajectory being projected onto a manifold (surface): >  >  (5) >  >  QUESTION: We wish instead to project the ODE  (rather than algebraic) solution curves onto the surface.   Please someone provide MAPLE code to do that? >

 

Download Question_-_projection_of_ODE_solution_curves_onto_manifolds.mw

Hello everyone, here i Denmark we are using comma as a separator for decimal numbers. I have looked in the forum and found that if i changed my region in windows settings, but it is Denmark. When i use number formatting it tells me the output will be with comma as a separator,  but the output comes out with period as the separator. Is there a way to change it manually ?

 

 

I am interested in the root of the equation ((2*n-1)*x + 2*n + 1)*x^n = 1 - x, where 0 < x < 1 and n is a large positive integer.  I believe that the root converges to 1 as n goes to infinity.  How does one obtain an asymptotic estimate of the root for large n?

 

Hello fellows,

I am currently using Maple and surprised I can not solve the following problem. I want to use a function as the upper limit of an integral, for example

f(t) := int(1/x, x = 1 .. t) .

Could anyone tell me why Maple won't understand this?

Thank you for your help, I wish you a nice week.

Mac