Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi:

For many reasons, until recently, I have been always using Classic interface under Maple through Maple 16. When trying Maple 16 with Java interface I have noticed several quirks which were not present under the Classic interface of earlier versions through (at least) Maple 13. One of them is this:

The shortcut keys such as CTRL+Insert and SHIFT+Insert do not work and I am forced to use a mouse and copy/paste icons in Maple. This really slows down my keyboard...

 

Hi

 

I want to solve the equation of

F'''(x)+F(x)*F''(x)=0

with the following  boundary conditions

F(-5)=F'(-5)=0

F'(5)=0.14

for solving it numerically in the computational domain of [-5...5]. I have seprated the domain into two parts and the governing equation and boundary condition in each domain is assumed as follows


g'''(x)+g(x)*g''(x)=0
g(0)=alpha
g(-5)=0
g'(-5)=0

Here is my code:

restart:
with(LinearAlgebra):with(plots):with(ExcelTools):
n:=0.8:    e:=1:   h2:=0.01:
for N from 20 by 10 to 100 do
for b from 0.01 by 0.05 to 0.26 do
for h from 0.02 by 0.02 to 0.2 do
for h1 from 0.015 by 0.015 to 0.105 do
  A1:=proc(N::posint,b1,e1,h0,h11,h21) local M;
      M:=Matrix(N,(i,j)->if i=j then e elif abs(i-j)=1 then h else 0 fi);
      M[1,2...

When using "plotsetup(png,...)" the third rotational angle is ignored.  Attached are a worksheet and the corresponding three plots that demonstrate this bug.

Dear everyone,

 

  Hello,

 

  I have a question about how to get the power 'n' in an expression "(x+y)^n". For instance, I have a Maple code

 

******

eq1:=4*x^n*y^m+(x+y)^n+x^3;

f1:=map(t -> `if`(match(t = c*y^d, y, 's'), subs(s,d), NULL),
      convert(eq1, list));
f2:=map(t -> `if`(match(t = c*x^d, x, 's'), subs(s,d), NULL),
      convert(eq1, list));

It will be possible to manipulate this behavior with options "dsolve / event" or "dsolve / Error Control"?

> restart;
> Digits;                                10
> ode := sin(x(t))*(diff(x(t), t, t))+x(t)^3 = t*sin(x(t));

> ics := x(0) = 0, (D(x))(0) = 0;

> s := 250;

The explosion of cloud this - cloud that, is rather obvious now.

 

SO, I decided to check-out the cloud feature to my Maple? One problem, I cannot see - what is in it?

  • Yes, I understand, you use the cloud palette in Maple.
  • It shows the group, who the document / notebook is from, what it's title is? And such.
  • _
  • BUT, is there a web page - listing? hidden somewhere, on MaplePrimes or someplace on Maplesoft's service? (

Hi Guys,

I was wondering if there is a way to get data from an array solution.

The system of odes, I am trying to solve is

eq2:= diff(f(eta),`$`(eta,3))+f(eta)*diff(f(eta),`$`(eta,2))+beta*(1-diff(f(eta),eta)^2) = 0;

eq2 := diff(theta(eta),`$`(eta,2))+pr(f(eta)+(beta-1)*
eta*diff(f(eta),eta))*diff(theta(eta),eta) = 0;

bcs:=f(0)=gama,D(f)(0)=lambda,D(f)(N)=1,theta(0)=0,theta(N)=1:

gama:=2;lambda:=1:pr:=10:N:=10:beta:=3:

I would like to draw some areas between curves. I tried to use ?plots/implicitplot but I often need to use gridrefine option which consumes so much memory. When I want to draw e.g. the area between f and g:

f := x-> piecewise(x < 5, 2, x-3):
g := x-> piecewise(x < 4, x, 4):

I can do:

p:=plot({f, g}, color...

I am making a 2D plot and need to mark certain irregular points on the x-axis, so I managed to put tickmarks (and gridlines) at those points, and somehow Maple even labels each tickmark, which is what I want. 

Except, there is no room along the x-axis for all the labels, so I need the labels to appear vertically at each tickmark.

Is  there any way to do this? The only instructions I can find for vertical labelling just apply to axis labels, not tickmark or gridline labels.


Hola,
How can leave a trail/trace of black point in the yellow and purple screens? Something very similar to the option "animate / trace."

SVD-filtering.mw

Make it possible to take the numerical data of each screen from the "plottools [getdata]"

Gracias

 

seriestoratpoly(series(simplify(LaguerreL(5, x)),x,6),['egf']);

seriestoratpoly(series(simplify(LaguerreL(5, a, x)),x,6),['egf']);

inverse z transform of sqrt(1-a^2)/(1-a/z)

restart;
alpha := 0;
s := x;
tau := alpha + 1 - x;
source := s*Diff(y(x),x$2)+tau*Diff(y(x),x)+n*y(x)=0;
sol := dsolve(source,y(x),method=laplace);
zsol := ztrans(sol,x,z);
convert(zsol, ratpoly);

sol := y(x) = _C1*invlaplace((_s1-1)^n*_s1^(-1-n), _s1, x)
zsol := ztrans(y(x), x, z) = _C1*invlaplace((_s1-1)^n*_s1^(-1-n)*z/(z-1), _s1, x)
Error, (in convert/ratpoly) expecting a Laurent series or a Chebyshev series

what are _C1 and _s1?

Maple 16:

Given Matrix A of size nRow,nCol, and a vector b of size nCol, I wanted to divide each value in column 
K by b(K) where k=1..nCol and the result will be a matrix B also of size nRow,nCol ofcourse.

For example, I want to do A[..,1]/b[1],  A[..,2]/b[2], etc.. for each column of A

The way I do it now is by using seq(), then convert the result to matrix, like this:


nCol := LinearAlgebra:-ColumnDimension(A):
B    := convert([seq(A[..,i]/b[i],i=1..nCol)],Matrix);

And it works.  But was wondering if there is a better way to do this. I tried to use apply map but did 
not know how to
Here is below the complete Maple 16 code I have:

------------------------
A := LinearAlgebra:-RandomMatrix(3,4);
b := MTM:-sum(A,1);
nCol := LinearAlgebra:-ColumnDimension(A):
B := convert([seq(A[..,i]/b[i],i=1..nCol)],Matrix);
MTM:-sum(B,1);  #verify the sum is 1

When I execute several times my maple code the program completely freezes, even though in the beginning of the code teh command "restart" is inserted. I work with matrices, perform integrations and derivations for their components and calculate eigenvalue problems. The freeze happens at no particulat block of execution.

 

Any thoughts?

First 1594 1595 1596 1597 1598 1599 1600 Last Page 1596 of 2256