Hi everyone

I have solved an ODE in Maple with the following command:

soln := dsolve({ode1, f(0) = 0, f(1) = 1, (D(f))(0) = 0, (D(f))(1) = 0}, f(eta), numeric)

where ode1 is:

ode1 := diff(f(eta), [`$`(eta, 4)])+R*(eta*(diff(f(eta), [`$`(eta, 3)]))+3*(diff(f(eta), [`$`(eta, 2)]))-f(eta)*(diff(f(eta), [`$`(eta, 3)]))+(diff(f(eta), eta))*(diff(f(eta), [`$`(eta, 2)]))) = 0

then I saved the soln plot in p1 as:

p1 := odeplot(soln, [eta, f(eta)...

Hi everyone,

Can somebody help me to make this code faster? Or if somebody have fast computer and can execute it and give me the result.

When I execute it, I don't have the result of the last command and get the error of not enough memory allocation or something.

Here is the code:

>restart;with(LinearAlgebra):with(VectorCalculus):with(linalg):
> a1:=1:
> a2:=1:
> T:=<x1,x2,x3,x4,y1,y2,y3,y4,u1,u2,u3,u4>:
> 

I have created the following:

a := proc (b, r)
local x, n;
x := [seq(b, x = 0 .. 99)];
for n from 1 to 99 do x[n+1] := r*x[n]*(1-x[n])
end do;
return x
end proc;

> a(.5, 2.5);

This returns a list of points in a list representing a population at each generation from x_0 to x_99.  How can I get Maple to plot these points as a time series?

I have the following procedure that I am working with. I would like to store the result of the procedure in a variable called "result".

[code]
f := proc(result, x)
     if x <= -7 then
         result :=  1:
     elif x > -7 and x <= 7 then
         result := 2:
     elif x > 7 then
 ...

The new VolumeOfRevolution command in the Student[Calculus1] package has lots of nice features for controlling the output (Riemann sum, 3D plot, animation, inert integral, numerical approximation, ...).

The default colors for the plot are not very appealing. Not only is the pink easily washed out by a projector in a moderately light room, but the detail is missing. To illustrate, consider

with(Student[Calculus1]):
VolumeOfRevolution(2-x^2, x, x=0..1, output=plot );

Hi,

I want to put some code into a Maple library. I'm new to producing packages. The following code runs OK, but when I try to use one of the procedures added by the LibraryTools[Save] command, Maple returns

input: MySum(a, b);
output: MyMapleLibrary:-MySum(a, b)

Why does Maple return an unevaluated call to MySum?

The code is below:

restart:
with(LibraryTools):

MyMapleLibrary:=module()
  description "A sample Maple library.":

Hi,

In chapter 2 of the Advanced Programming Guide, it says:

"Important Always ensure that the standard Maple library is write-protected to avoid saving expressions in it. If you accidentally save something to the standard Maple library, you may need to restore the original from media on which you obtained the Maple software".

How do I make the standard Maple library write-protected? Is it by marking as read-only the following folder?

How can I see the steps of solving a PDE ( partial differential equation ) in maple .
for example see the sepration variables and other equations that lead to final Eq .

a=5/6

then how to plot the function in Maple 12 

restart:with(plots):

eq := diff(f(eta), eta, eta, eta)+.11-(diff(f(eta), eta));

bc2 := f(0) = 0, (D(f))(0) = 1, (D(f))(10) = .1;

dsolve({eq,bc2}):

F1:=evalf(rhs(%));

contourplot(x*F1,x=-3..3,eta=-0.2..0.2);

I need to get the data from the above plot in a .txt file.

Thanks  

I have been given the logistics map:   x_n+1 = rx_n(1-x_n) where x_n is a number between zero and one that represents the size of the population in the nth generation as a ratio of the environment’s maximum capacity.
The parameter r is a positive number related to both the reproduction rate and the rate at which the
population dies due to limited resources.  Create a procedure (taking arguments r and initial population size x

Is this really the answer we want?  

> sum(log(n)/n,n=1..infinity);  

Error, (in Zeta) numeric exception: division by zero  

rsolve(a(n)=2/a(n-1)+1,a(1)=1),

thx in advance.

Sorry, this seems like a silly question. But is there an easy way to convert trig funcitons, or even non trig functions to orthogonal (in this case Legendre) polynomials? 

Hello,

I am trying to solve 9 non linear equations in 3 do loops(program attached without loop) and I am using digits:=20; for computing 10 values its takes around 5 minutes. I need to compute around 18000 such values for a plot. It might take several days for this. Is there any mistake in my code which is taking it so long to run or is it just the complexity of the equations.

Secondly, I am trying to compute some calculations to...