I tried to get the maximum and minimum values of the following function. From the plot I get them but its not accurate. Please advise me to get them accurate.

F:=0.85:B:=0.5:

K:=N->(N*(1+F*N/(N^2+B^2-F*N)));

implicitplot(((N^2+B^2-F*N)*K=N*(N^2+B^2-F*N+F*N),K=0..10,N=0..10,view=[0..5,0..4],numpoints=90000,axes=boxed,thickness=2,color=black,font=[1,1,20],tickmarks=[3, 3],linestyle=1));

I gave this function:

g := -2-k[1](lambda*alpha[2]*k[2]+alpha[1](-2-2*k[2]+k[2]*lambda^2))+k[1](lambda*alpha[2]*k[2]+alpha[1](-2-2*k[2]+k[2]*lambda^2))*lambda(lambda*alpha[2]*k[2]+alpha[1](-2-2*k[2]+k[2]*lambda^2))^2

I would like to factor out or to collect:   (lambda*alpha[2]*k[2]+alpha[1](-2-2*k[2]+k[2]*lambda^2))

I use following command: collect(g, lambda*alpha[2]*k[2]+alpha[1](-2-2*k[2]+k[2]*lambda^2))

Nevertheless, I receive an error ...

Error, (in collect) cannot collect lambda*alpha[2]*k[2]+alpha[1](-2-2*k[2]+k[2]*lambda^2)

Could you help me to solve my issue ? 

What am I doing wrong ? Do I have to use another command for this ?

I would like to thank you in advance.

Best regards,

D.L.

q[1] = sqrt(x)*alpha-lambda(sqrt(x)*alpha-lambda*q[1]*q[3]-p[2])*(sqrt(x)*alpha-lambda*q[1]*q[2]-p[3])-p[1]

I am looking for q[1] solution. When I solve for q[1], maple gives me following answer:

q[1] = RootOf(-_Z+sqrt(x)*alpha-lambda(sqrt(x)*alpha-lambda*_Z*q[3]-p[2])*sqrt(x)*alpha+lambda(sqrt(x)*alpha-lambda*_Z*q[3]-p[2])*lambda*_Z*q[2]+lambda(sqrt(x)*alpha-lambda*_Z*q[3]-p[2])*p[3]-p[1])

Is it possible to obtain a classical solution for the calculations above. (can not understand the meaning of: RootOf and _Z. I need q[1] in order to solve further in my system of eqautions for  q[2],  q[3]

could you help me please to find a solution for this issue...
I would like to thank you in advance 
Best regards,

D.L.

((-1+lambda(lambda*k[2]*alpha[2]+alpha[1](k[2]*lambda^2-2*k[2]-2))^2)*k[1](lambda*k[2]*alpha[2]+alpha[1](k[2]*lambda^2-2*k[2]-2))+(-lambda*lambda(lambda*alpha[1]*k[1]+alpha[2](k[1]*lambda^2-2*k[1]-2))^2+lambda)*k[2](lambda*alpha[1]*k[1]+alpha[2](k[1]*lambda^2-2*k[1]-2))-k[1]*(-alpha[1]+lambda*alpha[2])*lambda(2+5*k[2])^2+lambda^5*alpha[2]*k[1]*k[2]-alpha[1]*k[1]*lambda^4*k[2]-2*lambda^3*alpha[2]*k[2]+2*alpha[1]*k[2]*lambda^2+(4*alpha[2]*(1+k[2])*k[1]+4*alpha[2]+4*k[2]*alpha[2]+2)*lambda-4*alpha[1]*(1+k[2])*k[1]-4*alpha[1]-4*alpha[1]*k[2]-2)*(-(lambda(lambda*k[2]*alpha[2]+alpha[1](k[2]*lambda^2-2*k[2]-2))-1)*(lambda(lambda*k[2]*alpha[2]+alpha[1](k[2]*lambda^2-2*k[2]-2))+1)*(-1+k[1]*lambda^2-k[1])*k[1](lambda*k[2]*alpha[2]+alpha[1](k[2]*lambda^2-2*k[2]-2))+(-lambda*lambda(lambda*alpha[1]*k[1]+alpha[2](k[1]*lambda^2-2*k[1]-2))^2+lambda)*k[2](lambda*alpha[1]*k[1]+alpha[2](k[1]*lambda^2-2*k[1]-2))-k[1]*(-alpha[1]+lambda*alpha[2])*lambda(2+5*k[2])^2+lambda^5*alpha[2]*k[1]*k[2]-alpha[1]*k[1]*lambda^4*k[2]-2*lambda^3*alpha[2]*k[2]+(2*k[1]+2*alpha[1]*k[2])*lambda^2+(4*alpha[2]*(1+k[2])*k[1]+4*alpha[2]+4*k[2]*alpha[2]+2)*lambda-(4*(1+k[1]))*(alpha[1]+alpha[1]*k[2]+1/2))*x

I am looking to simplify this long term...

could you help me please to find a solution for this issue...

I would like to thank you in advance

Best regards,

D.L.

Hi,

1-Triying to plot a function divided by its maximum value,sometimes it works with some parameters that means, the max.value of the plot is 1.

But when i change the data the max. value in the plot in graeter than 1 which is wrong!! should be 1.

dont know why??

2- Changing different data in the parameters, the programme takes long long time then i stop it?

please help me with these two problems.

restart:
>
------------------------- Defining the nature of the variables used ----------------------
assume(T,real):Digits:=25:n:=1:tau:=Pi:
theta:=0:phi:=0:
lambda:=n;Omega:=1:Gamma:=0.01:
--------------------- Input---------------------------------
1

J1

term1:=(exp((Gamma+I*d)*tau)-1)/(2*(Gamma+I*d)):
Ak1:=d->(exp((Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))+(exp((Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma)):
Ak2:=d->(exp((Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))+(exp((Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma)):
term2:=(evalf(-0.25*sum(BesselJ(k,Omega*Gamma/(4*n))*Ak1(d)+BesselJ(k,-(Omega*Gamma)/(4*n))*Ak2(d),k=0..50))):
J1:=(term1+term2):
J1mod:=(Re(J1))^2+(Im(J1))^2:
###### J2#########################
Ak1:=d->(exp((Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))-(exp((Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma)):
Ak2:=d->(exp((Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))-(exp((Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma)):

J2:=(evalf(-0.25*sum(BesselJ(k,Omega*Gamma/(4*n))*Ak1(d)+BesselJ(k,-Omega*Gamma/(4*n))*Ak2(d),k=0..100))):
######################

J2mod:=(Re(J2))^2+(Im(J2))^2:
J3 same as J1differ in sign
term1:=(exp((Gamma+I*d)*tau)-1)/(2*(Gamma+I*d)):
Ak1:=d->(exp((Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega-2*lambda*k/Gamma))+(exp((Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega+2*lambda*k/Gamma)):
Ak2:=d->(exp((Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d-0.5*Omega-2*lambda*k/Gamma))+(exp((Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma))*tau)-1)/(Gamma+I*(d+0.5*Omega+2*lambda*k/Gamma)):
term2:=(evalf(0.25*sum(BesselJ(k,Omega*Gamma/(4*n))*Ak1(d)+BesselJ(k,-Omega*Gamma/(4*n))*Ak2(d),k=0..100))):
J3:=term1+term2:
J3mod:=(Re(J3))^2+(Im(J3))^2:
J4 same as J2 but -0.25-->2

J4:=-2*J2:
######################

J4mod:=(Re(J4))^2+(Im(J4))^2:

calculate the spectrum

Spec:=d->(exp(-2*Gamma*tau)*(J1mod*cos(theta/2)^2+J2mod+J3mod*sin(theta/2)^2-0.5*Re(J3*J4*sin(theta)*exp(I*phi))+0.5*Re(J1*J4*sin(theta)*exp(-I*phi)))):

with(plots):

tit:=sprintf("l=%g,W=%g,G=%g",lambda,Omega,Gamma):
Smax1:=max(seq(evalf(Spec(d)),d=-100..100)):
plot(evalf(Spec(d)/Smax1),d=-15..15,axes=boxed,title=tit,color=black,font=[2,3,18],thickness=2,tickmarks=[3,3],titlefont=[SYMBOL,14],font=[1,1,18],linestyle=1);

I have two equations:

and I need to extract the coefficients of x0 and xp0 (to build a matrix; but that is not the issue). So I use coeffs:

coeffs(op(2,X),x0);

coeffs(op(2,XP),x0);

xp0

While I can live with the first result (for X) (not that I like it), the fact that the 2nd one (for XP) has a completely different structure prevents any kind of algorithmic extraction of the coefficients (and the second one of course is 0) for further use. I tried the form of coeffs with a third argument (a name which gets the result assigned) but the same result. This example is a real case, obviously a trivial one and others will be much more involved so I really would like this to work. And yes, I did "collect" before using coeffs (not that it was needed here).

Any ideas out there?

Mac Dude.

solve({[(alpha[1]-alpha[2]*lambda)*sqrt(x)+p[2]*lambda]*[k[1]*(1-lambda^2)+2] = p[1]*(2*k[1]*(1-lambda^2)+2), [(alpha[2]-alpha[1]*lambda)*sqrt(x)+p[1]*lambda]*[k[2]*(1-lambda^2)+2] = p[2]*(2*k[2]*(1- lambda^2)+2)}, [p[1], p[2]])

Warning, solutions may have been lost

could you help me please to find a solution for this issue...

I would like to thank you in advance 

Best regards,

D.L.

I have Maple 15 and I am transferring it to a different computer but I cannot find the download link anywhere. Where would I find it? Thank you!

Hello,

I just installed maple 15 as well as the latest update. However, when I was trying to plot, it gives me this:

Which had never happened before I reinstall my operation system.

Could anyone help me please? Thanks a lot!

Hi. I have a question it goes.

Find the equation f(x) which satisfies the following conditions. f'(x)= 1/(1+cot(x)) and f(pi/4)=1.

So I enter the following in Maple and this is what I end up with.

http://s5.postimg.org/eriq99rs7/Maple.jpg

This should work as far as I can tell. The final diff command should return the original f'(x), but it doesn't. Can anyone tell me why?

I have used this exact method before on similair problems and it...

Himmelblau.mw

     On the basis of Dragнilev method…

     Is there anyone interested in the algorithm to reduce the distance between the points of the given constraints? The algorithm is adapted for use in R ^ n. This is an example of its work on the surface:  
f = - (x1 ^ 2 x2-.3) ^ 2 - (x1 x2 ^ 2-.7) ^ 2 - 5;  

     Approximate description of the algorithm in pictures.

Greetings to all.

I have been using the numtheory package for quite some time now and it has helped me advance on a number of problems. Recently an issue came to my attention that I have known about for a long time but somehow never realized that it can be fixed. This is the fact that the numtheory package does not know about Dirichlet series, finite and infinite. Here are two links:

Hello,

it is more or less my first time using maple for my research. During my mesurements I collect data and store those in an excel sheet. The I made calculations row by row.

I figured out how to import the excel sheet and calculate row by row. Maple calculates the correct results, but how can I easily store the results again, for example row by row in an excel sheet or csv?

This is my code:

> restart;#Range of the Matrix A1:KXXXXprint(??);

Hi, I want to know how to change the rectangular form to polar form, and how to manipulate the operation of them, for example the picture below. The angle of polar form contain radian and degree, how to represent them in Maple?Thank you.

1.

   Suppose thT the population y(t) of a certain kind of fish is given by (D(y))(t) = Ay(t)-By(t)^2, and fish are caught at a rate Hy proportional to y.Solve this so-called schaefer model. Find the euilibrium solutions y1,y2(>0) when H

2.< Intermittent harvesting>

 In 1. assume that you fish for 3 years, then fishing is banned for the next 3 years. Thereafter you start again. And so on. This is called intermittent harvesting....