Maple 2018 Questions and Posts

These are Posts and Questions associated with the product, Maple 2018

Please check this code and fix the error. And help me with the below two questions.

1- In the above image, I'm getting the error when I'm to plot my solutions.

2- How to get information from this system?

      The values of x_1 and corresponding y_1 and so on and errors etc

This is a downloadable link to my worksheet [Practice.mw]

Thanks!


>restart;
> ODEtools[declare]((S, L, B)(t), prime = t);
         ODEtools[declare](S(t), L(t), B(t), prime = t)
> DE1 := diff(S(t), t) = -L*S*beta-`μS`+mu;
                    d                               
            DE1 := --- S(t) = -L S beta - μS + mu
                    dt                              
> DE2 := diff(L(t), t) = L*S*beta-(gamma+mu)*S;
                  d                                  
          DE2 := --- L(t) = beta S L - (gamma + mu) S
                  dt                                 
> DE3 := diff(B(t), t) = L*gamma-gamma*mu;
                      d                           
              DE3 := --- B(t) = L gamma - gamma mu
                      dt                          
> init_conds := S(0) >= 0, L(0) >= 0, B(0) >= 0;
         init_conds := 0 <= S(0), 0 <= L(0), 0 <= B(0)
> sys := {init_conds, diff(B(t), t) = L*gamma-gamma*mu, diff(L(t), t) = L*S*beta-(gamma+mu)*S, diff(S(t), t) = -L*S*beta-`&mu;S`+mu};

> sol := dsolve(sys, numeric, parameters = [mu, beta, gamma, S(t), L(t), B(t)], method = rkf45);

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations


> sol(parameters = [mu = .234, beta = 2.345, gamma = 5.678, S(t) = 6.678, L(t) = 6.789, B(t) = 7.123]);
           sol(parameters = [mu = 0.234, beta = 2.345, gamma = 5.678, S(t) = 6.678, L(t) = 6.789, B(t) = 7.123])

On CBS tv show Survivor S38E4, an interesting scenario arose. There was initially 2 teams comprising 6 persons (blue team) and 9 persons (yellow team). Then they decided to mix it up, and split it into 3 teams. Each player was randomly assigned a buff, (red, green or orange =3 colored buffs *5).= 3 teams of 5.

What was interesting was 5 out of the original 6 blue team are now in the new red team. 5 of the original yellow team are in the new green team and 4 of the original yellow team are in the new orange team. 

What is the chances of that happening?
A couple of attempts:

survivor_rand.mw

I'm using Maple on a compute server (via Putty from my computer which has windows system). It is an only-text mode server so it is a black screen and text and no menu bar to choose file and then save by clicking etc. You only can write commands and press Enter.

The IT section suggested to use writeto("filename.txt"), but I can't see any file on my drive, I used writeto("\drvename\filename.txt"), but still nothing.

Thanks to @acer@Joe Riel and the IT section. The problem now is solved. It was about the way that I was writing the directory path. I had to add something before the "drivename" in my case. Plus that it is a Unix path, so I should use / not \.

Hello, i am experiencing some problems when trying to open the maple 2018 software*
I have tried unistalling and download it again.
I have tried to search for sollution but there is very ittle intel
When i open Maple 2018 it just lingeres on the start up (pic below) and just disappears after 10 seconds

Can someone please help i have a very important examination upcoming

restart;
F0:=proc(sigma__xx,N)
local  x,y,Fx,Fy,:
assume (w,real,w>0):assume (h,real,h>0):

for n from 0 by 2 to N do Fx:=integrate(sigma__xx*cos(w*Zeta*h),Zeta=0..infinity):
end do;
for n from 1 by 2 to N do Fx:=integrate(sigma__xx*sin(w*Zeta*h),Zeta=0..infinity):
end do;
return [Fx]:
end proc;

sigma__xx := -(sqrt(Zeta^2*h^2+h^2)^(-n+2)*cos(n*arctan(h, Zeta*h))*n^2+sqrt(Zeta^2*h^2+h^2)^(-n+2)*cos(n*arctan(h, Zeta*h))*n-2*sqrt(Zeta^2*h^2+h^2)^(-n+2)*cos(n*arctan(h, Zeta*h)))*Zeta^2*h^2/(Zeta^2*h^2+h^2)^2+(sqrt(Zeta^2*h^2+h^2)^(-n+2)*cos(n*arctan(h, Zeta*h))*n^2-3*sqrt(Zeta^2*h^2+h^2)^(-n+2)*cos(n*arctan(h, Zeta*h))*n;

F0(sigma__xx,N);

Hello everyone,

I am a studen and have an astrophysics class where I programmed a basic space probe to go from a planet A to a Planet B.

I have a 2D grap and a 3D graph already done however I would like to animate them.

I have included both a picture and my whole program. Also see some code below. 

Thanks to all.

 

Download SondeII_(1).mwSondeII_(1).mwSondeII_(1).mw

restart:

with(linalg):
with(DEtools):

#Sonde
Position := [x(t), y(t)]:

#Terre
omega1 :=2*Pi:
r1 := [cos(omega1*t), sin(omega1*t)]:
x1(t) := innerprod([1, 0], r1):
y1(t) := innerprod([0, 1], r1):

#Mars
omega2 := sqrt(2)*Pi/2:
phi2 := Pi/2:
r2 := [2*cos(omega2*t + phi2), 2*sin(omega2*t+ phi2)]:
x2(t):= innerprod([1, 0], r2):
y2(t):= innerprod([0, 1], r2):

#Les couplages gravitationnelles (masse).
c1 := 1.0:
c2 := 0.2:
c0 := 100:

#Les forces appliqués sur la sonde
ForceGravitationnelle1 := -c1*(Position-r1)/(sqrt((x(t)-x1(t))^2+(y(t)-y1(t))^2))^3:
ForceGravitationnelle2 := -c2*(Position-r2)/(sqrt((x(t)-x2(t))^2+(y(t)-y2(t))^2))^3:
ForceGravitationnelle0 := -c0*(Position)/(sqrt((x(t))^2+(y(t))^2))^3:

#La somme des forces.
Force := ForceGravitationnelle1 + ForceGravitationnelle2 + ForceGravitationnelle0:

Fx := innerprod([1, 0], Force):
Fy := innerprod([0, 1], Force):

#L'interval de temps.
TempsInit := 0:
TempsFinal := 3:

#Les équations différentielles de deuxieme ordre.
eq1x := (D(D(x)))(t) = Fx:
eq1y := (D(D(y)))(t) = Fy:

#Les conditions initiales..
phi0 :=(Pi)/2:
V0 := 12.946802:
x0 := 1:
y0 := 0.1:
Vx0 := V0*cos(phi0):
Vy0 := V0*sin(phi0):

ConditionsInit := x(0) = x0, y(0) = y0, D(x)(0) = Vx0, D(y)(0) = Vy0:

#La trajectoire de la sonde.
Trajectoire := dsolve({eq1x, eq1y, ConditionsInit}, {x(t), y(t)}, numeric, range = TempsInit..TempsFinal, maxfun=0):

#Tracage du graphique de la trajectoire en 2D
plots[odeplot](Trajectoire, [[0,0],[x1(t),y1(t)],[x2(t), y2(t)], [x(t), y(t)]],
TempsInit..TempsFinal, numpoints = 1000, axes = boxed, scaling = constrained, thickness = [2],
color = ["Black", "Green", "Blue", "Red"],
labels = ["X (L)", "Y (L)"],
labelfont = ["Times", 14], title = "Mouvement de la sonde dans le plan",
titlefont = ["Helvetica", 14], style=[point,line,line,line], symbol = solidcircle);

#Tracage du graphique en 3D:
plots[odeplot](Trajectoire, [[0,0,t],[x1(t),y1(t), t],[x2(t), y2(t), t], [x(t), y(t), t]],
TempsInit..TempsFinal, numpoints = 1000, axes = boxed, scaling = constrained, thickness = [3],
color = ["Black", "Green", "Blue", "Red"],
labels = ["X (L)", "Y (L)", "t"],
labelfont = ["Times", 14], title = "Mouvement de la sonde dans le plan",
titlefont = ["Helvetica", 14], style=[point,line,line,line], symbol = solidcircle);
 

 

 

 

 

 


 

Hello,

How I can take Laplace inverse?

Thank you

LAPLACE

Maple Worksheet - Error

Failed to load the worksheet //convert/LAPLACE
 

Download LAPLACE

 

Give the following functions find Domain, Range, Possible Asymptotes, Intercepts, Critical Points, Intervals of Increase, Decrease, Relative and Absolute Extrema, and Concavity.

A) f(x)=x(x^2-6x+8)

B) f(x) =x^3/4 -3x

After two days of trying to take the Jacobian of a multi-variable mapping in Maple, I'm still finding no way. Here is an example of what I mean by "multi-variable mapping":

(x,y) -> x^2-y^2

Every AP high-school student knows how to take the partial derivatives of that, but Maple's "Jacobian" feature doesn't work on mappings. It only works on expressions.

So I spent today trying to leverage off the "D" feature, only to find it rendered useless by bugs: https://mapleprimes.com/questions/227007-How-To-Take-Partial-Derivative-Of-A-Mapping

It took me less than two days to do it in Mathematica, by way of comparison.

On the website, they show a photo of a man elated to be doing math in Maple, but that's not how it's gone for me. I'm cursing outloud, and my blood pressure is high. I'm angry that I lost money on a product that can't take the Jacobian of a mapping, and I'm angry that I lost two days of my life doing tedious debugging while paying money instead of getting paid for my work.

Maple pays math PhD's to release a commercial CAS that can't take the Jacobian of a mapping. So much for the value of a PhD.

I need to find the intersection of 3 planes.

x+3*y-5*z=0
x+4*y-8*z=0
-2*x-7*y+13*z=0

I managed to get the intersection using LinearSolve but I keep getting an error when I try to plot out the planes using plot3d. What can I change in the commands for it to work?
 

with(LinearAlgebra)

[`&x`, Add, Adjoint, BackwardSubstitute, BandMatrix, Basis, BezoutMatrix, BidiagonalForm, BilinearForm, CARE, CharacteristicMatrix, CharacteristicPolynomial, Column, ColumnDimension, ColumnOperation, ColumnSpace, CompanionMatrix, CompressedSparseForm, ConditionNumber, ConstantMatrix, ConstantVector, Copy, CreatePermutation, CrossProduct, DARE, DeleteColumn, DeleteRow, Determinant, Diagonal, DiagonalMatrix, Dimension, Dimensions, DotProduct, EigenConditionNumbers, Eigenvalues, Eigenvectors, Equal, ForwardSubstitute, FrobeniusForm, FromCompressedSparseForm, FromSplitForm, GaussianElimination, GenerateEquations, GenerateMatrix, Generic, GetResultDataType, GetResultShape, GivensRotationMatrix, GramSchmidt, HankelMatrix, HermiteForm, HermitianTranspose, HessenbergForm, HilbertMatrix, HouseholderMatrix, IdentityMatrix, IntersectionBasis, IsDefinite, IsOrthogonal, IsSimilar, IsUnitary, JordanBlockMatrix, JordanForm, KroneckerProduct, LA_Main, LUDecomposition, LeastSquares, LinearSolve, LyapunovSolve, Map, Map2, MatrixAdd, MatrixExponential, MatrixFunction, MatrixInverse, MatrixMatrixMultiply, MatrixNorm, MatrixPower, MatrixScalarMultiply, MatrixVectorMultiply, MinimalPolynomial, Minor, Modular, Multiply, NoUserValue, Norm, Normalize, NullSpace, OuterProductMatrix, Permanent, Pivot, PopovForm, ProjectionMatrix, QRDecomposition, RandomMatrix, RandomVector, Rank, RationalCanonicalForm, ReducedRowEchelonForm, Row, RowDimension, RowOperation, RowSpace, ScalarMatrix, ScalarMultiply, ScalarVector, SchurForm, SingularValues, SmithForm, SplitForm, StronglyConnectedBlocks, SubMatrix, SubVector, SumBasis, SylvesterMatrix, SylvesterSolve, ToeplitzMatrix, Trace, Transpose, TridiagonalForm, UnitVector, VandermondeMatrix, VectorAdd, VectorAngle, VectorMatrixMultiply, VectorNorm, VectorScalarMultiply, ZeroMatrix, ZeroVector, Zip]

(1)

A := `<|>`(`<,>`(1, 1, -2), `<,>`(3, 4, -7), `<,>`(-5, -8, 13))

Matrix(%id = 18446745896624469046)

(2)

b := `<,>`(0, 0, 0)

Vector[column](%id = 18446745896624462294)

(3)

x := LinearSolve(A, b)

Vector[column](%id = 18446745896652349678)

(4)

P1 := x+3*y

P2 := x+4*y

P3 := -2*x-7*y

plot3d([P1, P2, P3], x = -8 .. 8, y = -20 .. 20, plotlist = true, color = [blue, red, green])

Error, (in plot3d) unexpected options: [(Vector(3, {(1) = -4*_t[3], (2) = 3*_t[3], (3) = _t[3]})) = -8 .. 8, y = -20 .. 20]

 

``


 

Download problem.mw

Here's how to break Maple's 'D' function:

~$ cmaple mapleBug.mpl
    |\^/|     Maple 2018 (X86 64 WINDOWS)
._|\|   |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc. 2018
 \  MAPLE  /  All rights reserved. Maple is a trademark of
 <____ ____>  Waterloo Maple Inc.
      |       Type ? for help.
> D[1]((x, y) ->x^2-y^2);
                                 (x, y) -> 2 x

> D[2]((x, y) ->x^2-y^2);
                                 (x, y) -> -2 y

> with(VectorCalculus):
> assume(x, 'real');
> assume(y, 'real');
> D[1]((x, y) ->x^2-y^2);
                                  (x, y) -> 0

> D[2]((x, y) ->x^2-y^2);
                                  (x, y) -> 0

> quit
memory used=2.4MB, alloc=8.3MB, time=0.06

~$

What if a person needs to have these lines:

> with(VectorCalculus):
> assume(x, 'real');
> assume(y, 'real');

and then he wants to take partial derivatives? Not possible? Once you load VectorCalculus and declare x and y real, Maple stops doing partial derivatives? Why.

Sorry,

I have been away from Maple for a year.
Then, when I used op command, I am puzzed to notice the results were different from those I know.

>op((x+5)^2*(x+y));
      _EXPSEQ((x+5)^2, x+y)

Result I know is 

        (x+5)^2, x+y

Has a modifire such as _EXPSEQ   automatically come to be attatched?
Or, can I have maple express it  in the form I know:  (x+5)^2, x+y?

Thank you in advance.

taro

 

 

 

 

 


 

Given the following functions and respective intervals graph them and determine all values of in the interval (a,b) such that

f'(c) =f(b)-f(a)/b-a (apply the Mean Value Theorem) 

 Question 1: f(x)=x^3-2*x      [0,2]

 Question 2 : g(x)= cuberoot(x-3)^2     [-3,4]

 

Please HELP!!! 

 

how I can write a program code for newmark method.

in this method time has 3 order derivation

First 6 7 8 9 10 11 12 Last Page 8 of 41