Maple 18 Questions and Posts

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

We have the system with one discrete variable along x-axis (i.e. 'i' is discrete in the attached file) and other variable 't' is continuous. But maple return error.

CD.mw

Is there any command that allows me to extract the input number of any function f?

exmple:

f: (x1; x2;...;xn)-> y

command (f)=n

thanks for the help!

I need to define a simple recursive algorithm (i'm not a programmer) such that:

xi=xi-1+b 

 

with i=0,1,...,n; 

and with all xi's elements of a Set A

how can i achieve this?

hallo every body 

please can you help me 

how do i solve this differential linear system with respect lambda is positive number  

i use maple 18

Let the differential system with $\lambda>0$

\begin{equation}
\begin{array}{ccc}
\dot{x}=y(t)\\
\dot{y}=z(t)\\
\dot{z}=-\lambda y(t)-h(t)
\end{array}
\end{equation}

prob.pdf

Why maple return trivial solution after integration (see (11) in the attached .mw))? The result should be some non-trivial solution.

solutionsg.mw

I want to find the relation between w and k by putting Determinant(A) is equal to zero. Since the determinant is large (A is a matrix of size 6x6) and I am not able to find the simplified expression. How to do that?

 

determinant.mw

Hi,

Please can someone help me with a sample code for bifurcation? You can use parameter values for the parameters. I'm using maple 18. Below is my model:

restart:

f__1 := Delta -(psi + mu)*S(t);

Delta-(psi+mu)*S(t)

(1)

f__2 := psi*S(t) -(delta + mu)*E(t);

psi*S(t)-(delta+mu)*E(t)

(2)

f__3 := Delta*E(t) -(gamma+gamma__1 + mu)*X(t);

Delta*E(t)-(gamma+gamma__1+mu)*X(t)

(3)

f__4 := gamma__1*X(t)-(eta + xi + mu)*H(t);

gamma__1*X(t)-(eta+xi+mu)*H(t)

(4)

f__5 := xi*H(t) - mu*R(t);

xi*H(t)-mu*R(t)

(5)

f__6 := gamma*X(t)-eta*H(t) - d*D(t);

gamma*X(t)-eta*H(t)-d*D(t)

(6)

f__7 := b*D(t) - b*B(t);

b*D(t)-b*B(t)

(7)

f__8 := phi__p + sigma*X(t)+eta__1*H(t) +d__1*D(t)+ b__1*B(t) - alpha*P(t);

phi__p+sigma*X(t)+eta__1*H(t)+d__1*D(t)+b__1*B(t)-alpha*P(t)

(8)

 

NULL

Download Bifurcation.mw

Non-Linear.mw

Hi, I have here a interesting non-linear system.

If I attempt to solve it using some specific form of the non-linear equations (form X*Y=Z) of the system, Maple (Verison 18) finds a solution.

But, if I replace some of them by some other forms (like form Y=Z/X), fsolve fails.

I usually use the non-quotient form. But is there any way to guide or configure fsolve to reach a solution?
I set up some of the regular options: placing a seed close to the solution, indicating intervals of possible solutions; but none of that works if I do not set up the non-quotient form of the equations. In some cases, fsolve does not reach a solution at all, no matter the form of the equations.

In the file, the equations that are causing the isssue are the last 3, those who start with the variable f1,f2 and f3.
I ran the system twice with both cases: non-quotient form and quotient form.

Thanks for your attention! 

Given a set (list) of PDE, is there a way to search all possible solution sets? For instance, pdsolve will output the solution

{_eta[0](t, x) = 0, _xi[t](t, x, u) = _C1, _xi[x](t, x, u) = _C2, eta[1](t, x) = 0}

for the list of PDEs below. But I am aware that there is another solution different from the above, is there way to seek these other solutions?

 

[alpha*u^2*(diff(eta[1](t, x), t))+alpha*u*(diff(_eta[0](t, x), t))-u*(diff(eta[1](t, x), t))-u*(diff(eta[1](t, x), x))+u*(diff(eta[1](t, x), x, x, t))-(diff(_eta[0](t, x), x))-(diff(_eta[0](t, x), t))+diff(_eta[0](t, x), x, x, t), -(diff(_xi[x](t, x, u), u, u, u)), -(diff(_xi[t](t, x, u), x))-(diff(_xi[t](t, x, u), t)), -(diff(_xi[t](t, x, u), x, x)), -2*(diff(_xi[t](t, x, u), x, x))-(diff(_xi[x](t, x, u), x, x))+2*(diff(eta[1](t, x), x)), diff(eta[1](t, x), t)-2*(diff(_xi[x](t, x, u), x, x)), -(diff(_xi[x](t, x, u), t))*alpha*u-(diff(_xi[x](t, x, u), x, x, x))-(diff(_xi[t](t, x, u), x))+2*(diff(eta[1](t, x), x, t)), -(diff(_xi[x](t, x, u), x)), -(diff(_xi[t](t, x, u), t))*alpha*u+(diff(_xi[t](t, x, u), x))*alpha*u+(diff(_xi[x](t, x, u), x))*alpha*u+(diff(_xi[x](t, x, u), t))*alpha*u+eta[1](t, x)*alpha*u+alpha*_eta[0](t, x)+diff(_xi[t](t, x, u), t)-(diff(_xi[t](t, x, u), x, x, x))-(diff(_xi[x](t, x, u), x))-(diff(_xi[x](t, x, u), t))+diff(eta[1](t, x), x, x), -2*(diff(_xi[x](t, x, u), u)), -2*(diff(_xi[t](t, x, u), x, u)), -2*(diff(_xi[t](t, x, u), u))-(diff(_xi[x](t, x, u), u)), -(diff(_xi[t](t, x, u), u)), -(diff(_xi[t](t, x, u), u)), -5*(diff(_xi[x](t, x, u), u, u)), -(diff(_xi[t](t, x, u), u, u)), -3*(diff(_xi[t](t, x, u), x, u)), -2*(diff(_xi[x](t, x, u), u)), -(diff(_xi[t](t, x, u), u)), -(diff(_xi[t](t, x, u), u)), -2*(diff(_xi[t](t, x, u), u, u))-(diff(_xi[x](t, x, u), u, u)), -(diff(_xi[t](t, x, u), u, u, u)), -3*(diff(_xi[t](t, x, u), u, u)), -3*(diff(_xi[t](t, x, u), x, u, u)), (diff(_xi[t](t, x, u), u))*alpha*u-(diff(_xi[x](t, x, u), u))-3*(diff(_xi[t](t, x, u), x, x, u)), -2*(diff(_xi[x](t, x, u), x, u))-4*(diff(_xi[t](t, x, u), x, u)), -(diff(_xi[t](t, x, u), u))*alpha*u+(diff(_xi[x](t, x, u), u))*alpha*u+diff(_xi[t](t, x, u), u)-(diff(_xi[x](t, x, u), u)), -3*(diff(_xi[x](t, x, u), x, x, u))-(diff(_xi[t](t, x, u), u)), -7*(diff(_xi[x](t, x, u), x, u)), -3*(diff(_xi[x](t, x, u), x, u, u))]

Could you Please Help me,the Maple code for Plot this equations in any Numerical Method

restart

with(PDETools, declare):

with(LinearAlgebra):

with(PolynomialIdeals):

with(plots):

declare(f(x));

f(x)*`will now be displayed as`*f

(1)

eq1 := diff(f(x), x, x, x)+(1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[fl])*(eta*`cosω`+f(x)*`sinω`)*(diff(f(x), x, x))+(1-phi)^2.5*M*sin^2*alpha*(1-(diff(f(x), x)))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[fl])*lambda[T]*theta = 0;

diff(diff(diff(f(x), x), x), x)+(1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[fl])*(eta*`cosω`+f(x)*`sinω`)*(diff(diff(f(x), x), x))+(1-phi)^2.5*M*sin^2*alpha*(1-(diff(f(x), x)))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[fl])*lambda[T]*theta = 0

(2)

eq2 := K[nf]*(diff(theta(x), x, x))/K[f]+(1/2)*Pr*(eta*`cosω`+f*`sinω`)*(diff(theta(x), x)) = 0;

K[nf]*(diff(diff(theta(x), x), x))/K[f]+(1/2)*Pr*(eta*`cosω`+f*`sinω`)*(diff(theta(x), x)) = 0

(3)

bcs := F(0) = 0, F(1) = 0, F(10) = 1, Theta(0) = 1, Theta(10) = 0;

F(0) = 0, F(1) = 0, F(10) = 1, Theta(0) = 1, Theta(10) = 0

(4)

a1 := [M = 1, alpha = 0, phi = 0.5e-1, `cosω` = 1, `sinω` = 1, sin(alpha) = 0, lambda[T] = 0, Pr = 6.2, rho[s] = 5200, rho[fl] = 997.1, `ρβ`[fl] = 20939.1, K[nf] = .6842, eta = 0];

[M = 1, alpha = 0, phi = 0.5e-1, `cosω` = 1, `sinω` = 1, sin(alpha) = 0, lambda[T] = 0, Pr = 6.2, rho[s] = 5200, rho[fl] = 997.1, `ρβ`[fl] = 20939.1, K[nf] = .6842, eta = 0]

 

[Pr = 6.2, M = 1, phi = 0.5e-1, `cosω` = 1, `sinω` = 1, sin(alpha) = 0, lambda[T] = 0, `ρβ`[fl] = 20939.1, K[nf] = .6842, K[f] = .613, eta = 0]

(5)

b1 := subs(a1, eq1);

diff(diff(diff(f(x), x), x), x)+.5325197465*f(x)*(diff(diff(f(x), x), x)) = 0

(6)

c1 := subs(a2, eq2);

1.116150082*(diff(diff(theta(x), x), x))+3.100000000*f*(diff(theta(x), x)) = 0

(7)

d1 := dsolve({b1, bcs}, {bcs, c1}, numeric); d1(0)

Error, (in dsolve/numeric/process_input) invalid argument: {1.116150082*(diff(diff(theta(x), x), x))+3.100000000*f*(diff(theta(x), x)) = 0, F(0) = 0, F(1) = 0, F(10) = 1, Theta(0) = 1, Theta(10) = 0}

 

d1(0)

(8)

``

Download  Numerical Method.mw

Why solutions (5) and (6) are different? The solution (5) is obtained by putting f=0 of the series, while (6) is the result by taking limit at f->0. 

check.mw

hallo evert body please how i do calculate this integral

in maple 18

ttp.pdf

int_{0}^{2\pi}(\cos^{i}(\theta)-\cos^{i+2}(\theta))d\theta

Please I will like to know if there exists any command to isolate the known regular mathematical functions from an expression. For example, given the equation 

I need a means of isolating (ex & sin(v)) from the above expression

Thanks

Please, how do plot the contour plot of im(qq) vs Re(qq) from the list below: 

qq= [2.106333379+.6286420119*I, 2.106333379-.6286420119*I, 4.654463885, 7.843624703, 10.99193295, 14.13546782, 17.27782732, 20.41978346, 23.56157073, 26.70327712, 29.84494078, 32.98658013, 36.12820481, 39.26982019, 42.41142944, 45.55303453, 48.69463669, 51.83623675, 54.97783528, 58.11943264, 61.26102914]

How can we export image with minimum size? I plot the solution and then export the image as an .eps file. But when I try to generate the generate the pdf file on Latex, the file size is too large. But when I plot the solutions on Matlab and export the images as .eps file, then Latex generate pdf file with samll size. Why Maple generate images of large sizes?

solution.mw

1 2 3 4 5 6 7 Last Page 1 of 79