MaplePrimes Questions

if just expect to find whether equation is like

(a-b)*(c-d) + (e-f)*(g-h) = i

if use equal , it will care the variable name, is it possible treat the variable name as a wildcard to match the equation structure only 

the structure i mean is *, +, -, =, ()

I followed the steps given in the example of Maple help "Example 5: Adding a MathMLViewer Region" but I get the next error window.

Do you know what is my mistake?

thank you.

I am experimenting with rotating a group of pionts, lines cylinders etc. the only aminations that  will display are where I put the transformations directly as the point coordinates. This is very tedious. Is there a way around this problem efficiently? Attached worksheet shows problem(s). Ultimately I will have a whole group of objects to rotate.
 

restart

with(plots)

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, shadebetween, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

(1)

with(plottools)

[annulus, arc, arrow, circle, cone, cuboid, curve, cutin, cutout, cylinder, disk, dodecahedron, ellipse, ellipticArc, exportplot, extrude, getdata, hemisphere, hexahedron, homothety, hyperbola, icosahedron, importplot, line, octahedron, parallelepiped, pieslice, point, polygon, prism, project, rectangle, reflect, rotate, scale, sector, semitorus, sphere, stellate, tetrahedron, torus, transform, translate]

(2)

``

``

``

a1 := proc (t) options operator, arrow; display(line([0, 0, 0], [10*sin(t), 0, 0], colour = red), point([10*sin(t), 0, 0], colour = green, symbolsize = 50, symbol = solidsphere), line([0, 0, 0], [0, 10*sin(t+1.6), 0], colour = orange), point([0, 10*sin(t+1.6), 0], colour = green, symbolsize = 40, symbol = solidsphere), line([0, 0, 0], [0, 0, 10*sin(t+.5)^2], colour = yellow), point([0, 0, 10*sin(t+.5)^2], colour = green, symbolsize = 30, symbol = solidsphere)) end proc

animate(a1, [t], t = 0 .. 2*Pi, frames = 100)

 

c1 := Vector(3, [1, -.2, 0]):

c2 := Vector(3, [-1, -.2, 0]):

c3 := Vector(3, [0, 1, 0]):

c4 := Vector(3, [1, -.2, 0]):

c5 := Vector(3, [-1, -.2, 0]):

c6 := Vector(3, [0, 1, 0]):

c7 := Vector(3, [.5, 0, 0]):

c8 := Vector(3, [0, .5, 0]):

c9 := Vector(3, [0, 0, .5]):

c := [seq(c || n, n = 1 .. 9)]

[Vector[column](%id = 18446744074349586662), Vector[column](%id = 18446744074349586782), Vector[column](%id = 18446744074349586902), Vector[column](%id = 18446744074349587022), Vector[column](%id = 18446744074349587142), Vector[column](%id = 18446744074349587382), Vector[column](%id = 18446744074349587502), Vector[column](%id = 18446744074349587622), Vector[column](%id = 18446744074349587742)]

(3)

Rot := Matrix([[cos(w*t), -sin(w*t), 0], [sin(w*t), cos(w*t), 0], [0, 0, 1]])

Matrix(%id = 18446744074349582806)

(4)

d := seq(Rot.c[n], n = 1 .. 9)

Vector[column](%id = 18446744074349577870), Vector[column](%id = 18446744074349577990), Vector[column](%id = 18446744074349578110), Vector[column](%id = 18446744074349578230), Vector[column](%id = 18446744074349578350), Vector[column](%id = 18446744074349578590), Vector[column](%id = 18446744074349578710), Vector[column](%id = 18446744074349578830), Vector[column](%id = 18446744074349578950)

(5)

``

w := .25

.25

(6)

``

``

[d[1][1], d[1][2], d[1][3]]

[cos(.25*t)+.2*sin(.25*t), sin(.25*t)-.2*cos(.25*t), 0.]

(7)

ff1 := proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], [d[1][1], d[1][2], d[1][3]], colour = red, thickness = 4)]) end proc

proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], [d[1][1], d[1][2], d[1][3]], colour = red, thickness = 4)]) end proc

(8)

plots:-animate(ff1, [t], t = 0 .. 5, frames = 100)

 

 

Lend1 := [cos(w*t)+.2*sin(w*t), sin(w*t)-.2*cos(w*t), 0.]

[cos(.25*t)+.2*sin(.25*t), sin(.25*t)-.2*cos(.25*t), 0.]

(9)

ff2 := proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], Lend1, colour = red, thickness = 4)]) end proc

proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], Lend1, colour = red, thickness = 4)]) end proc

(10)

plots:-animate(ff2, [t], t = 0 .. 5, frames = 100)

 

``

ff3 := proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], [cos(.35*t)+.2*sin(.35*t), sin(.35*t)+(-1)*.2*cos(.35*t), 0.], colour = red, thickness = 4)]) end proc

proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], [cos(.35*t)+.2*sin(.35*t), sin(.35*t)-.2*cos(.35*t), 0.], colour = red, thickness = 4)]) end proc

(11)

plots:-animate(ff3, [t], t = 0 .. 5, frames = 100)

 

``


 

Download Anim_How-1.mw

restart

with(plots)

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, shadebetween, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

(1)

with(plottools)

[annulus, arc, arrow, circle, cone, cuboid, curve, cutin, cutout, cylinder, disk, dodecahedron, ellipse, ellipticArc, exportplot, extrude, getdata, hemisphere, hexahedron, homothety, hyperbola, icosahedron, importplot, line, octahedron, parallelepiped, pieslice, point, polygon, prism, project, rectangle, reflect, rotate, scale, sector, semitorus, sphere, stellate, tetrahedron, torus, transform, translate]

(2)

``

``

``

a1 := proc (t) options operator, arrow; display(line([0, 0, 0], [10*sin(t), 0, 0], colour = red), point([10*sin(t), 0, 0], colour = green, symbolsize = 50, symbol = solidsphere), line([0, 0, 0], [0, 10*sin(t+1.6), 0], colour = orange), point([0, 10*sin(t+1.6), 0], colour = green, symbolsize = 40, symbol = solidsphere), line([0, 0, 0], [0, 0, 10*sin(t+.5)^2], colour = yellow), point([0, 0, 10*sin(t+.5)^2], colour = green, symbolsize = 30, symbol = solidsphere)) end proc

animate(a1, [t], t = 0 .. 2*Pi, frames = 100)

 

c1 := Vector(3, [1, -.2, 0]):

c2 := Vector(3, [-1, -.2, 0]):

c3 := Vector(3, [0, 1, 0]):

c4 := Vector(3, [1, -.2, 0]):

c5 := Vector(3, [-1, -.2, 0]):

c6 := Vector(3, [0, 1, 0]):

c7 := Vector(3, [.5, 0, 0]):

c8 := Vector(3, [0, .5, 0]):

c9 := Vector(3, [0, 0, .5]):

c := [seq(c || n, n = 1 .. 9)]

[Vector[column](%id = 18446744074349586662), Vector[column](%id = 18446744074349586782), Vector[column](%id = 18446744074349586902), Vector[column](%id = 18446744074349587022), Vector[column](%id = 18446744074349587142), Vector[column](%id = 18446744074349587382), Vector[column](%id = 18446744074349587502), Vector[column](%id = 18446744074349587622), Vector[column](%id = 18446744074349587742)]

(3)

Rot := Matrix([[cos(w*t), -sin(w*t), 0], [sin(w*t), cos(w*t), 0], [0, 0, 1]])

Matrix(%id = 18446744074349582806)

(4)

d := seq(Rot.c[n], n = 1 .. 9)

Vector[column](%id = 18446744074349577870), Vector[column](%id = 18446744074349577990), Vector[column](%id = 18446744074349578110), Vector[column](%id = 18446744074349578230), Vector[column](%id = 18446744074349578350), Vector[column](%id = 18446744074349578590), Vector[column](%id = 18446744074349578710), Vector[column](%id = 18446744074349578830), Vector[column](%id = 18446744074349578950)

(5)

``

w := .25

.25

(6)

``

``

[d[1][1], d[1][2], d[1][3]]

[cos(.25*t)+.2*sin(.25*t), sin(.25*t)-.2*cos(.25*t), 0.]

(7)

ff1 := proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], [d[1][1], d[1][2], d[1][3]], colour = red, thickness = 4)]) end proc

proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], [d[1][1], d[1][2], d[1][3]], colour = red, thickness = 4)]) end proc

(8)

plots:-animate(ff1, [t], t = 0 .. 5, frames = 100)

 

 

Lend1 := [cos(w*t)+.2*sin(w*t), sin(w*t)-.2*cos(w*t), 0.]

[cos(.25*t)+.2*sin(.25*t), sin(.25*t)-.2*cos(.25*t), 0.]

(9)

ff2 := proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], Lend1, colour = red, thickness = 4)]) end proc

proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], Lend1, colour = red, thickness = 4)]) end proc

(10)

plots:-animate(ff2, [t], t = 0 .. 5, frames = 100)

 

``

ff3 := proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], [cos(.35*t)+.2*sin(.35*t), sin(.35*t)+(-1)*.2*cos(.35*t), 0.], colour = red, thickness = 4)]) end proc

proc (t) options operator, arrow; plots:-display([plottools:-line([0, 0, 0], [cos(.35*t)+.2*sin(.35*t), sin(.35*t)-.2*cos(.35*t), 0.], colour = red, thickness = 4)]) end proc

(11)

plots:-animate(ff3, [t], t = 0 .. 5, frames = 100)

 

``


 

Download Anim_How-1.mw

Connecting to a linux x2go server with an x2go client (tried Mac and Linux) and try to start maple2016 in a terminal, I get this error:

$ /opt/maple2016/bin/xmaple 
Picked up JAVA_TOOL_OPTIONS: 
Exception in thread "Request id 1" java.lang.UnsupportedOperationException: PERPIXEL_TRANSLUCENT translucency is not supported
    at java.awt.Window.setBackground(Window.java:3842)
    at java.awt.Frame.setBackground(Frame.java:988)
    at com.maplesoft.worksheet.application.WmiSplashScreen.<init>(Unknown Source)
    at com.maplesoft.worksheet.application.WmiGenericStartupStrategy.showSplash(Unknown Source)
    at com.maplesoft.worksheet.application.WmiGenericStartupStrategy.doStartup(Unknown Source)
    at com.maplesoft.worksheet.application.WmiWorksheetStartupStrategy.doStartup(Unknown Source)
    at com.maplesoft.application.Maple.doStartup(Unknown Source)
    at com.maplesoft.application.Application.startup(Unknown Source)
    at com.maplesoft.application.ServerProtocol$StartApplicationHandler.processCommand(Unknown Source)
    at com.maplesoft.application.ServerProtocol.executeCommand(Unknown Source)
    at com.maplesoft.application.ServerProtocol.processNextStep(Unknown Source)
    at com.maplesoft.application.ExchangeProtocol.executeProtocol(Unknown Source)
    at com.maplesoft.application.ApplicationManager$Listener.run(Unknown Source)
    at java.lang.Thread.run(Thread.java:745)

The same xmaple will start graphically if I connect with a normal remote X connection. However a normal remote X connection is unusable on a slow network connection.

 

 

I have 5 equations mentioned as below,

a[0]=L+(2*beta*mu)/(3*alpha*omega^2);

a[1]=sqrt((4*mu)/(3*alpha*omega^2));

a[2]=(2*beta*mu)/(9*alpha*omega^2);

a[3]=sqrt(mu^3/(432*alpha*omega^2));

omega[s]=omega-(mu^2/16*omega)-((2*beta^2mu)/(9*alpha*omega));

 

Left-hand parameters i.e. a[0],a[1],a[2],a[3], omega[s] are known , and each equation on right hand contain unkonowns i.e alpha, mu, beta, omega.

I want expression for unknowns in terms of knowns.

Hello,

Maple 2016 used to work fine on my PC but now it won't even load. Whenever I try to click the program I can see the loading icon next to the cursor and then nothing happens next. It works fine on my laptop and I'm using the program under a univsersity license. Yesterday I was able to open it once for some reason and now it's the same as before. I've tried uninstalling it again and again several time but that doesn't work.

Any help would be much appreaciated.

I have a problem in this form:

a1x1+a1x1+...anx1 +(x1)= b1

a1x2+a1x2+...anx2 +(x2)= b2

..........................................

a1xn+a1xn+...anxn +(xn)= bn

Here, x1,x2,.......,xn  are unknowns, a1,a2,...,aand b1,b2,...,bare known contants.

Is there a possible way to find unknowns?

 

Hello,

How to build cuboid in pyramid like this animation

http://www.irem.univ-mrs.fr/article/optimisation-du-volume-d-un-pave.html

Thanks !

Hello!

I was creating this code for the "Müller Method" for Numerical Analysis. Everything works fine till the if (for getting the x3) goes on. The problem is, the Maple doesn't detects the if.
I tried checking every variable and every variable is calculated except "x3"; even "disc".

I would like to know what happened. Anyways, here's the code:
http://imgur.com/a/DGlgl
 

                     
x0 := 0; x1 := 8.4; x2 := 10; iter := 1000; tol := 10^(-8); f := proc (x) options operator, arrow; 3*x^3+7*x^2+x+2 end proc; f(x); plot(f(x), color = green); printf(" n          x0               x1               x2              x3                Error \n"); for i to iter do d0 := evalf((f(x1)-f(x0))/(x1-x0)); d1 := evalf((f(x2)-f(x1))/(x2-x1)); h1 := evalf(x2-x1); h0 := evalf(x1-x0); a := evalf((d1-d0)/(h1-h0)); b := evalf(a*h1+d1); c := evalf(f(x2)); disc := sqrt(-4*a*c+b^2); if abs(b+disc) > abs(b-disc) then x3 := x2+(-2*c)*(1/(b+sqrt(-4*a*c+b^2))); erry := abs((x3-x2)/x3) else x3 := x2+(-2*c)*(1/(b-sqrt(-4*a*c+b^2))); erry := abs((x3-x2)/x3) end if; if erry > tol then x0 := x1; x1 := x2; x2 := x3; printf("%2d     %2.8f      %2.8f       %2.8f         %2.5 f     %2.8 f \n", i, a, b, c, x3, erry) else printf("una raiz es: %2.8f ", x3); break end if end do;

I put the image (as it looks on my maple) and the "code" so you can copy-paste it in Maple.

Dear hope you will be fine I want to represent the data in term of solid line not dotted, asterik etc of the following code. Please help me to fix this problem

 

restart; epsilon := 0; Pr := 1; beta := .1; Sc := 1; S := 0; L := 15;
                               15
for i from -L while i <= L do a[i] := 1.0*i/L end do;
for i2 from -L while i2 <= L do fw := a[i2]; Eq1[i2] := eval(diff(F(eta), eta, eta, eta)+F(eta)*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2+S*(epsilon-(diff(F(eta), eta)))+epsilon^2); Eq2[i2] := eval((diff(G(eta), eta, eta))/Pr-G(eta)*(diff(F(eta), eta))+F(eta)*(diff(G(eta), eta))); Eq3[i2] := eval(diff(H(eta), eta, eta)+Sc*(F(eta)*(diff(H(eta), eta))-beta*H(eta))); IC[i2] := F(0) = a[i2], (D(F))(0) = 1, (D(F))(L) = epsilon, G(0) = 1, G(L) = 0, H(0) = 1, H(L) = 0; dsys1[i2] := {Eq1[i2], Eq2[i2], Eq3[i2], IC[i2]}; dsol1[i2] := dsolve(dsys1[i2], numeric, output = listprocedure, range = 0 .. L); dsol1x[i2] := subs(dsol1[i2], diff(F(eta), eta, eta)); dsol1y[i2] := subs(dsol1[i2], G(eta)); dsol1z[i2] := subs(dsol1[i2], H(eta)) end do;
for j from -L while j <= L do g[j] := eval(-dsol1x[j](0)) end do;
with(plots);

g1 := pointplot({seq([n/L, g[n]], n = -L .. L)}, symbol = asterisk, symbolsize = 15, color = blue);
display(g1);
 

Hello everyone,

I have 5 equations (fa,fb,fc,fk,fv) and 5 variables (V0,A0,A1,A2,k1) and I want to solve them numerically. Problem is my inability to set properly intervals and starting values. 

Here is my command: fsolve({fa,fb,fc,fk,fv},{V0,A0,A1,A2,k1});

If I run it just like this, I get some values, lets say V0=0.00045 etc.But when I set range for V0 like this:

fsolve({fa,fb,fc,fk,fv},{V0=0.0004..0.0005,A0,A1,A2,k1});

I get this:

Error, (in fsolve) fsolve cannot solve on V0 = 0.4e-3 .. 0.5e-3

I used interval 0.0004..0.0005 only to prove that there must be the wrong syntax, because obviously 0.0004<0.00045<0.0005. In reality i need interval 0.0001..0.0002 and it is also necessary.to set some initial values (A0=0.0023) but first I need to solve my problem with syntax.

Any advice ?

hello i have the following set of ode's:

ode_sub := diff(S(t), t) = -k1*S(t)-S(t)/T1_s;
ode_P1 := diff(P1(t), t) = k1*S(t)-k2*(P1(t)-P2(t)/keq)-P1(t)/T1_p1;
ode_P2 := diff(P2(t), t) = -k2*(-keq*P1(t)+P2(t))/keq-k4*P2(t)-P2(t)/T1_p2;
ode_P2e := diff(P2_e(t), t) = k4*P2(t)-P2_e(t)/T1_p2_e;

ode_system := ode_sub, ode_P1, ode_P2, ode_P2e;

with these parameters:
s0 := 10000;
k2 := 1000; T1_s := 14; T1_p2_e := 35; T1_p2 := T1_p1;
 

i want to find the unkown parameters : T1_p1, k1, keq and k4

my idea was this:

init:=S(0)=s0,P1(0)=0,P2(0)=0,P2_e(0)=0

dsolve({ode_system,init})

sol := combine(expand(%));
PS := subs(sol, [S(t), P1(t), P2(t), P2_e(t)]);
 

P1fu := unapply(PS[2],t);
Sfu := unapply(PS[1],t);
P2fu := unapply(PS[3],t);
P2e_fu := unapply(PS[4],t);
P2_total := unapply(P2fu+P2e_fu, t);
 

the following data is given:

T:=<0,2,4,6,8>

S:=<9999.99913146527,8328.870587730016,6937.009129218748,5777.745632133724,4812.209983843559>

P1:=<0.0,67.86790056712294,114.88787098501874,145.95438088662502,164.85650644237887>

P2_P2e:=<0.0,271.68492651947497,461.9130396605823,589.3710176125417,668.9967533337124> # data from P2(t)+P2_e(t)

 

making the rediduals:

RP1 := convert(P1-P1fu~(T), list);
RS := convert(S-Sfu~(T), list);
RP2_P2e := convert(P2_P2_e-P2_total~(T), list);
 

RPs := [op(RS), op(RP2_P2_e), op(RP1)]

res := Optimization:-LSSolve(RPs, k1 = 0 .. 1, keq = 0 .. 10, k4 = 0 .. 1, T1_p1 = 0 .. 100)

i dont know wheter or not the last step work to get the parameters becuase it takes to long to compute. is there a smarter way to obtain the parameters of the ode's? a numeric approch ?

i tried with dsolve({ode_sysytem,init},numeric,'parameters'=[k1,keq,k4,T1_p1]) however it doesnt seem to get my anywhere since i need to know the parameters to use this (i think)

hope someone can help:)

 

 

Hello everyone,

I am making a mobile robot simulation in the MapleSim environment and I need to create varriable (uneven) terrain, do I need an additional library for that?

If not how can I do so?

Thanks for your time!

John

how to convert a^2*b+c to func2(func1(func1(abc[1],abc[1]),abc[2]),abc[3])

when i use custom function func2 to represent plus, func1 to represent multiply

input 3 parameters,

one is a^2*b + c one is [func1, func2] and second is [abc[1],abc[2],abc[3]] corresponding to a, b, c
a^2*b + c = func2(func1(func1(abc[1],abc[1]),abc[2]),abc[3]);

This is simplified from a larger example.  I never saw Maple exit from a numercal calculation.

> 123456789 ^ 987654321:
Execution stopped: Stack limit reached.

Process Maple R2016 exited abnormally with code 158

> kernelopts(maxdigits);
                                  38654705646

> lprint(123456789. ^ 987654321.);
.4339566080e7991619731

kernelopts(memusage) shows about 2~3 meg bytes used.

The system has 16G Ram.  16G swap.

Seems like something maple should have caught.

First 1010 1011 1012 1013 1014 1015 1016 Last Page 1012 of 2429