Given the sequence defined by the recursive relation a[n+1] = r*a[n](1-a[n])
You need to use the procedure iterate.
Throughout this problem you should choose initial values in the interval 0<a0<1.
(a) Let r=3/2. Calculate a moderate number of terms in the sequence (between 10 and 20). Does the sequence appear to be converging? If so to what value? Does the limit depend upon your choice of initial value? Plot the terms you have calculated
(b) Let r=2.8. Calculate a moderate number of terms in the sequence (between 10 and 20). Does the sequence appear to be converging? If so to what value? Does the limit
depend upon your choice of initial value? Plot the terms you have calculated How does this sequence differ from that in part (a).
(c) Let r=3.2. Calculate a moderate number of terms in the sequence (between 10 and 20). Show that the sequence does not appear to converging. Plot the terms you have calculated and describe how the sequence behaves in this case.
(d) Consider intermediate values between 2.8 and 3.2 to determine more precisely where the transition in behaviour takes place. Provide a few plots (no more than 4) showing the values you have investigated.
(e) Consider the values of r in the range 3.43<r<3.46. Determine as accurately as you can the value of r for which the period of oscillation doubles.
(f) As r increase further period doubling occurs. Try to find the when the sequence appears to oscillate between 8 values.
(g) Let r =3.65 and calculate a considerable number of terms (at least a few hundred) and plot your values.
(h) For r=3.65 choose a0=0.3 and then a0=0.301. Find and plot some terms in the sequence for each initial value. Determine how long the terms in the two sequences remain close together and when they begin to depart significantly from each other.

I have used data and the curve fitting procedure to find a curve and I now have to put error bars on the graph to complete my task? I know that there is an errorplot function but I cant seem to find a way to do it with the curvefitting thing?

I have attached the whole worksheet, I hope someone can help! Thanks Detectors_Tutorial.mwDetectors_Tutorial.mw

Even if there was some method of exporting the graph to another medium that would also be very helpful

Hi, I'm trying to create a procedure like this;

f:=proc(a);
if a in RealRange(0,1) then a else no;
end if;
end proc;

 But when I try with f(0.5); for instance, I get the error: 

Error, (in f) cannot determine if this expression is true or false: .5 in RealRange(0, 1)

What am I doing wrong here, why can Maple not determine this?

Thank you! 

Hello, 

I am pretty new to Maple and Im trying to find the parameters of an equation using some kind of fit routine but I can only find such a routine to fit an expression to data and not vice versa.

my equation is as follows

epsilon(E):= a0 + a1*ln(E) + a2*ln(E)^2 + a3*ln(E)^3

I do have data for E which I imagine I need. 

E:= {121, 244, 344, 411, 444, 778, 867, 964, 1085, 1112, 1212, 1299, 1408} all in keV :)

Any suggestions/help would be much appreciated, although I am new to the program so go easy on me :) 

Cheers

Hello everyone, I'm a new one to Maple, I've just learnt some basic tools.

I want to creat a command that can animate the graph of line y=ax+b by the parameter a, and b will be subscribe later. For example, I can plot y=x+b by:

F:=b->plot(x+b,x=-1..1):

display(F(1));

It did work.

However, applying this with animation didn't seem to work. 

F:=b->animate(plot,[a*x+b,x=-1..1],a=0..10,frames=5):

display(F(1));

It did not create an animation, instead 5 frames of this graph for a=0, 2.5, 5, 7.5, 10

Please show me a solution for this problem, thank you

Hi

f := x->x^2:

AB(f, 0, 5);

But my proc is not true for below example:

f := x -> x^3-2*x;
AB(f, -1, 4);
I think c has two value and it is not true.

How can I add a condition to my proc that c be between a and b (a<c<b)?

Hi everyone, I'd like to know if it's possible to use the startup code editor as a general editor when writing large programs. It nicely recognizes keywords such as "if" and "proc" and you can indent code using TAB. However I'm not sure whether you can/should use it as a general editor to run code. I've read something about user-profiles related to the startup code but I just don't understand how it works.

Hi all,

I need to solve det[M]=0 for omega.

M is:

M := Matrix(8, 8, {(1, 1) = BesselJ(n, tp*a), (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = -BesselJ(n, tg*a), (1, 6) = -BesselY(n, tg*a), (1, 7) = 0, (1, 8) = 0, (2, 1) = k*n*BesselJ(n, tp*a)/(tp^2*a), (2, 2) = I*k*n*mu0*omega*(diff(BesselJ(n, tp*a), a))/(tp^2*a), (2, 3) = 0, (2, 4) = 0, (2, 5) = -k*n*BesselJ(n, tg*a)/(tg^2*a), (2, 6) = -k*n*BesselY(n, tg*a)/(tg^2*a), (2, 7) = -I*`&mu;g`*omega*(diff(BesselJ(n, tg*a), a))/tg^2, (2, 8) = -I*`&mu;g`*omega*(diff(BesselY(n, tg*a), a))/tg^2, (3, 1) = 0, (3, 2) = BesselJ(n, tp*a), (3, 3) = 0, (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (3, 7) = -BesselJ(n, tg*a), (3, 8) = -BesselY(n, tg*a), (4, 1) = -I*omega*`&varepsilon;p`*(diff(BesselJ(n, tp*a), a))/tp^2, (4, 2) = k*n*BesselJ(n, tp*a)/(tp^2*a), (4, 3) = 0, (4, 4) = 0, (4, 5) = I*`&varepsilon;g`*omega*(diff(BesselJ(n, tg*a), a))/tg^2, (4, 6) = I*`&varepsilon;g`*omega*(diff(BesselY(n, tg*a), a))/tg^2, (4, 7) = -k*n*BesselJ(n, tg*a)/(tg^2*a), (4, 8) = -k*n*BesselY(n, tg*a)/(tg^2*a), (5, 1) = 0, (5, 2) = 0, (5, 3) = k*n*BesselY(n, t0*b)/(t0^2*b), (5, 4) = I*mu0*omega*(diff(BesselY(n, t0*b), b))/t0^2, (5, 5) = -k*n*BesselJ(n, tg*b)/(tg^2*b), (5, 6) = -k*n*BesselY(n, tg*b)/(tg^2*b), (5, 7) = -I*`&mu;g`*omega*(diff(BesselJ(n, tg*b), b))/tg^2, (5, 8) = -I*`&mu;g`*omega*(diff(BesselY(n, tg*b), b))/tg^2, (6, 1) = 0, (6, 2) = 0, (6, 3) = BesselY(n, t0*b), (6, 4) = 0, (6, 5) = -BesselJ(n, tg*b), (6, 6) = -BesselY(n, tg*b), (6, 7) = 0, (6, 8) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = BesselY(n, t0*b), (7, 5) = 0, (7, 6) = 0, (7, 7) = -BesselJ(n, tg*b), (7, 8) = -BesselY(n, tg*b), (8, 1) = 0, (8, 2) = 0, (8, 3) = -I*epsilon0*omega*(diff(BesselY(n, t0*b), b))/t0^2, (8, 4) = -k*n*BesselY(n, t0*b)/(t0^2*b), (8, 5) = I*`&varepsilon;g`*omega*(diff(BesselJ(n, tg*b), b))/tg^2, (8, 6) = I*`&varepsilon;g`*omega*(diff(BesselY(n, tg*b), b))/tg^2, (8, 7) = -k*n*BesselJ(n, tg*b)/(tg^2*b), (8, 8) = -k*n*BesselY(n, tg*b)/(tg^2*b)});

and n=0,1,2.

Except of omega and k ,other parameters is canstant.

After using
with(LinearAlgebra):
detM := Determinant(M):

I used solve(detM=0,omega) and fsolve() but it dosnt work. how can i solve it?

Thanks alot.

Hi Mr Preben Alsholm .

i before asked you about a problem.now i have another question...

for determine unknown parameter oemega(or eigenvalue) according to line ''newsys2 := subs(omega^2 = omega2*10^5,newsys)'' and also line ''approxsoln = [omega2 = 1, f(x) = x^2*(1-x)^2]''

power (5) in term'' 10^5''  in  (omega2*10^5,newsys) and initial guess (1)  (omega2 = 1) are very effective on final result eigenvalues .sometimes gain different eigenvalues ,which it is not impossible recognize that

which one is correct and convege.however i need Positive

 eigenvalue that is minimum between them,If and only if , it converge after
some iteration in two section of maple file which is attached as
Eigenvalue.mw

  According to my code converge occuerd until  eigenvalue gained in first section (before line `11` := ((1+6*alpha2)*(1/12))*(int(fy11^2, x = 0 .. 1))...........) are equal with those obtain in second section .(please see below file for example for_example.mw)

it is necessary mention that between section 1 and 2 is relations.for example amont of first obtained eigenvalue which obtain from first section must be repleaced in line ''approxsoln = [omega2 = 0.661514014001420, h(theta) = theta^2*(1-theta)^2]...........''in second section,

and this procedure should be continued between section 1 and 2 until convergence with desired accuracy occurred.Another relate is that,the first ODE system  can be solved using the first set of boundary data to obtain first estimate for  second section (fy11). At the next step by repeating in the same manner, 

this time by obtained function (g) at the end of maple code (Solution of the second ODE system  with second set of boundary conditions leads to the first estimate of function g),that at this stage, the first iteration is completed.

Now by replacing omega2 which is determined in this section in to the first section ,fy11 is updated and gained Further.

Next, the updated function , by continuing the iterative procedure.

in matlab bvp4c rule is used for this purpose.is this impossible in maple software ?if not please help me for solve and gain correct omega..

thanks alot

Hello,

I use Maple 18. How can I produce a table of 50 values of the fanction f(x) in the interval [a,b]? size step is (b-a)/n

hi, is there any maple command that will allow me to test the speed of a procedure? something like the ctime in c++ or time in java.thanks

and why are there decimal points after certain numbers?

I'm taking 2nd semester calculus and in this maple lab, I can't seem to get Maple to evaluate this correctly. Any help is appreciated.

My friend did this and got the first constant, 392..., and just a², so it was 392.07a²

I using Maple 18.

I write below procedure,

Eqplot:=proc(a,b)
local y,x;

y:=a*x+b;

y,plot(y,x);

end proc:

when I use Eqplot(2,3) , then I will expect the result for both of equation 2x+3 and sketch y=2x+3, but unfortunately the result is only sketch y=2x+3.

please hint me about that how I can write a proc that result be both of 2x+3 and sketch 2x+3 in maple 18?

Perhaps the question is trivial, but I could not find the solution.

I am solving numerically an ODE (e.g., a simple harmonic oscillator) with the righthand side that contains a random part. For example it is

eq:=diff(a(t),t,t) + a(t) = (1+0.01*R(t))*cos(5*t);

where R(t) is a random function of t. How can I make such a function?

The naive attempt: 

eq:=diff(a(t),t,t) + a(t) = 3.*(1+0.1*rand()/1000000000000.)*cos(5.*t);

gave me a fixed (while random) value, e.g.

...

But, I need this coefficeint to be random each time step.

Any suggestions are very welcome!