MaplePrimes Questions

Nt+1 = RNt/ [1+ (R-1/K)Nt] where R> 1, K>0. When No >0 , lim t approaching infinityNt = K for all values of R>0.

Find Nt for t=1,2,3...10 for K= 100 and No= 20 when a) R=2, b)R=5 and c)R= 10.

I have managed to get the answers for a), b) and c) manually.

But how do I plot Nt as a function of t for the three choices of R in one coordinate system?

The pH levels of a lake controls the conc. of harmless ammonium ions(NH4+) and toxic ammonia (NH3) in the lake.

For pH levels <8, conc. of ammonium ions are little affected by pH changes, but decline over many orders of magnitude as pH levels increase beyond pH 8.

Toxic ammonia are negligible at low pH , increase over many orders of magnitude as the pH level increases and reach a high plateau at about pH =10 (after which, NH3 are little affected by pH changes).

 

How do I illustrate this graphically?

 

Percentage survivorship is a function of time; there ia negligible decline in survivorship at pH=6, marked decline in survivorship at pH=3.5, resulting in no survivors after just 8 hours.

What is the function for this problem and how do I plot it in Maple?

 

Hi all,

I have the following expression,

Typesetting[delayDotProduct](b, K[P], true)+k*K[D]-(Typesetting[delayDotProduct](b, m[1], true)+Typesetting[delayDotProduct](b, m[2], true)+K[D]*m[2])*k*K[P]/(Typesetting[delayDotProduct](b, K[D], true)+k*m[1]+k*m[2]+K[P]*m[2]-m[1]*m[2]*(Typesetting[delayDotProduct](b, K[P], true)+k*K[D])/(Typesetting[delayDotProduct](b, m[1], true)+Typesetting[delayDotProduct](b, m[2], true)+K[D]*m[2]))

where K and K_P are the controller gains, k is the stiffness, b is the damper, m1 and m2 are masses.

How can I find the condition on variables such that the numerator of this expression is greater that zero?

The conditions should appear as inequalities.


 

b*K[P]+k*K[D]-(b*m[1]+b*m[2]+K[D]*m[2])*k*K[P]/(b*K[D]+k*m[1]+k*m[2]+K[P]*m[2]-m[1]*m[2]*(b*K[P]+k*K[D])/(b*m[1]+b*m[2]+K[D]*m[2]))

b*K[P]+k*K[D]-(b*m[1]+b*m[2]+K[D]*m[2])*k*K[P]/(b*K[D]+k*m[1]+k*m[2]+K[P]*m[2]-m[1]*m[2]*(b*K[P]+k*K[D])/(b*m[1]+b*m[2]+K[D]*m[2]))

(1)

``


 

Download Routh-Hurwitz_stability.mw

 

The rate of resource consumption by organisms,v depends on resource concentration, S. The Blackman model of resource consumption assumes a linear relationship between v and S. Below a threshold concentration, Sk , consumption rate increases linearly with S=0 when v=0; when S=Sk, v reaches its max value, vmax. For S> Sk, v stays at vmax.

The function like this is defined piecewise,

  v= { g(S)   for 0<=S< Sk

         vmax  for S>=Sk

g(S)=( vmax/Sk).S

How do I graph v as a function of S with the graph showing the point where Sk=vmax?

Egg development

It takes 3 days at 20oC, but almost 20 days at 5oC for an egg to develop and hatch. When graphed on log -log plot, it says that egg development time (in days) as a function of temperature (in oC) is a straight line.

 

How do I sketch the log-log graph in Maple?

The chemical reaction is given as

A+2B ---> AB2

The reaction occurs ina closed vessel and the initial concentration of A=[a] and B=[b]

The reaction rate is given by

R(x)= k(a-x)(b-2x)2

where x=[AB2]

Graph R(x) for the relevant values of x when a=5,b=6 and k=0.3.

 

Need help creating a proc() for Kaprekar's routine. Ive been able to convert a number (x,n) to ascending/descending order and then defining f to be the sum of them subtracted, but i cant figure out how to put it into a proc so that it continues to loop until kaprekar's number 6174 is reached and stop.

If you have any idea please advise! thanks

Im trying to create a groebner basis to solve a series of equations to make a basis. The values can only be 1 to 4 for a..p which im also unsure if/how to define which may be causing the problem.

 

 My equations are as follows:

 

e1 := (((a-1)(a-2))(a-3))(a-4) = 0;
e2 := (((b-1)(b-2))(b-3))(b-4) = 0;
e3 := (((c-1)(c-2))(c-3))(c-4) = 0;
e4 := (((d-1)(d-2))(d-3))(d-4) = 0;
e5 := (((e-1)(e-2))(e-3))(e-4) = 0;
e6 := (((f-1)(f-2))(f-3))(f-4) = 0;
e7 := (((g-1)(g-2))(g-3))(g-4) = 0;
e8 := (((h-1)(h-2))(h-3))(h-4) = 0;
e9 := (((i-1)(i-2))(i-3))(i-4) = 0;
e10 := (((j-1)(j-2))(j-3))(j-4) = 0;
e11 := (((k-1)(k-2))(k-3))(k-4) = 0;
e12 := (((l-1)(l-2))(l-3))(l-4) = 0;
e13 := (((m-1)(m-2))(m-3))(m-4) = 0;
e14 := (((-1+n)(n-2))(n-3))(n-4) = 0;
e15 := (((-1+o)(o-2))(o-3))(o-4) = 0;
e16 := (((p-1)(p-2))(p-3))(p-4) = 0;
e17 := a+b+c+d-10 = 0;
e18 := a*b*c*d-24 = 0;
e19 := e+f+g+h-10 = 0;
e20 := e*f*g*h-24 = 0;
e21 := i+j+k+l-10 = 0;
e22 := i*j*k*l = 24;
e23 := m+n+o+p-10 = 0;
e24 := m*n*o*p-24 = 0;
e25 := a+b+e+f-10 = 0;
e26 := a*b*e*f-24 = 0;
e27 := c+d+g+h-10 = 0;
e28 := c*d*g*h-24 = 0;
e29 := i+j+m+n-10 = 0;
e30 := i*j*m*n-24 = 0;
e31 := k+l+o+p-10 = 0;
e32 := k*l*o*p-24 = 0;
e33 := a+e+i+m-10 = 0;
e34 := a*e*i*m-24 = 0;
e35 := b+f+j+n-10 = 0;
e36 := b*f*j*n-24 = 0;
e37 := c+g+k+o-10 = 0;
e38 := c*g*k*o-24 = 0;
e39 := d+h+l+p-10 = 0;
e40 := d*h*l*p-24 = 0;
e41 := c = 4;
e42 := e = 1;
e43 := j = 2;                  
e44 := p = 3;
 

Hi

This question may have been asked before but I cannot find it!

I am currently have 3 lists:

L:= [1,3,5,2,3]

M:=[4,2,4,6,2]

R:= [5,3,2,1,4]

I would like to be able to come up with a way (I assume a procedure) where I can order L and M based on the numerical order of R. So in this case that would be:

R:=[1,2,3,4,5]

L:=[2,5,3,3,1]

M:=[6,4,2,2,4]

I hope this makes sense? Any help would be greatly appreciated!

 

Hi All,

I am ploting with the command of:

logplot([1-R(Pp, tp, tf, 'qN(Pp, tp, tf)'), 1/(1+1+tp), .5/(.5+.5+tp), (1-Pp)*Beta(2, 1+tp)/(Pp*Beta(1, 1)+(1-Pp)*Beta(1, 1+tp))], Pp = 0 .. 1, style = pointline, numpoints = 1, gridlines = true, symbolsize = 10, legend = ["Our worst-case prior with Pp", "Flat prior Beta (1,1)", "Jeffrey's prior Beta (0.5,0.5)", "Normalised flat prior Beta (1,1) with Pp"], labels = ["Pp", "Posterior expected pfd"], symbol = [asterisk, diagonalcross, box, circle], title = "Given fixed t_past=100")

and get:

 

which is ok but the legends are wrong. where are the symbols in the legends? Am I missing something....

 

how to solve differential equations with dependent functions in maple? 

with(PDEtools):
V:=-exp(I*(k*x+omega*t+theta))*sqrt((k^2+omega)/(k*sigma))*sech(sqrt(k^2+omega)*(-2*k*t+x));
pde[1] := simplify(I*(diff(V, t, t))+diff(V, x, x)-I*sigma*V*(conjugate(V)*(diff(V, x))-V*(diff(conjugate(V), x)))) = 0;

 

Dear sir, can we get a plot of graph with line style like  line+symbol so on.. 

I am attaching codes and sample graph .

 

h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=(4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4:
lambda:=(F,Nb,delta2)->Int(K1,z=0..1):

L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
A:=plot( [seq(seq(lambda(F,Nb,delta2), Nb=[0.1,0.2,0.3]), F=[0.1,0.2,0.5])],
          delta2=0.02..0.1, linestyle = [solid,longdash,dashdot],'
          thickness = 2',color=[red$3,blue$3,black$3],
          legend=["Nb=0.1","Nb=0.2","Nb=0.3",""$6]):
B:=plots:-textplot([[0.05,-1,F=0.1],[0.05,1.5,F=0.2], [0.05,3,F=0.5]], font=[times, 14]):

plots:-display(A, B);

Hi, I want to convert differential operator (D) to diff in below file, but "convert" doesn't work properly. I think beacuse both "D" and "diff" exist in the answer or mabey there is an other reason. what should I do?

I attached file here.1.mw

thanks in advance.

First 844 845 846 847 848 849 850 Last Page 846 of 2434