> restart;
> with(LinearAlgebra);
> tolerance = 10^(-3); x1km1 := 0; x2km1 := 0; x3km1 := 0; for k to 2 do k := k; x1k := 1/4*(5-x2km1+x3km1); x2k := 1/3*(-4+x1km1-x3km1); x3k := 1/5*(1-2*x1km1-2*x2km1); evalf(x1k); evalf(x2k); evalf(x3k) end do;
                               1
                               5
                               -
                               4
                               -4
                               --

q1.mw

This equations should have numerical solution, but they can’t be calculated. The result is too long, and include ‘RootOf’

It just like this:

> t := 15000; S := solve({n = b+2*c+f+2*g, (a+b+c)/(a+b+c+p+q+f+g) = 5*(1/10), q^2-3.47*p*10^27*t^(3/2)*e^(-10.97*10^4/t) = 0, n*b-(2.407*a*10^21*8)*t^(3/2)*exp(-18.27*10^4/t) = 0, n*c-(1/2)*(2.407*b*10^21*5)*t^(3/2)*exp(-32.035*10^4/t...

Regards,

I need help with solving an equation. I have the following equation:

24*b[0]+24*b[1]*xi+24*b[2]*xi^2+(2*(-12+24*xi))*(b[1]+2*b[2]*xi)+(2*(-12*xi+12*xi^2))*b[2]+P[cl]*L^2*(-12*xi+12*xi^2)

to this equation I applied the comand: series(%,xi) and as a result obtained:

24*b[0]-24*b[1]+(72*b[1]-72*b[2]-12*P[cl]*L^2)*xi+(144*b[2]+12*P[cl]*L^2)*xi^2

My question is: the way to solve this equation in the manner I need is setting the coefficients...

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.