Hello everybody,
I have a question please. I would like to resolve an optimization program.
The objective functions are not linear and the constraints are the inequality functions.
Can somebody tell me if we can resolve this type of programs with maple.
Thanks a lot,
Sarah
Do you have some suggestions on how to do this?
> restart;
> K := exp(-b*x+c*exp(x))*(log(x))^a-exp(-x);
> # for this identity, I want to get: {a=0, b=1, c=0} .
> solve(identity(K,x),{a,b,c}) assuming x>10;
I really seem to be having a lot of trouble with Maple that I don't have with Mathcad and I'm hoping you guys can show me what I"m doing wrong. the issue is that I'm trying to graph more than one curve on a grap which is very easy to do in Mathcad but not obviously so with Maple. My teacher told me to setup variables for each function and then use the display property to show the graphs and this works most of the time but then sometimes it does nothing at all but print the functions verbally.
here is an example:
restart;
with(plots):
g1 := implicitplot( 4-(x-2)^2):
g2 := implicitplot( x/2 ):
The CDF function permits you to get a probability, given a random variable x.
But how do I get the random variable x, given the probability?
I have tried using the solver to do this:
df := 29:
chiSquare := RandomVariable(ChiSquare(df)):
solve(CDF(chiSquare, x) = .95, x)
But this code only returns the following warning: "Warning, solutions may have been lost"
But even if the code did work, using a numerical solver for this kind of thing sucks. Surely Maple has a routine that does a revesre lookup?
Can someone please give me some pointers on this?
Jim Cicon
Hello,
Does anybody know if there's any way to display a 2-D plot & a 3-D plot on the same axes? I tried using display and I got the message:
Error, (in plots:-display) cannot display 2-D and 3-D plots together
Regards,
Drew
In the example below, I find integers that solve the equation. That solution happens to be nonnegative. But in general, just plugging in zeros may not give me a nonnegative solution. What is a good way to find a *nonnegative* integer solution, if one exists?
> aa := x^(-4)*exp(-2*x);
> pmu:=product(muval[j]^a[j],j=1..kmax);
Hello,
This is related to my recent posts at
http://www.mapleprimes.com/forum/integral-equations#comment-8339
I have a function of three variables one of which is inside a (difficult) integral. I want to calculate & store the array of points which satisfy the function, then use transform to change those points and then plot them. The function looks like this:
P:=(x,m)->sqrt(x)/(sqrt((x+m)^2-1)*(x+m+sqrt((x+m)^2-1)));
IntegralP:=(m)->Int(P(x,m),x=0..infinity);
f4A:=(beta,Omega,m)->beta^2+(3*(beta*Omega)^(3/2)*evalf(IntegralP(m)))/(4*2^(3/4))-1;
In my earlier case, m=1 (and there was another variable, but it wasn't under the integral) so the integral was not a problem. I tried the same method;
Hello,
I hope someone can help me with this.
1) In Worksheet mode, if I type a:=b+c; and press Enter, Maple echoes
a:=b+c
2) In Document mode, if I type a:=b+c; and press Enter, Maple echoes
b+c
I actually want Document mode to give me the same form of output as Worksheet mode (i.e. I want Document mode to give me a:=b+c as the output). Can this be done?
Thank you for any assistance.
Ty
Still being a newbie to Maple, I am stuck on this one. I am trying to create a general system of changing my equations in an [x,y] coordinate system to a [u,v] system.
As a specific example: I have two substitution equations, u=2x-3y and v=-x+y. I have four equations: x=0, x=-3, y=x, y=x+1. I have tried MapToBasis with both static and procedure statements and I am not having any luck. The equations do not completely change to [u,v]. Here is one example:
> proc (u, v) u = 2*x-3*y end proc;
> proc (u, v) v = -x+y end proc;
> with*VectorCalculus;
> SetCoordinates('cartesian', [u, v]);
I am trying to plot the equations shown in Figures 5 and 6, which are solutions to Equations (14), (15), and (16). These equations are developed in the Lo paper (attached). I am working through the symbolic solution by hand and would greatly appreciate any suggestions on how to set these equations and the graphs of the same in a worksheet. Thank you for the assistance anyone may offer.
Wayne Bell
Download 4865_Lo_p691.pdfView file details
Download 4865_Lo_p692.pdfView file details
I'm having a problem to run this Maple.
Thank you for your help.
> solve({e1⊢e1=e1,e2⊢e2=e2, e1⊣e1=alpha1*e1+alpha2*e2,e1⊣e2=alpha3*e1+alpha4*e2,e2⊣e1=alpha5*e1+alpha6*e2,e2⊣e2=alpha7*e1+alpha8*e2, (e1⊣e1)⊣e1=e1⊣(e1⊣e1),(e1⊣e1)⊣e1=e1⊣(e1⊢e1),(e1⊢e1)⊣e1=e1⊢(e1⊣e1),(e1⊣e1)⊢e1=e1⊢(e1⊢e1),(e1⊢e1)⊢e1=e1⊢(e1⊢e1), (e1⊣e1)⊣e2=e1⊣(e1⊣e2),(e1⊣e1)⊣e2=e1⊣(e1⊢e2),(e1⊢e1)⊣e2=e1⊢(e1⊣e2),(e1⊣e1)⊢e2=e1⊢(e1⊢e2),(e1⊢e1)⊢e2=e1⊢(e1⊢e2), (e1⊣e2)⊣e1=e1⊣(e2⊣e1),(e1⊣e2)⊣e1=e1⊣(e2⊢e1),(e1⊢e2)⊣e1=e1⊢(e2⊣e1),(e1⊣e2)⊢e1=e1⊢(e2⊢e1),(e1⊢e2)⊢e1=e1⊢(e2⊢e1), (e1⊣e2)⊣e2=e1⊣(e2⊣e2),(e1⊣e2)⊣e2=e1⊣(e2⊢e2),(e1⊢e2)⊣e2=e1⊢(e2⊣e2),(e1⊣e2)⊢e2=e1⊢(e2⊢e2),(e1⊢e2)⊢e2=e1⊢(e2⊢e2), (e2⊣e1)⊣e1=e2⊣(e1⊣e1),(e2⊣e1)⊣e1=e2⊣(e1⊢e1),(e2⊢e1)⊣e1=e2⊢(e1⊣e1),(e2⊣e1)⊢e1=e2⊢(e1⊢e1),(e2⊢e1)⊢e1=e2⊢(e1⊢e1), (e2⊣e1)⊣e2=e2⊣(e1⊣e2),(e2⊣e1)⊣e2=e2⊣(e1⊢e2),(e2⊢e1)⊣e2=e2⊢(e1⊣e2),(e2⊣e1)⊢e2=e2⊢(e1⊢e2),(e2⊢e1)⊢e2=e2⊢(e1⊢e2), (e2⊣e2)⊣e1=e2⊣(e2⊣e1),(e2⊣e2)⊣e1=e2⊣(e2⊢e1),(e2⊢e2)⊣e1=e2⊢(e2⊣e1),(e2⊣e2)⊢e1=e2⊢(e2⊢e1),(e2⊢e2)⊢e1=e2⊢(e2⊢e1), (e2⊣e2)⊣e2=e2⊣(e2⊣e2),(e2⊣e2)⊣e2=e2⊣(e2⊢e2),(e2⊢e2)⊣e2=e2⊢(e2⊣e2),(e2⊣e2)⊢e2=e2⊢(e2⊢e2),(e2⊢e2)⊢e2=e2⊢(e2⊢e2)},{alpha1,alpha2,alpha3,alpha4,alpha5,alpha6,alpha7,alpha8});
I wonder if this is the right way to plot following function.
> restart:
> huur:=t->400*(1.01)^ceil(t);
floor(t)
huur := t -> 400 1.01
> plot(huur(t),t=0..50,discont=true);
As you can see I have a the function huur(t) with t a positive integer.
I can't find another way to plot this function except by making use of the floor function.
Isn't there a better way?
Asumming that t is a positive integer didn't worked.
Hello,
in a differential equation (in Maple 10) a variable should be replaced by a constant K so that all differentiations of this variable should evaluate to zero.
This is my code:
constants:= constants, K:
myDG1:= u(t) + diff(u(t),t)^2 + v(t)=0:
subs(u(t)=K,myDG1);
K+(diff(K, t))^2+v(t) = 0
Although I have already assigned K to be constant, Maple doesn't replace d/dt K by 0. Any help would be appreciated.
Consider a sum:
Sum((-1)^(i+1)*Sum(1/j,j=1..i)*x^i,i=1..infinity);
value(%);
The result given by Maple 11 is:
sum((-1)^(i+1)*(Psi(i+1)+gamma)*x^i, i = 1 .. infinity)
However, from an calculus book I know the answer is:
log(1+x)/(1+x).
How can I get this result from Maple? Whether should I give some assumptions?
Best wishes.
Here I difine a subfunction:
mandelbrot := proc(x, y)
local c, z, m;
c := evalf(x+y*I);
z := c;
for m to 30 while abs(z) < 2 do
z := z^2+c
od;
abs(z);
end:
Then I want to plot this function using densityplot command:
plots[densityplot](mandelbrot, -2 .. 0.7, -1.2 .. 1.2, grid=[50, 50], style=patchnogrid, scaling=constrained);#This command is right.
But...
plots[densityplot](mandelbrot(x,y), x=-2 .. 0.7, y=-1.2 .. 1.2, grid=[50, 50], style=patchnogrid, scaling=constrained);#Error, the proc I defined could not determine whether abs(x+1.*I*y) < 2.
