Maple 11 Questions and Posts Maple 11 Questions and Posts Feed

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

In point of fact, most people are stupid when it is on par with Follinique. We'll look at the Follinique strategies you are using. Follinique has had lasting success. I was just thinking touching on Follinique and there's no moss growing on Follinique. You need a quality Follinique. I came across a Follinique that legalizes an aura for a Follinique. Knock 'em dead, tiger! Clearly, happily, this is not the complete story. As always, here is my uncomplicated solution to that problem. I'm willing to show you some information into what I've learned relevant to Follinique since they see that Follinique has made them happy. I know it is difficult to tell us all something that points out Follinique so poorly. This will be good for something. You just made my day. Follinique Reviews @ http://www.revommerce.com/follinique/

what should i do to solve this ?!

exit:restart;

grtw();

`GRTensorII Version 1.79 (R4)`

`6 February 2001`

`Developed by Peter Musgrave, Denis Pollney and Kayll Lake`

`Copyright 1994-2001 by the authors.`

`Latest version available from: http://grtensor.phy.queensu.ca/`

`E:/Gravitation/Grtii(6)/Metrics`

(1)

makeg(Einstein5);

 

 

Makeg 2.0: GRTensor metric/basis entry utility

 

To quit makeg, type 'exit' at any prompt.

 

Do you wish to enter a 1) metric [g(dn,dn)],

                       2) line element [ds],

                       3) non-holonomic basis [e(1)...e(n)], or

                       4) NP tetrad [l,n,m,mbar]?

 

 
makeg>

2;

Enter coordinates as a LIST (eg. [t,r,theta,phi]):

 
makeg>

[t,r,theta,phi,psi];

Enter the line element using d[coord] to indicate differentials.

(for example,  r^2*(d[theta]^2 + sin(theta)^2*d[phi]^2)

[Type 'exit' to quit makeg]

 ds^2 =

 
makeg>

f(r)*d[t]^2+d[r]^2/f(r)+r^2*(d[theta]^2+sin(theta)^2*(d[phi]^2+sin(phi)^2*d[psi]^2));

 

If there are any complex valued coordinates, constants or functions

for this spacetime, please enter them as a SET ( eg. { z, psi } ).

 

Complex quantities [default={}]:

 
makeg>

{};

`The values you have entered are:`

Coordinates = [t, r, theta, phi, psi]

`Metric:`

"g[a] [b]=[[[f(r),0,0,0,0],[0,(1)/(f(r)),0,0,0],[0,0,r^2,0,0],[0,0,0,r^2 (sin(theta))^2,0],[0,0,0,0,r^2 (sin(theta))^2 (sin(phi))^2]]]"

You may choose to 0) Use the metric WITHOUT saving it,

                  1) Save the metric as it is,

                  2) Correct an element of the metric,

                  3) Re-enter the metric,

                  4) Add/change constraint equations,

                  5) Add a text description, or

                  6) Abandon this metric and return to Maple.

 

(2)
makeg>

1;

Information written to: `E:/Gravitation/Grtii(6)/Metrics/Einstein5.mpl`

Do you wish to use this spacetime in the current session?

(1=yes [default], other=no):

 
makeg>

1;

Initializing: Einstein5

`Default spacetime` = Einstein5

`For the Einstein5 spacetime:`

Coordinates

x(up)

`x `^a = (Vector[row](5, {(1) = t, (2) = r, (3) = theta, (4) = phi, (5) = psi}))

`Line element`

` ds`^2 = f(r)*` d`*t^`2 `+` d`*r^`2 `/f(r)+r^2*` d`*theta^`2 `+r^2*sin(theta)^2*` d`*phi^`2 `+r^2*sin(theta)^2*sin(phi)^2*` d`*psi^`2 `

makeg() completed.

(3)

grdef(`G2{e f}:=g{e f}*(R{a b c d}*R{^c ^d ^a ^b}-4*R{^a ^b}*R{a b}+Ricciscalar^2)+4*(R{f c a b}*R{^a ^b ^c e}+2*R{f a e b}*R{^b ^a}+2*R{f c}*R{^c e}-Ricciscalar*R{f e})`);

Created definition for R(up,dn)

Created definition for R(up,up)

Created definition for R(up,up,up,up)

Created definition for R(up,up,up,dn)

Created definition for G2(dn,dn)

 

grdef(`G1{e f}:=G{e f}+Lambda*g{e f}`);

Created definition for G1(dn,dn)

 

grcalc(G1(dn,dn));

`CPU Time ` = .156

(4)

grcalc(G2(dn,dn));

`CPU Time ` = .141

(5)

grdef(`A{a}:=[h(r),0,0,0,h(r)*N*6*sin(theta2)^2*sin(theta1)^2,h(r)*N*6*sin(theta2)^2]`);

Error, (in grF_grdef) Dimension of default metric and number of components are not equal

 

 

 

Download c.mws

How do i make 3d kink kink collition soliton  animation  graph in maple.

Hi 

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 http://www.maplesoft.com/support/help/Maple/view.aspx?path=examples/pdsolve_boundaryconditions , 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]];
pdsolve(sys[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.

 

Hello

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?

Hello,

 

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.

 

Thanks.

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!

maple.mws


restart;
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 :=
evalf(dsolve({
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;


od:
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

1 2 3 4 5 Page 1 of 5