Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

how I can determined time period?

thank you

period.mw
 

d := (10+20*cos(Omega*t)+30*cos(9*sqrt(2)*t))^2

(10+20*cos(Omega*t)+30*cos(9*2^(1/2)*t))^2

(1)

with(StringTools)

period(d)

period((10+20*cos(Omega*t)+30*cos(9*2^(1/2)*t))^2)

(2)

``


 

Download period.mw

 


I get the following errors when attempting to use the Sockets package to interface with the serial input and output for a USB device connected and reported to have no known problems by Windows 10:

with(Sockets):
LookupService("busboy");
Error, could not determine determine port number for service "busboy"

LookupService(998);

Errror, cannot  determine "tcp"  service on port 998




server :=
proc (sid)
Sockets:-Write(sid, sprintf("Hello %s on port %d, from %s\r\n", Sockets:-GetPeerHost(sid), Sockets:-GetPeerPort(sid), Sockets:-GetHostName()))
end proc;

Sockets:-Serve(GetPeerPort(sid), server);
 
Error, (in Sockets:-GetPeerHost) Unknown error

sid := Open("localhost", "echo");

Sockets:-Serve(GetPeerPort(sid), server);

Error, (in Sockets:-Serve) cannot bind address: Unknown error
 

 

How can we solve the following pde by Maple? 

where v is velocity, v with dot is acceleration. (So, I think we will assume that acceleration is fixed.) And \delta is Dirac distribution.  E,I,m, M , g are fixed numbers.

Boundary conditions are:

Initial conditions are:

 

You can find the equation in the code: question.mw

Hi,

this "sum(1/(1+x)^t, t=1..infinity)" is (in my opionon) one of the most standard infinity summation and has the closed form 1/x. 
i used maple 18 and it was executed and i got the closed form 1/x.

with maple 2018 i get the non executed form, also the same inert form as "Sum(1/(1+x)^t, t=1..infinity)".

if i try "sum(1/(1+2)^t, t=1..infinity)" i get 1/2 as result.

why does the version above not working? any ideas?

thank you.

Dear friends,

I'm trying to solve a linear system of PDEs. After applying the casesplit command Maple returns: 

diff(U, z, y, x, x) = 0; diff(G, y,y,z) = diff(U,x,y,z);  diff(G,x,y,z) = diff(U,x,x,z)     

(with U(x,y,z), G(x,y,z)) 

The solution to the first equation is U = F4(x, y)+ F3(x, z) + F2(y, z) + x*F1(y, z). 

However, given this solution, I cannot satisfy diff(G,x,y,z) = diff(U,x,x,z).  What could I be doing wrong? 

Many thanks for your help.  

Hi

I have a problem with my maple. 

 

corrupt file 

Is there anybody who can help me solve this problem. 

Casper 

sys:={x+y=10,x^2-y^2+z^2=1};

{x+y = 10, x^2-y^2+z^2 = 1}

(1)

isolve(sys, n);

{x = -(1/20)*(1+20*n)^2+101/20, y = (1/20)*(1+20*n)^2+99/20, z = 1+20*n}

(2)

_SolutionsMayBeLost;  # It was set in Maple 10 or so!  E.g. x=5,y=5,z=-1 is not found.

_SolutionsMayBeLost

(3)

# Workaround

yz:=eliminate(sys,x);

[{x = -y+10}, {z^2-20*y+99}]

(4)

S:=isolve(yz[2], n);

{y = 20*n^2-38*n+23, z = 19-20*n}, {y = 20*n^2-22*n+11, z = 11-20*n}, {y = 20*n^2-18*n+9, z = 9-20*n}, {y = 20*n^2-2*n+5, z = 1-20*n}

(5)

map(u -> ( u union eval(yz[1],u) ), [S])[];  # The correct solution

{x = -20*n^2+38*n-13, y = 20*n^2-38*n+23, z = 19-20*n}, {x = -20*n^2+22*n-1, y = 20*n^2-22*n+11, z = 11-20*n}, {x = -20*n^2+18*n+1, y = 20*n^2-18*n+9, z = 9-20*n}, {x = -20*n^2+2*n+5, y = 20*n^2-2*n+5, z = 1-20*n}

(6)

 

This may seem a bit trivial, but I prefer f'(x) to writing diff(f(x),x) in 1D input. How to achieve?

differential.mw


 

15

 

"maple init loaded..."

(1)

In Document mode, this works fine.

f := proc (x) options operator, arrow; x^2 end proc

proc (x) options operator, arrow; x^2 end proc

(2)

diff(f(x), x)

2*x

(3)

But I mainly use Worksheet (1D) mode, and I can't seem to acheve the same, without using diff(f(x),x)

``

f:x->x^2

proc (x) options operator, arrow; x^2 end proc

(4)

f'(x)

Error, unexpected single forward quote

 

``


 

Download differential.mw

 

Hello,

Simple question:

Why does

restart;

(exp(I*t))^q;

`assuming`([simplify(%)], [t > 0, t < 2*Pi, q > 0, q < 1])

 

 

not simplify to exp(I*t*q)

?

does anyone khow to force maple to automatically simplify expressions that look like this:

restart:

cc:=arctan(sin(phi)/cos(phi))-arccos(cos(phi))      
simplify(cc) assuming phi::real; #should be zero

i've played with convert/expand/simplify and in certain cases of course i can force maple to do this... (kind of manually).. but was hoping this could be more automatic.

thanks!

 

The round(x) statement rounds x to the nearest integer but if x has unirs it appears not to work. Known bug?

RoundUnits.mw
 

restart

a := 15; b := 5.3

c := (1/100)*a+b

5.450000000

(1)

c := round(c)

5

(2)

``

A := 15*Unit('cm'); B := 5.3*Unit('m')

C := A+B

15*Units:-Unit(cm)+5.3*Units:-Unit(m)

(3)

C := simplify(C)

5.450000000*Units:-Unit(m)

(4)

Cr := round(C)

round(5.450000000*Units:-Unit(m))

(5)

``


 

Download RoundUnits.mw

 

Hello, I am getting the following output from maple: (-ln(lambda)-gamma-ln(k+b))/(k+b) . I have all variables (lambda, k, b) but not gamma and I am not sure what actually it is. I believe it is some kind of Gamma function but I cannot find any expressions for that. Ussually for gamma function I get something like GAMMA(x). Does someone know what this lower case gamma is?

I wish to write a simple procedure to evaluate the Poisson quantile function, F, for many possible parameter values, lambda.

The Maple commands to evaluate F for individual lambda values works just fine, however, I have tried to write a simple procedure to evaluate F for prescribed lambda values (imported from Excel) but to no avail. I'm missing something quite basic, I'm sure.

Can anybody offer a suggestion please? Thanks.

 

Inverse_Poisson_Procedure.mw

i want to gain diff(p(t), t) and diff(q(t), t) and Jacobian matrix
 according to the attached pdf file.

please help me.

thanks

simplify.mw
 

k := diff(a(t), t) = -mu*a(t)-(1/4)*alpha6*a(t)*sin(gamma(t))

diff(a(t), t) = -mu*a(t)-(1/4)*alpha6*a(t)*sin(gamma(t))

(1)

j := a(t)*(diff(gamma(t), t)) = 2*a(t)*sigma-(6*(1/8))*(alpha1-alpha2+(1/3)*alpha3)*a(t)^3-(1/2)*alpha6*a(t)*cos(gamma(t))

a(t)*(diff(gamma(t), t)) = 2*a(t)*sigma-(3/4)*(alpha1-alpha2+(1/3)*alpha3)*a(t)^3-(1/2)*alpha6*a(t)*cos(gamma(t))

(2)

"p(t):=a(t)*cos(gamma(t))"

proc (t) options operator, arrow, function_assign; a(t)*cos(gamma(t)) end proc

(3)

"q(t):=a(t)*sin(gamma(t))"

proc (t) options operator, arrow, function_assign; a(t)*sin(gamma(t)) end proc

(4)

diff(p(t), t)

(diff(a(t), t))*cos(gamma(t))-a(t)*(diff(gamma(t), t))*sin(gamma(t))

(5)

(-mu*a(t)-(1/4)*alpha6*a(t)*sin(gamma(t)))*cos(gamma(t))-a(t)*(2*sigma-(6*(1/8))*(alpha1-alpha2+(1/3)*alpha3)*a(t)^2-(1/2)*alpha6*cos(gamma(t)))*sin(gamma(t))

(-mu*a(t)-(1/4)*alpha6*a(t)*sin(gamma(t)))*cos(gamma(t))-a(t)*(2*sigma-(3/4)*(alpha1-alpha2+(1/3)*alpha3)*a(t)^2-(1/2)*alpha6*cos(gamma(t)))*sin(gamma(t))

(6)

diff(p(t), t)

2*t

(7)

``


subs.pdf

Download simplify.mw

 

 

First 749 750 751 752 753 754 755 Last Page 751 of 2224