in harmonic ocillator version, why it has error after it dsolve?

with(Physics[Vectors]):
with(SumTools):
Setup(mathematicalnotation = true):

Qs is Qss*cos(m*t+phi);

H := 1/(8*Pi*c^2)*Summation(Commutator(diff(Qs*cos(w*t+phi),t)*diff(conjugate(Qs*cos(w*t+phi)),t)+w^2*Qs*cos(w*t+phi)*conjugate(Qs*cos(w*t+phi))),ks=1..3)-1/c*(Summation(conjugate(Qs*cos(w*t+phi)),ks=1..3));

or

H := 1/(8*Pi*c^2)*Summation(Commutator(diff(Qs,t)*diff(conjugate(Qs),t)+w^2*Qs*conjugate(Qs)),ks=1..3);

#Qs = p, Ps := Pss*sin(m*t+phi) Ps = qs

i use

H := 1/(8*Pi*c^2)*Summation(Commutator(diff(Qs*cos(w*t+phi),t)*diff(conjugate(Qs*cos(w*t+phi)),t)+w^2*Qs*cos(w*t+phi)*conjugate(Qs*cos(w*t+phi))),ks=1..3)-1/c*(Summation(conjugate(Qs*cos(w*t+phi)),ks=1..3));

 HJ := subs(Qs= diff(f(q,P,t), q), H);

H:=subs( f(q,P,t) = f1(q) + f2(t), HJ);
sol:=dsolve({rhs(H)=E,lhs(H)=E});
S:=rhs(sol[1][1]+sol[1][2]);
p:=diff(S,q);
Q:=diff(S,E);

Hello

when i am trying to isolate for Rev in the following expression

the program doesnt respond.do you know why?

I have a vector of dimension n with each component being an equation of a linear system.

Can maple convert this Vector to a Matrix-Vector form with the matrix being constant coefficients?

Let B is a list of polynomial conditions such that  are none zero. Consider one polynomial f. How can I detect that f is none zero w.r.t. B? For example if B=[a-1,b+2,b-c,ac-1] and f=a^2c-ac-a+1. From B we can conclude that a<>1 and b<>2 and b<>c and ac<>1. How can I deduce that f<>0 w.r.t. B automatically?

Hi all,

I have been trying to plot Riemann Zeta Function on certain range, but it takes a lot of time(several minutes ,maybe hours).These are the basic things.

How to speed up ploting Zeta?

The same problem with Optimization.

Zeta-Problems.mw

Mathematica don't have this such problems.

Mathematica ploting Zeta with 0.125 second and optimizing 3.42 second.Maple few hours maybe......

How would I find inflection points? I believe it would the same as if I was finding critical points for f '.

Function is:(7-x)*sin(x^2-7)

So would it be:

a:=fsolve(f2,-2.4..2.4);
b:=fsolve(f2,-2.4..-0.06);
c:=fsolve(f2,-0.06..2.4);
d:=fsolve(f2,1.49..3);
e:=fsolve(f2,-3 ..-1.52);

Plz and thanks!

Hi all,

I am a new use of Maple. I have a expression :

2.560000000*10^(-30)*k1*(2.309486127*10^38*n-1.154743064*10^38-1.186994552*10^37*k2^2*n^2+2.373989104*10^37*k2^2*n+8.541613702*10^37*k2*n^2-1.154743064*10^38*n^2+7.119519043*10^37*n^2*k1^2+6.495962587*10^37*k1*n^2-1.293058266*10^38*k2*n^2*k1-1.186994552*10^37*k2^2-1.708322740*10^38*k2*n+1.293058266*10^38*k2*n*k1+8.541613702*10^37*k2-3.529088273*10^37*k1-2.966874313*10^37*k1*n)/(1.620847396*10^9*n-1.620847396*10^9)

I want to simplify it. Obviously, it can be simplifed by simplifying the coeffecient of each terms. 

When I use 'simplify' command in Maple 18, I couldn't got what I want. 

Any suggetstions?

Thanks in advance.

Hi, 

  I have the following input

***

restart;
with( Statistics ):

a:=2;c:=0.3;
g:= exp(-a*x) + c*a*exp(-a*x);
#f := x -> piecewise( x < 0, 0, x>0, g );
 f :=x -> piecewise( x < 0, 0, x>0, exp(-a*x) + c*a*exp(-a*x));

norm_factor:=int( f(x), x=0..infinity );
print(norm_factor);

randomize():
F := Distribution( PDF = 1/norm_factor*f ):
X := RandomVariable( F ):

N := 20;
S := convert( Sample(X,N), list );

print(`cc`,S[1]);

***

The code works. However, if I comment out 

 f :=x -> piecewise( x < 0, 0, x>0, exp(-a*x) + c*a*exp(-a*x));

 , then use

f := x -> piecewise( x < 0, 0, x>0, g );

i.e.

f := x -> piecewise( x < 0, 0, x>0, g );
 #f :=x -> piecewise( x < 0, 0, x>0, exp(-a*x) + c*a*exp(-a*x));

It is said "

Error, (in Statistics:-Sample) unable to construct the envelopes for _R, try to specify the initial range"

The norm_factors are actually the same for both inputs. What is the reason for the error message?  Suppose I still want to use something like

f := x -> piecewise( x < 0, 0, x>0, g );

,how to fix the problem?

Thank you very much

I'm drawing a complete blank on this. 

The probability of a branch breaking in a particular tree as someone climbs is 2%.  After the person has climbed over 40 branches, did any of them break?  I'm trying to set that up as a simulation in Maple but I'm not sure about how to do it. 

I guess if I generate a random number and it falls under .02 it breaks, is that how I determine that?

randomize():

n := 0;
for i to 40 do
ran := evalf(rand()*10^(-12));
if ran < 0.2e-1 then n := n+1: end if:
end do:
print("Branch broke", n, "times")

I'm not sure I've got the probability right.

Shading is not valid as an option in the command defining a plottools object e.g. a parallelepiped, and yet shading can be selected in the actual display of the object. Is there a way to code for shading when creating the object?

Hi,

Ive been trying to solve critical points for maple but i keep getting this werid equation:

(1/7)*RootOf(tan(_Z)^2+4*_Z*tan(_Z)+4*_Z^2+28*tan(_Z)-140*_Z-1176)+1+(1/14)*tan(RootOf(tan(_Z)^2+4*_Z*tan(_Z)+4*_Z^2+28*tan(_Z)-140*_Z-1176)).

What I have typed

f:=x->(7-x)*sin(x^2-7);

f1:=D(f);

When I try putting:

solve(f1(x)=0);

I get that werid equation

Hello

Can I make a button with the code so it will execute the entire worksheet?

I have tried without sucess to run the

Equation manipulator.

I get the following error:

"

Error, (in Maplets:-Tools:-CleanElements) invalid input: op expects 1 or 2 arguments, but received 0

"

Do I have to load some package for it to run?

I have searched for a solution but can not find one.

Getting desperate.. anyone?

TIA

LWc

I want write loop for this code:

V := 24;

eq := {Eq2, Eq3};

ans := fsolve(eq, {a1 = 0, a3 = 0});

r := 0;

W1 := rhs(ans[1])*phi11(r)+rhs(ans[2])*phi21(r)

V and W1 is Variable, I change V and solve eq and determinen a1 and a3 after i calculated W1

I want write Loop for this,

can you help me????????/ 

%piecewise(``, x+y-2*z = 1, ``, 2*x+y-3*z = 5, ``, -x+y+z = 1); What does % here do?

I can find no help page to explain the functionality of the percent sign in this command and removing it produces an error. Please explain this.