Download compre1try3.mw

I am facing the same problem with Maple 17. Which is also discussed here

https://www.mapleprimes.com/questions/222168-Maple-Not-Evaluating-Symbolic-Variables?reply=reply

But the link provided below is dead
https://www.maplesoft.com/support/faqs/detail.aspx?sid=34802

Can you recall how to solve it?  I have tried your solution below, but it did not solve it.  

"@peacess If the file isn't there you can just create a new file with a text editor (e.g. Notepad) and add the required option.
Or you can use the maplejava.l4j.ini file inside the .zip file I've attached:  maplejava.l4j.zip"

Thank you.  

This is a very general, basic question.  

When I supply an input is the calculuation/computation performed on my computer?  Is it uploaded to a Maple server for calculating? 

I am assuming it takes place on my computer but just wanted to double check.  

Also, is there a way to see the progress of a calculation?  If it is something that might take days or even weeks..

Hi there

I am trying to solve the ODE below

schechter_guo_v2.mw
 

Download schechter_guo_v2.mw

My questions are:

1. Why does the solution for the longer time span (t<10) looks different from the shorter time span (t<0.1)? I have read about stiff ODEs and probably this has something to do with it. I am trying to figure out what is going on with the solutions.

2. I tried dsolve with stiff methods (lsode and rosenbrock) and both gave me the same solutions as above. I have not tried the advance options yet. How do I set dsolve so that the solution for the longer span (t>10) looks similar to shorter span?

Many thanks for your answers/suggestions.

I am trying to use standard notation for machine constants such as Xd' and Xd'', the d being a subscript. The prims are interperted as differential operators. How can I disable the d/dx.   

I've followed all the recommendations, but I can't see the mistake yet...:

z := proc(n::integer)

local A,E,S,k,s1,s2;

if n = 0 then

matrix(2,1,[0,0]);

else if n = 1 then

matrix(2,1,[-2,epsilon^()+sqrt(epsilon^(2)+4)]);

else

A(n) := matrix([[1+epsilon^(2)-lambda^(n),epsilon^()],[epsilon^(), 1-lambda^(n)]]);

E := matrix(2,1,[epsilon^(2),epsilon^()]);

S := 0;

for k  from 1 by 1 to n-1 do

s1 := thisproc(n-k);

s1 := s1[1];

s2 := thisproc(k);

s2 := s2[1];

S := S+s1*s2;

end do;

A(n)^(-1)*E*S ;

end if;

end proc

I need this to do a Dynamical systems problem called Hénon map, which has by function formula:

f(x,y) = (x+epsilon*(y+epsilon*x*(1-x)), y+epsilon*x*(1-x))

so the problem consists in finding the expression of function P such that f(P(t)) = P(lambda*t), where lambda = 1+(epsilon/2)*(epsilon-sqrt(epsilon^2+4)) is one of the eigenvalues of the differential matrix of f valued on (0,0). In Maple, I've put z = (x,y) to let the problem be a vector problem. And as I made the calculations with the help of a LaTex-presentation PDF, I took the formula A(n)*z(n) = (epsilon^2;epsilon)*S, where s is the following recursive sum in the code. If I don't get the result, I won't be able to explain why Matlab plotting goes too far, even when n = 2!!! Bref, trying to put things much easier in a difficult problem, but Maple returns me the error like I did some mistake, and I tried to debug these code lines!!! Along the code, please note that lambda and epsilon are mere parametres, not local variables. If the error comes from there, please put epsilon = 0.2 and lambda as the formula just below the code.

Dear!

Hope everyone is fine. I want to write the values of "phi" with their line colors in legend or in title. I know how to make legends but the problem is here there are nine graphs and we need only three legends. Please fix my problem. I am waiting the positive response. Thanks in advance

Legend.mw

I want to solve this integer-linear programming problem. There are four variables (xp,xq,xifp,xifq)
Two variables are nonnegative and integer (xp,xq)

the other two variables are binary(xifp,xifq)
 

When executing this command, maple says that the problem appears to be unbounded

LPSolve(750*xp+250*xq-50000*xifq,maximize=true,
{20*xp+30*xq<=2600,
30*xp+20*xq<=1000,
20*xp+30*xq<=2000,
xq<=99999*xifq,
xq>=20*xifq,
xp<=99999*xifp,
xp>=30*xifp,
xp>=0,
xq>=0},
integervariables={xp,xq},
binaryvariables={xifp,xifq}
);

and gives this solution[0, [xifp = 0, xifq = 0, xp = 0, xq = 0]]

however if i set a range for the variables xp,xq from 0 to 99999 for example

LPSolve(750*xp+250*xq-50000*xifq,maximize=true,
{20*xp+30*xq<=2600,
30*xp+20*xq<=1000,
20*xp+30*xq<=2000,
xq<=99999*xifq,
xq>=20*xifq,
xp<=99999*xifp,
xp>=30*xifp,
xp>=0,
xq>=0},
integervariables={xp,xq},
binaryvariables={xifp,xifq},
xp=0..99999,
xq=0..99999
);
maple returns this solution 

[49950000, [xifp = 1, xifq = 1, xp = 50000, xq = 50000]]

but if xp=50000 and xq=50000 this doesnt satisfy the first constrain
because 20*xp+30*xq</=2600
why does maple gives this mistaken solution, how can i solve it ?
please help,thanks for reading

hi every body: 

i want to solve the below ODE with boundry conditions ics, can anybody help me?

The ODE:

ode := x^3*(diff(f(x), x, x, x, x))+alpha*(x^4*(diff(f(x), x, x))+x^3*(diff(f(x), x, x))-x^2*(diff(f(x), x)))-2*x^2*(diff(f(x), x, x, x))+3*x*(diff(f(x), x, x))-3*(diff(f(x), x))+R*x*f(x)^2-R*x^2*(diff(f(x), x, x))*(diff(f(x), x))-3*R*x*f(x)*(diff(f(x), x, x))+3*R*f(x)*(diff(f(x), x))+x^2*R*f(x)*(diff(f(x), x, x, x))-M^2*(x^3*(diff(f(x), x, x, x))-x^2*(diff(f(x), x))) = 0

The boundry conditions:
ics := f(0) = 0, f(1) = 1, (D(f))(1) = 0, limit(diff((diff(f(x), x))/x, x), x = 0) = 0

M,R,alpha are constants.

with regards... 

I would like procedure P to perform inside the map statement exactly as it does outside the map statement. Can this be done?

Map_Procedure.mw

Hi,

I have problem with my first DataTable component ( no interaction with my slider) . ( See attachement)
 

Download TabValeursCombin.mw
 

Download TabValeursCombin.mw
 

Download TabValeursCombin.mw

Some ideas ? 

Thanks

Say I want to use multiples commands on an expression

simplify(expand(eval(expr,x=y+1)))

Now similar I could use

(simplify@expand)(expr)

but how do I do it with sth like eval which needs extra arguments

(simplify@expand@eval(x->y, x=y+1))(f)

??

Hello everybody . I have a problem with plotting a trig function : arctan(x)

I would live to plot 

In the plot ,the angle jump from -pi/2(-1,57) to pi/2 (1,57)and for u-> there is an asymptot in delta=0

But instead this the asymptote that I would like to see is in delta= -3pi.

How can I remove the discontinuity of the arctan function and plot continuos negative angle?

Thank's and best regards

t := 5;
                               5
u := exp((1/2)*x-2*t*(1/3));
                            /1     10\
                         exp|- x - --|
                            \2     3 /
plot(u(x), x = -4 .. 4);

this is the exact solution graph

while below is the approximate solution graph

h := -1;
                               -1
t := 5;
                               5
u := exp((1/2)*x)*(1+t*(27*h^4+108*h^3+240*h^2+256*h)/2^7+t^2*(27*h^4+24*h^3+168*h^2)/2^7+(1/3)*t^3*(9*h^4+16*h^3)/2^6+(1/4)*t^4*h^4/2^4);
                         6877    /1  \
                         ---- exp|- x|
                         192     \2  /
plot(u(x), x = -4 .. 4);

with(Statistics):

Variance(Distribution(ProbabilityFunction = (k -> piecewise(k > 0, 2^(-k)))));
Error, invalid input: Statistics:-Variance expects its 1st argument, X, to be of type
{algebraic, {array, list, rtable, DataFrame, DataSeries}}, but received _m1085629680096

dd := Distribution(ProbabilityFunction = (k -> piecewise(k > 0, 2^(-k))));

Variance(dd);
                               6

The first input works in neither 2D nor 1D input mode.

Another issue is that the value of the variance is in fact wrong; Maple takes the probability function to have value 1 at zero.