Maple 11 Questions and Posts

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

With the following command I can plot two spheres and plot them.

f1 := x^2+y^2+z^2 = 1

f2 := x+y+z = 1



S1 := implicitplot3d(f1, x = -1 .. 1, y = -1 .. 1, z = -1 .. 1, style = patchnogrid, color = blue, scaling = constrained, axes = boxed)

S2 := implicitplot3d(f2, x = -1 .. 1, y = -1 .. 1, z = -1 .. 1, style = patchnogrid, color = gold, scaling = constrained, axes = boxed)


My questions are:

1- How can I display (highlight) the circle which is the intersection between these two sphere on the same figure?

2- How can I find the equation of this circle?

Thank you.

ifactor(8650368) produces

                                                        2503  *  2^7 * 3^3


how can I produce the output with factors in order from small to large and printed f rom left


2^7 * 3^3 * 2503



I have 3 equations

eq1 := vs = (Rs+Z)*i1+Z*i3

eq2 := A*vi = Z*i1+(Z+Rf+ro)*i3

eq3:= vo = (Rf+Z)*i3+Z*i1


and I want to solve for vo/vs . How to do that ?

the expected solution is 

Please help me to solve this integration

restart; with(LinearAlgebra); int(exp(-(ln(z/(snr*B^2))+4*sigma^2)^2/(32*sigma^2))*eta^2*(y/z)^((1/2)*eta^2-1)/(z*sqrt(32*Pi*sigma^2)*(2*sqrt(y*z))*(2*A[o]^(eta^2))), z, z = y/A[o]^2 .. infinity);




restart; with(LinearAlgebra); int(exp(-(ln(z/(snr*B^2))+4*sigma^2)^2/(32*sigma^2))*eta^2*(y/z)^((1/2)*eta^2-1)/(z*sqrt(32*Pi*sigma^2)*(2*sqrt(y*z))*(2*A[o]^(eta^2))), z, z = y/A[o]^2 .. infinity)

I'm trying to execute the program, which can be found here , but it does not work. I copied exactly what is written there:

restart; with(PDEtools):
U := diff_table(u(x,t)):
pde[1] := U[t]+c*U[x]=-lambda*U[];
bc[1] := eval(U[], t=0) = phi(x);
sys[1] := [pde[1], bc[1]];

But after last command it just sais that

Error, (in pdsolve/sys/info) found functions with same name but depending on different arguments in the given DE system: [u(x,t), u(x,0)]

What's wrong?

Dear Community Members,


We have problem with calculation in Maple v11 and v18. when we make a calculation by using maple v11 and v18, we was not able to get the solution as you see enclosed. when we clicked to "enter + ; ", programme does not run.



I have a loop with the do structure but there is an error in the loop .

how can I continue the loop by error or disregard it?



could you help me solve this error ? I don't understand what it means.


> eq3:=diff(x(t),t,t)+Gamma*diff(x(t),t)+omega[0]^2*(x(t)-(diff(x(t),t,t)+Gamma*diff(x(t),t)+omega[0]^2*x(t)+omega[0]^2*X[0])/omega[0]^2) = -omega[0]^2*X[0]:
> dsolve(eq3);
Warning, it is required that the numerator of the given ODE depends on the highest derivative. Returning NULL.



I have a great problem with this integral and Maple gives two answers completely different:


int(x^-5/3*cos((x-1)*h), x = 0..infinity)

so I get two different results :


-(27/8)*h^2+3/2+(27/8)*h^(7/6)*LommelS2(11/6, 1/2, h)


or this:


-(27/8)*h^2+3/2+(27/8)*h^(7/6)*LommelS1(11/6, 1/2, h)

In the first integral A get Lommels2 and If I get the Integral by using Taylor of cos((x-1)*h) and after that I resum I get Lommels1.


Thank you.




Greetings everyone. I tried to get the roots of this polynomial by using allvalues command. But I got the roots in indeces instead of the value of the roots.

what should i do?

Hello! I wrote a program in Maple but it doesn't work. It crashes with the following error: "Error, (in r_nach_1[3]) too many levels of recursion". I am a newbie in Maple and don't know how to solve this problem. A part of my program that causes the error is given below. Many thanks!


t1F := 20: A11 := 2.5*10^(-3): c1 := 4*10^2: ss1 := 0.5*10^(-2):
pI1 := 0.4*10^(-1): r1max := 0.6*10^(-2): procent:=0.2:
A12 := procent*A11: B1:=-1: alpha:=0.000001:

beta1:=0.5: r_nach_1[1]:=0: r_nach_1[2]:=0.2*10^(-2):

iter_psi_c1_1:=proc(t) 0 end proc: iter_psi_c2_1:=proc(t) 0 end proc:
s1[1]:=proc(t) 0 end proc:

for i from 1 to 2 do

_p11 := evalf(dsolve({
diff(p1(t), t) = A11*(p1(t)-'s1'[i](t))+B1*(r_nach_1[i+1](t)-r_nach_1[i](t))+A12*(p1(t)-'s1'[i](t))^2+ss1,
p1(0) = pI1}, numeric,method = dverk78, abserr = 1.*10^(-8), relerr = 1.*10^(-8),optimize,output = listprocedure, known=[s1[i],iter_psi_c1_1,iter_psi_c2_1])):

p1F := rhs(_p11(t1F)[2]):
s1[i+1] := subs(_p11, p1(t)):

q11 := evalf(dsolve({diff(q1(t), t) = c1*r_nach_1[i+1](t), q1(0) = 0},numeric,range=0..t1F,known=[iter_psi_c1_1,iter_psi_c2_1])):

q1F := rhs(q11(t1F)[2]):

F[i]:= beta1*q1F;

_psi_c1_1 :=
diff(psi_c1_1(t), t) = A11+2*A12*('s1'[i+1](t)-'s1'[i](t))*psi_c1_1(t), psi_c1_1(t1F) =  -1+beta1},numeric,method = dverk78, abserr = 1.*10^(-8), relerr = 1.*10^(-8),optimize,known=[s1[i],s1[i+1]],output=listprocedure)):

_psi_c2_1 :=evalf(dsolve({
diff(psi_c2_1(t), t) = 0, psi_c2_1(t1F) =  -beta1},numeric,output=listprocedure)):
iter_psi_c1_1 := op([2,2],_psi_c1_1):
iter_psi_c2_1 := op([2,2],_psi_c2_1):

r_nach_1[i+2] := t-> r_nach_1[i+1](t) +  alpha*(B1*(iter_psi_c1_1(t)) + (c1*iter_psi_c2_1(t))):

_p11 := evalf(dsolve({
diff(p1(t), t) = A11*(p1(t)-'s1'[i+1](t))+B1*(r_nach_1[i+2](t)-r_nach_1[i+1](t))+A12*(p1(t)-'s1'[i+1](t))^2+ss1,
p1(0) = pI1}, numeric,method = dverk78, abserr = 1.*10^(-8), relerr = 1.*10^(-8),optimize, output = listprocedure, known=[s1[i+1],iter_psi_c1_1,iter_psi_c2_1])):

p1F := rhs(_p11(t1F)[2]):

s1[i+2] := subs(_p11, p1(t)):

q11 := evalf(dsolve({diff(q1(t), t) = c1*r_nach_1[i+2](t), q1(0) = 0},numeric,range=0..t1F,known=[iter_psi_c1_1,iter_psi_c2_1])):

q1F := rhs(q11(t1F)[2]):

F[i+1]:= beta1*q1F;

Error, (in r_nach_1[3]) too many levels of recursion

1. nonlinear ODE with parameter "epsilon"

(x^n +epsilon*y(x))dy/dx + n*x^(n-1) * y(x) =m*x^(m-1) ; y(1)=b>1

where n=2,3,4,.. and m=0,1,2,3,...


2. Duffing equation with parameter "epsilon"

d^2 y(x)/dx^2 + y(x) + epsilon*y(x)^3=0 ; y(0)=A ; y'(0)=0

Hi, im new with maple

so I am still a bit confused with how to use it..

i have an nonlinear ODE with parameter "b"

(x+b*y(x))dy/dx + y(x) =0 ; y(1)=1

if i want b=0.1 and x=0, 0.1, 0.2, ..,1

how do i

Useing the following procedure I'd like to collect a set of  roots in a list of lists, so they can be manipulated and presented in various plot options. Some of which could be very interesting (to me only perhaps).

One option is to create a loop which produces values of newton roots x0 for a given function. For example:              f:=x->(x^7)-5

for x from -1 to 1 by 0.05 do

To highlight the problem I am haveing let us observe the following example. If we solve the following numerically,

d2/dx2 (s(x)) +s(x) -1=0,

ic: s(0)=1/2,s'(0)=0

How would I go about createing a plot of (s(x),x) for  x a member of [0,2*Pi) (useing polar coordinates clearly)?

The probelm I am haveing is plotting for x a member of [0,2*Pi). I know how to do this for x a member of [0,2*Pi] and (0,2*PI), but not  [0,2*Pi).

1 2 3 4 5 Page 1 of 5