The version 16 of Maple does not know the function "pow", it is replaced by "power". Nevertheless, the codegen[optimize] procedure still generates it!  Moreower, it's not optimal even to write poser(..., 1/2)  as, most probably, the function sqrt must run faster. 

  Here's an example: (the function being optimized just contains some nested square roots)

> iang10 := [codegen[optimize](Re(int11), tryhard)];
[       2        t17                                   

Hi, experts

I have an array which is empty, say

f:= []:

is there away to check if it's empty? I want to do something like in Matlab:

if isempty(f)

......

end

Many thanks

Kyle

Hi together,

i wanted to assume the sequence term A(n) of the first Perrin-Pseude-Prime n = 271441.

It´s not that difficult, but i do not understand the Maple-output, and i´ve never found an explenation for it.

the input was: > MatrixVectorMultiply(M^90480, v);  --> its the Matrix formula for the recurrence.

the output gives me a vector (3x1) with the following entry in the second line (witch is the line i want to know):

Hello colleagues! :)

I have some problems after create two objects, with array row in object definition...

The result - two same object...

Use the Gram-Schmidt process to transform the basis
 
           {Vector(3, {(1) = 0, (2) = 1, (3) = 1}), 

             Vector(3, {(1) = 1, (2) = 1, (3) = 1}), 

             Vector(3, {(1) = 1, (2) = 2, (3) = 3})}

for the Euclidean space R3 into an orthonormal basis for R3

                               
Show 3D result.

I have a multivariate polynomial equation, in that somehow I know the coeffcients, using this information, I want to extract the variables. This will be the opposite of coeffs function.

for e.g. I have 3*x3 + 5*x4

Given 3 and 5, I want to extract x3 and x4.

Thanks in advance.

Satya

I believe this is a result of my ignorance rather then the flaws of the system. Can anyone teach me how to do this question?

Given 

>sum(1/n, n = 1 .. N)

I want to find the value of N such that sum(1/n, n = 1 .. N) is greater then 20.12.

>solve(sum(1/n, n = 1 .. N) > 20.12, N)

Warning, solutions may have been lost

Help please! Thanks.

Hi, experts,

I have an equations:

How to solve a, and b? It's easy to know that a = 2, b =5 though.

Many thanks! 

In the name of God

Any one knows how I can summerize this:

v1:=Vector(3):v2:=Vector(3):v3:=Vector(3):v4:=vector(3):v5:=Vector(3):

these 5 are independent.

Maybe some help page of maple helps me.

Hossayni

How can i solve for G(p,q,r,t) in the following pde

diff(G(p, q, r, t), t) = {lambda*(p-1)+gamma*(r-1)}*G(p, q, r, t)+delta*(1-p)*(diff(G(p, q, r, t), p))+kappa*(N*r-q)*(diff(G(p, q, r, t), q))+mu*(1-r)*(diff(G(p, q, r, t), r))+beta*(q*e^(-`ρτ`)-p*r)*(diff(G(p, q, r, t), r, p))

Dear All

I trying to solve the following in MAPLE, the code works fine in MATLAB but I am not sure why it doesn't give me values in MAPLE

n:=3: alpha_p:=1: p:=2: mu0:=4*Pi*1e-7: Br:=1.12:

A1:=unapply(sin((n*p+1)*alpha_p*Pi/(2*p))/((n*p+1)*alpha_p*Pi/(2*p)),n):

A2:=unapply(sin((n*p-1)*alpha_p*Pi/(2*p))/((n*p-1)*alpha_p*Pi/(2*p)),n):

M1:=unapply((Br/mu0)*alpha_p*(A1(n)+A2(n)),n):

M2:=unapply((Br/mu0)*alpha_p*(A1(n)-A2(n)),n):

How can i solve for X(t) in the following equation

diff(mean(X(t)), t) = lambda-delta*mean(X(t))-beta*mean(X(t)*V(t))

and how can i plot X(t) against t.

Here is what I tried.

I need to get the sum of all the real roots of 3*sin(x^2) = x

So I tried this. 

>solve(3*sin(x^2)-x = 0, x, AllSolutions)

Which gaved me this 3*sin(RootOf(_Z-9*sin(_Z)^2))

Then >sum(3*sin(RootOf(_Z-9*sin(_Z)^2)), _Z = -infinity .. infinity)

Which gives me this answer which I cannot make sense of.  signum(sin(RootOf(_Z-9*sin(_Z)^2)))*infinity

I have a question regarding on how to graph and prove that csc ^2(theta)-cot ^2(theta)=1 our professor provided us with a very poor example, i have to graph the left side of the equation first with a 0 to 2*Pi x-axis and then graph the right side on the same plot. 

Could anyone help with solving this integral?

> evalf((∫)[-infinity]^(infinity)((∫)[-infinity]^(v)( uv)/(sqrt(2 Pi))(1/2+1/2 erf((1 )/2 sqrt(2) u))(1/2-1/2 erf((1 )/2 sqrt(2) v))(e)^(-((u^(2+v^(2))))/(2)) ⅆu )dv);

thanks,