Hi there,

I would like to have an operator (in this case, the natural logarithm) applied to a list/array of points defined as:

ydata := [0.572594976618e-1, 0.327865007249e-1, 0.280821589546e-1, 0.114365745192e-1, 0.578537931608e-2, 0.139154661062e-2, 0.641467839994e-3, 0.18013801847e-3];

How can I apply Maple's ln() operator to the whole array (i.e. avoid to apply it to ydata [1], ydata [2], etc.)?

Thank you,

jon

Hello friends!

I 'm a student and I don't know a lot about Maple, so I would be really grateful if anyone could help me.

I want to solve a system of two equations and I have two unknowns, which are k and εα. However I don't know what I am doing wrong and I can't solve it.

I have attached my file.

Thank you very much in advance!agogos_2.mw

Hi there,
I have the following set of equation which I want to solve using Maple's dsolve command:

d[V(t)*C(t)]/dt = G - K *C(t)
dV(t)/dt = alpha - beta

where V is the volume, C is the concentration, and t is the time variable. G, K, alpha and beta are constant parameters of the problem.

The solution for V(t) is easy to find operating the second equation:
V(t) = V_0 + (alpha - beta)*t

But solving for C(t) is a bit harder.

I would like Maple to solve the system, but the result I get does not really make sense to me. My attempt is this: Maple_artificial_kidney.mw

Any thoughts about how to introduce the equations successfully?

Thank you,
jon

Could anyone assist in rectifying this error ''Error, (in fsolve) {f[1], f[2], f[3], f[4], f[5], f[6], f[7], f[8], f[9], f[10], f[11], theta[11]} are in the equation, and are not solved for''. Here is the worksheet FDM_Revisit_1.mw

I got a problem with solving a second order ODE. 

The ODE is :

-V(xi)+(1/2)*xi*(diff(V(xi), xi))+(1/4)*(diff(V(xi), xi, xi))=-(1/2)*k2*(diff(H(xi), xi))-k1*n*X/E+1+k2

where k1,k2,n,X,E  all are constant.

the condition is :

V(xi) tends to 2*xi^2 as xi tends to infinity.

I used 'dsolve' to solve the equation firstly, and got a solution with two constant C1 and C2, I want to use the condition to elimilate C2, so I used limit(sol,xi=infinity)=2*xi^2. But when I used the command 'limit', I can't get the answer.

Could any one help me? 

Many thanks!!!

Hi there,

I would like to know whether there is some way to print non-integer axis numbers for a given plot.

Let's say I have an array for an independent variable such as [0.5 1 1.5 2]. I would like Maple to use these values in the corresponding axis instead of printing the array indices [1, 2, 3, 4].

I would like to know how this is done for a surface plot.

Attached is the worksheet where "b" is such a parameter/independent variable: MapleSimulation_test11.mw

Thank you,

jon

Hey all, I am new to the Maple software.

I have a question: how to find sum of all elements in a sequence?

For example , A:=[1,2,3,4]  the stupid way to find the sum is add them together by : A[1]+A[2]+A[3]+A[4]=10

Is there a better way to sum them together?

How do I get values out of solved augmented matric in 2D math input? I cannot simply assign the solutions to a value by x:=A[1,5}, because it's an augmented matrix.

Hello,

Since I was working in Matlab with Galerkin method which implies periodic boundary conditions I was wondering how to implement this in maple.

I tried this:

restart;

pde2 := diff(u(x, t), t)+3*(diff(u(x, t)^2, x))+diff(u(x, t),x$3) = 0

IBC := {u(0, t) = u(2, t), u(x, 0) = sech(50*(x-1/2))^2+2*sech(30*(x-1))^2, (D[1](u))(0, t) = (D[1](u))(2, t), (D[2](u))(0, t) = (D[2](u))(2, t)}

pds := pdsolve(pde2, IBC, numeric, time = t, range = 0 .. 2)

But it's telling me: 

Error, (in pdsolve/numeric/process_IBCs) initial/boundary conditions can only contain derivatives which are normal to the boundary, got (D[2](u))(0, t)

So what's wrong?

I don't know what mean is T and V0x?

Please.

I'm trying to write the cumulative distribution function into maple but I seem to get an error. The code is:

>with(Statistics):

>N:=RandomVariable(Normal(0,1)):

>CDF(N,t,inert=true);

I get an error message "Error, (in Statistics:-CDF) unexpected parameters: inert = true".

It used to work, maybe because I was using a newer version of Maple.

Can maple solve maximization problem like

q := proc (a, b, c) options operator, arrow; .2*b+.1*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.2 b + 0.1 c
w := proc (a, b, c) options operator, arrow; .7*a+.1*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.7 ab + 0.1 c
e := proc (a, b, c) options operator, arrow; .7*a+.2*b*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.7 a + 0.2 b c

with(Optimization)

Maxmize(int(min(100+(.7*a+q)*(1/2), a), q)+int(min(100+(.2*b+w)*(1/2), b), w)+int(min(100+(.1*c+e)*(1/2), c), e)-a-b-c-ab-ac-bc)

Error, (in Optimization:-NLPSolve) cannot convert procedures to piecewise

Thanks alot if you can help me.Urgent! Really appreciate.

I am trying to perform the following manipulation (This is a minimum working example).

a < b  < c;
(1)*2;

Error, invalid terms in product: a < b and b < c

 Can anyone tell if it is possible to manipulate inequalities exactly as it is the case with equations?

Dear Maple users,

My problem is as follows:

I have a factor base [2,3,5,7,11,33,34,35,36,37,38,39,40]

The numbers from 2 till 11 are primes, the rest is not. 

Then I have to factor (H+c1)(H+c2) in numbers of the factor base , where c1 and c2 go from 1 to some pre-defined limit. H=32 in my case.
And then I have to put the powers of the numbers of the factor base in a matrix. For example: (H+1)(H+1)=33² but also (H+1)(H+1)=3²*11².

That will become in matrix form [0 , 2, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0 ] but also (!) [0 , 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0 ].

This is not what I want! I want no double representations....

What I want is that (H+c1)(H+c2) should be represented in primes in the matrix if possible and else just represented as the other numbers.

hope you guys can help me!

Hey guys,

I have a problem with this code:

d := Array(1 .. 15);

for i to s1 do

for c2 to 15 do

if `mod`(c2+(J+H)*(`mod`(1/(H+1), A[i])), A[i]) = 0 then d[c2] := d[c2]+l(A[i]) else d[c2] := d[c2]end if

end do

end do;

d;

The problem is that my code stops running when I encounter that `mod`(1/(H+1), A[i]) doesn`t have an modular inverse. What I would like to have is that the loop just continues and does nothing when my array when it encounters an error.

Hope you can help me!