## Sockets package issues - Serial interfacing ...

Asked by:

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

## Solving the following PDE...

Asked by:

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

## Corrupted file ...

Asked by:

Hi

I have a problem with my maple.

corrupt file

Is there anybody who can help me solve this problem.

Casper

## how i can Substituting......

Asked by:

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

 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 (7)
 >

Download simplify.mw

## How do I run coding with goto on Maple? ...

Asked by:

Hi, I'm using Maple 2018 and I tried to run coding from https://www.maplesoft.com/applications/view.aspx?sid=4194&view=html

however, it said : unable to parse. I figured out that the problem maybe is in the if loop. though it seems perfectly fine, but it has some goto commands that i cannot search on maple website. does this mean that the goto cannot be used here and should be replaced? if yes, then how?

i am still learning on how to use maple. any help would be much appreciated. thank you!

this is the coding for if loop:

label_7;

rv:=vector([p1(x1pt,x2pt),p2(x1pt,x2pt)]):

numgeval:=numgeval+1;

printf("%5d (%8.4f,%8.4f)",numIter,rv[1],rv[2]);

max:=n;

mg:=convert(sqrt(dotprod(rv,rv)),float);

printf("%12.4f",mg);

if(mg<tol or numIter>=max) then

goto(label_6);

else

numIter:=numIter+1;

fi;

v1:=x1pt+t*rv[1];

v2:=x2pt+t*rv[2];

newt:=evalf(subs({x1=v1,x2=v2},f1));

numfeval:=numfeval+1;

lam:=fsolve(diff(newt,t)=0,t,maxsols=1);

nv1:=evalf(subs({t=lam},v1));

nv2:=evalf(subs({t=lam},v2));

printf(" (%8.4f,%8.4f)%13.4f\n",x1pt,x2pt,lam);

x1pt:=nv1;

x2pt:=nv2;

goto(label_7);

label_6;

printf("\n\n-----------------------------------------");

printf("---------------------------------------------");

printf("\n\n Approximate Solution: ");

printf(" (%8.4f,%8.4f)\n",x1pt,x2pt);

Fvalue:=evalf(subs(x1=x1pt,x2=x2pt,f));

printf(" Maximum Functional Value: ");

printf("%21.4f",Fvalue);

printf("\n Number gradient evaluations:");

printf("%22d",numgeval);

printf("\n Number function evaluations:");

printf("%22d",numfeval);

printf("\n\n-----------------------------------------");

printf("---------------------------------------------");

end:

## Can someone explain what DESol mean as a solution ...

Asked by:

Hello,

I am fairly new to using the Maple software, so I apologize if my question is completely idiotic. Apologies, also, because I could not manage to enter my code as code. When I pressed the button it made the whole text as a code.

I run the following code to seek -if there are any- analytic solutions for the following differential equation.

odeplus := (r^2+L^2)^(5/2)*(diff(f(r), \$(r, 2)))+((15/4)*r*(r^2+L^2)^(1/2)+3*(r^2+L^2)^(5/2)/r)*(diff(f(r), r))+M^2*f(r)/(r^2+L^2)^(5/2)-((5/2)*((r^2+L^2)^(1/2))(l-1)+(55/64)*r^2/(r^2+L^2)^(3/2)+(r^2+L^2)^(5/2)*(l^2+3*l+3/2)/r^2)*f(r)+(((r^2+L^2)^(1/2))(5+(5/2)*l)+(5/8)*r^2/(r^2+L^2)^(3/2)-(r^2+L^2)^(5/2)*(3/2+l)/r^2)*f(r) = 0

and then I do

dsolve(odeplus, f(r))

The solutions that Maple returns is given in terms of DESol. Could anyone try and break it down for me? What is this telling me and if I can indeed from the output obtain analytic solutions? Is this some sort of operator acting on something?

Thank you in advance.

## Add IsPrime condition in a fonction SOLVE( EQNS,[x...

Asked by:

Hello ,

Iam new with maple, and iam trying to solve an equation using prime numbers. The equation is display as follow.

eqns := {G = x*q, 8695*y = 1341*q};

Iam trying to solve with the condition p and q are prime numbers

So is it posible to add condition prime in the solve fonction ? solve(eqns*[x, y]);

Thank you for your help

## Excessive "Recent" stacking...

Asked by:

I hate it when a user bumps every Question that they've ever asked over the past five years to the top of Recent (aka Active Conversations) stack! This is not the first time that this has happened. I rely heavily on that stack being truly the threads with the most-recent material. Can't something be done about this short of deleting the threads? That seems too drastic: While this user has occasionally been a pest in the past, they've also posted some decent material. There should be some way for a moderator to just reset the stack positions to what they were previously.

## simplify an equation ...

Asked by:

hi.

According to the fhgure attaceh how i can gain the equation (2-27) . I write the equation (2-26) in maple but I couldnot to gain that result.

If possible to reach equation via maple?

Thanks

 (1)

 (2)

Download diff.mw

## Strange Behavior of DynamicSystems[BodePlot]...

Asked by:

I tried the example in BodePlot help.

restart;
with(DynamicSystems):
sys := TransferFunction( 1/(s-10) ):
BodePlot(sys);

That works OK. But, if I invoke Syrup, the example no longer works.

restart;
with(Syrup);
with(DynamicSystems):
ckt := [V, Rsrc(50), C1(15e-9), L1(15e-6), C2(22e-9), L2(15e-6), C3(22e-9), L3(15e-6), C4(15e-9), 1, Rload(50)];

TF := subs(other, V=1, v[Rload]);
sys := TransferFunction(TF);

BodePlot(sys);
I get a message "not a valid plot structure".  OK, try the example, again.

sys := TransferFunction( 1/(s-10) ):
BodePlot(sys);
I also get the "not a valid plot structure" message.

What am I doing wrong?

## How to i solve this complex integration in Maple?...

Asked by:

Respected sir,

I have attached a file that contain integration . Sir I am unable to solve please give me some hint.

PLease note followings

1. I want to get final expression , want to get the value of constant from boundary condition.

2.I want to get the final expression (.......+c1)

3.I will put boundary condition to determine the value of c1

## Euler-method of forces on an airplane...

Asked by:

Can somebody please help me with this assignment?
Struggling a bit with the Euler-method from task d) and further.

## Fracsym package...

Asked by:

When i am running a code in maple worksheet , one error is shown by maple. My code and error (in bold) is below

Instructional workheet for the FracSym package
G. F. Jefferson and J. Carminati

Read in accompanying packages: ASP, DESOLVII and initialise using the with command:

read ASP v4.6.3.txt:

DESOLVII_V5R5 (March 2011)(c), by Dr. K. T. Vu, Dr. J. Carminati and Miss G.

Jefferson

The authors kindly request that this software be referenced, if it is used

in work eventuating in a publication, by citing the article:
K.T. Vu, G.F. Jefferson, J. Carminati, Finding generalised symmetries of

differential equations
using the MAPLE package DESOLVII,Comput. Phys. Commun. 183 (2012) 1044-1054.

-------------
ASP (November 2011), by Miss G. Jefferson and Dr. J. Carminati

The authors kindly request that this software be referenced, if it is used

in work eventuating in a publication, by citing the article:
G.F. Jefferson, J. Carminati, ASP: Automated Symbolic Computation of

Approximate Symmetries
of Differential Equations, Comput. Phys. Comm. 184 (2013) 1045-1063.

with(ASP);
[ApproximateSymmetry, applygenerator, commutator]
with(desolv);
[classify, comtab, defeqn, deteq_split, extgenerator, gendef, genvec,

icde_cons, liesolve, mod_eq, originalVar, pdesolv, reduceVar, reduceVargen,

symmetry, varchange]

Read in FracSym and initialise using the with command:
read FracSym.v1.16.txt;
FracSym (April 2013), by Miss G. Jefferson and Dr. J. Carminati

The authors kindly request that this software be referenced, if it is used

in work eventuating in a publication, by citing:
G.F. Jefferson, J. Carminati, FracSym: Automated symbolic computation of Lie

symmetries
of fractional differential equations, Comput. Phys. Comm. Submitted May 2013.

with(FracSym);
[Rfracdiff, TotalD, applyFracgen, evalTotalD, expandsum, fracDet, fracGen,

split]

BASIC OPERATORS

The Riemann-Liouville fractional derivatives is expressed in "inert" form using the FracSym routine Rfracdiff.
The explicit formula for the form of these fractional derivatives may be found in I. Podlubny, Fractional differential equations: An introduction to fractional derivatives, fractional differential equations, some methods of their solution and some of their applications, San Diego, 1999.)

Rfracdiff(u(x, t),t,alpha);
alpha
D[t     ](u(x, t))

If the fractional derivative is taken for a product, the generalised Leibnitz rule is used to express the result (the product operator used is &* and is non-commutative).
Rfracdiff(u(x, t)&*v(x,t),t,alpha);
infinity
-----
\
)                          (alpha - n)              n
/     binomial(alpha, n) D[t           ](u(x, t)) D[t ](v(x, t))
-----
n = 0
Rfracdiff(v(x, t)&*u(x,t),t,alpha);
infinity
-----
\
)                          (alpha - n)              n
/     binomial(alpha, n) D[t           ](v(x, t)) D[t ](u(x, t))
-----
n = 0

Fractional derivatives of integer order revert to the MAPLE diff routine.

Rfracdiff(u(x, t)&*v(x,t),t,2);
/ d  / d         \\             / d         \ / d         \
|--- |--- u(x, t)|| v(x, t) + 2 |--- u(x, t)| |--- v(x, t)|
\ dt \ dt        //             \ dt        / \ dt        /

/ d  / d         \\
+ u(x, t) |--- |--- v(x, t)||
\ dt \ dt        //

The FracSym rouine TotalD may also be used to find total derivatives. evalTotalD is then used to evaluate the result (in jet notation). For example,

TotalD(xi[x](x, y),x,2);
2
D[x ](xi[x](x, y))
evalTotalD([%],[y],[x]);
[     / d             \      2 / d  / d             \\
[y_xx |--- xi[x](x, y)| + y_x  |--- |--- xi[x](x, y)||
[     \ dy            /        \ dy \ dy            //

/ d  / d             \\       / d  / d             \\]
+ 2 |--- |--- xi[x](x, y)|| y_x + |--- |--- xi[x](x, y)||]
\ dy \ dx            //       \ dx \ dx            //]

EXAMPLE -  FINDING SYMMETRIES FOR A FRACTIONAL DE

Consider the fractional PDE from: R. Sahadevan, T. Bakkyaraj, Invariant analysis of time fractional generalized Burgers and Korteweg-de Vries equations, J. Math. Anal. Appl. 393 (2012) 341-347.

We use the Rfracdiff routine to express the
alpha
fractional derivative with respect to t:
fde1:=Rfracdiff(u(x, t),t,alpha) = (diff(u(x, t), x,x))+n*(u(x, t))^p*(diff(u(x, t),  x));
alpha             / d  / d         \\            p / d         \
D[t     ](u(x, t)) = |--- |--- u(x, t)|| + n u(x, t)  |--- u(x, t)|
\ dx \ dx        //              \ dx        /

sys1:=[Rfracdiff(u(x, t),t,alpha) = (diff(v(x, t), x)), Rfracdiff(v(x, t),t,alpha) = -u(x, t)*diff(u(x, t),x)];
[   alpha              d              alpha                      / d         \]
[D[t     ](u(x, t)) = --- v(x, t), D[t     ](v(x, t)) = -u(x, t) |--- u(x, t)|]
[                      dx                                        \ dx        /]

We use the the FracSym routine fracDet to find the determining equations for the symmetry for fde1.
NOTE: The fourth argument (some integer at least 1) corresponds to the number of terms to be "peeled off" from the sums which occur in the extended infintesimal function for the fractional derivative. A value of 2 provides a good balance between information for solution of determining equations and speed.

deteqs:=fracDet([sys1], [u, v],[x, t], 2, alpha=(0.1)..1);
Error, (in desolv/PickLHSDerivative) Cannot pick out the left hand side derivatives

Please suggest what problem it may be?

## can maple graph do "point to graph" to create feat...

Asked by:

would like to point to graph then it highlight graph with virtical line

and mark 1 in one of row in one of column in table like data

just like define feature manually for machine learning but using graph

## pdsolve does not work...

Asked by:

why pdsolve does not work?

Q and alpha are constant. I want to find f(x,y).

Thanks

 (1)

 (2)

Download PLATE.mw

 First 7 8 9 10 11 12 13 Last Page 9 of 241
﻿