Why won't simplify give the obvious result and Why does testeq give a Fail result in this case?

with(RealDomain);

assume(0<f,0<g);

simplify(f*sqrt(g/f));

testeq(f*sqrt(g/f)=sqrt(f)*sqrt(g))

Hi everyone,

I have 4 expression, each containing several parameters (3, I think).

I need to find the range of the values of the parameters for which one of the expressions is at least as great as any of the other expressions.

So the input of the program/function would be the expression and which one I want to be the (weakly) greatest, and the  output should be something like: [parameter a is a member of the set [0.1,0.5] and parameter b is a member of the set [1.2, 2.4] and parameters c is negative] (meaning if these are true, then the given exoression is the (weakly) greatest of all the given expressions.

Is there a way to do this in Maple?

(For the curious: I have a game-theoretic problem, where the payoff functions include several parameters. I am trying to find the best response functions to find the Nash.equilibria of the game for the different range of parameters)

Thank you in advance,

JTamas

I have a formula that essentially is a m-order sum over a totally ordered set(the integers).

sum_{n_1=1)^{n_0} ... m times ... sum_{n_m=1}^{n_{m-1}}

It may be a little confusing but essentially it is just summing over all monotonic lidyd of size m of [1..n_0]

e.g., if n_0 = 4 and m = 3 then

1 2 3

1 2 4

1 3 4

2 3 4

Does maple have a convient way to do this? Ideally some function like osum(..., i=1..n_0, m)

where i is a multi-index(just the ordered list here so I could access i[1] the smallest, value, i[2] is the next largest, etc).

How can i solve below equation?

diff(v,x,x)+diff(v,x)*cot(x)-v(cot²(x)+nu)=a²Q/D

that Q is function of x

Q=Q[x]

Download SC-plots.mw

I am having difficulty getting solutions to a pair of PDEs.  Would anyone like to cast an eye over the attached file, incase I am missing something.

Thanks

Melvin

Hi everybody, I’m Helena

I have a function with the following general characteristics:

y= const*M*[Int_def (from a to b) [f(x)*g(x, M)*p(x, t, M,K)]dx]

::::::::::::::::::::::::::::::::::::::::::::::::::

y= dependent variable

a= 1.76E-7

b= 1.76E-9

x= integration variable

M= unknown parameter

K= unknown parameter

t= independent variable

::::::::::::::::::::::::::::::::::::::::::::::::::

... and I'm trying to fit it to my experimental data (421 values), to determine the best values for M and K. In order to do that I’m using the package Statistic and the NonlinearFit command from maple, however, it seems that NonlinearFit cannot fit the model (y(t)) to the experimental data. It gives an error: “Error, (in Statistics:-NonlinearFit) integration range or variable must be specified in the second argument, got HFloat(1.0) = 0 .. 1” for which I could not find any explanation or solution on the internet. Could you please give me some advice or ideas about how can I solve this problem?

Thank you so much in advance for your attention and help.

Hi, I would like to represent the piecewise function as a Bar Graph. I can only get a line Graph. Thanks in advance.
 

Download Ch1_3_Functions_Exam8.mw

 Recently I often drew some graphs in Graph Theory  by Maple. I tried to use DrawGraph, but  I found that evey edge only meets straight line  style in this package.

For example: 

How do I draw some edges with curve syle in graph? like following:

Thanks!

In answer to a previous question (https://www.mapleprimes.com/questions/228965-How-Can-I-Save-A-Record) acer introduced me to the .mla archives as a possible store location for variables.

I have been working with this and now end up having 4 .mla files (called results1.mla...results4.mla) that each hold output from some lengthy batch running of maple. Each of these has two tables with values called Beams and Beams2 (same name in each of these archives). While named the same the content is different.

My question is: How do I get at each of these? As acer explained, the content (i.e. the tables Beams and Beams2) is mapped into the namespace of the Maple session once libname includes the directory where the .mla files sit. They are actually in the same directory as the associated Maple programs so they mask each other.

I tried to juggle libname prepending the .mla file I want to the rest; but that does not seem to work (I always get the same data). I do need the data all in one worksheet for collating and postprocessing and display.

TIA,

M.D.

Hi everybody:

How can I solve this system of nonlinear equations without use the fsolve command?

eq[1] := -3*c[0]+3*c[1] = c[0]*(1-259*d[0]*(1/192)+43*d[1]*(1/64)-11*d[2]*(1/64)+(1/64)*d[3])-23/6;
eq[2] := -3*d[0]+3*d[1] = d[0]*(-2+23*c[0]*(1/80)+139*c[1]*(1/320)-17*c[2]*(1/160)+3*c[3]*(1/320))-1;
eq[3] := -4*c[0]*(1/3)+c[2]+(1/3)*c[3] = (8*c[0]*(1/27)+4*c[1]*(1/9)+2*c[2]*(1/9)+(1/27)*c[3])*(1-2657*d[0]*(1/5184)-343*d[1]*(1/1728)-185*d[2]*(1/1728)-79*d[3]*(1/5184))-3;
eq[4] := -4*d[0]*(1/3)+d[2]+(1/3)*d[3] = (8*d[0]*(1/27)+4*d[1]*(1/9)+2*d[2]*(1/9)+(1/27)*d[3])*(-2+109*c[0]*(1/8640)+553*c[1]*(1/1440)+559*c[2]*(1/2880)+37*c[3]*(1/1080))-5/6;
eq[5] := -(1/3)*c[0]-c[1]+4*c[3]*(1/3) = ((1/27)*c[0]+2*c[1]*(1/9)+4*c[2]*(1/9)+8*c[3]*(1/27))*(1-673*d[0]*(1/5184)-455*d[1]*(1/1728)-505*d[2]*(1/1728)-767*d[3]*(1/5184))-49/27;
eq[6] := -(1/3)*d[0]-d[1]+4*d[3]*(1/3) = ((1/27)*d[0]+2*d[1]*(1/9)+4*d[2]*(1/9)+8*d[3]*(1/27))*(-2-173*c[0]*(1/4320)+241*c[1]*(1/2880)+59*c[2]*(1/180)+2191*c[3]*(1/8640))-11/36;
eq[7] := c[0] = 1;
eq[8] := d[0] = 0;

tnx...

I'm trying to Graph a Piecewise Defined Function.

j(x) = piecewise(x < 0, x^2 + 1, 0 <= x, x - 1);

It doesn't provide a plot option.

Download piecewisegraph.mw

Download funct_misrep.mw

data := [[0., 9.1300, 0.931e-1, 0.899e-1, .1000, 0.], [30.0000, 8.9300, .1270, .1230, .2270, 0.49e-2], [60.0000, 8.6000, .1510, .1390, .4920, 0.153e-1], [90.0000, 8.2800, .1540, .1490, .7780, 0.249e-1], [120.0000, 7.9700, .1540, .1570, 1.0700, 0.329e-1], [150.0000, 7.8600, .1540, .1600, 1.1700, 0.348e-1], [180.0000, 7.8100, .1530, .1530, 1.2100, 0.404e-1], [210.0000, 7.7700, .1400, .1420, 1.2800, 0.432e-1]];
des := [diff(y1(t), t) = -k1*y1(t)-k2*y1(t), diff(y2(t), t) = k2*y1(t)-k3*y2(t), diff(y3(t), t) = k1*y1(t)+k3*y2(t)-k4*y3(t), diff(y4(t), t) = k4*y3(t)-k5*y2(t)*y4(t)+k6*y5(t), diff(y5(t), t) = k5*y2(t)*y4(t)-k6*y5(t)];
ics := seq((y || i)(0) = data[1, i+1], i = 1 .. 5):
Error, unable to match delimiters
Typesetting:-mambiguous(Typesetting:-mambiguous(ics Assign seq

  lparlpary verbarverbar irpar(0) equals data(1comma i + 1)comma 

  i equals 1 periodperiod 5rparcolon, 

  Typesetting:-merror("unable to match delimiters")))
res := dsolve({ics, des[]}, numeric, parameters = [k1, k2, k3, k4, k5, k6]);
Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations
timeList := [0, 30, 60, 90, 120, 150, 180, 210];
              [0, 30, 60, 90, 120, 150, 180, 210]
sse := proc (k1, k2, k3, k4, k5, k6) res(parameters = [k1, k2, k3, k4, k5, k6]); add((rhs(select(has, res(timeList[i]), y1)[])-data[i, 2])^2+(rhs(select(has, res(timeList[i]), y2)[])-data[i, 3])^2+(rhs(select(has, res(timeList[i]), y3)[])-data[i, 4])^2+(rhs(select(has, res(timeList[i]), y4)[])-data[i, 5])^2+(rhs(select(has, res(timeList[i]), y5)[])-data[i, 6])^2, i = 2 .. 8) end proc;
proc(k1, k2, k3, k4, k5, k6)  ...  end;
c := GlobalOptimization:-GlobalSolve(('sse')(k1, k2, k3, k4, k5, k6), k1 = 0 .. 1, k2 = 0 .. 1, k3 = 0 .. 1, k4 = 0 .. 1, k5 = 0 .. 1, k6 = 0 .. 1, timelimit = 10); [.219132447080011505, [k1 = 0.852482740113834e-3, k2 = 0.52683998680924474e-4, k3 = 0., k4 = 0.5113239298267808e-1, k5 = 0.4363021255887466e-2, k6 = 0.]];
Error, `GlobalOptimization` does not evaluate to a module
[0.219132447080011505, [k1 = 0.000852482740113834, 

  k2 = 0.000052683998680924474, k3 = 0., 

  k4 = 0.05113239298267808, k5 = 0.004363021255887466, k6 = 0.]]

res(parameters = c[2]):  p:=Array(1..5):  for n from 1 to 5 do     p[n]:=plots:-display(plots:-odeplot(res, [t, (y || n)(t)], t = 0 .. 210),plots:-pointplot([seq([data[i, 1], data[i, n+1]], i = 1 .. 8)]));  end do; plots:-display(p) ;
Error, unable to match delimiters
Typesetting:-mambiguous(Typesetting:-mambiguous(res(parameters 

  equals c(2))colon  pAssignArray(1periodperiod5)colon  for n 

  from 1 to 5 do     p(n)Assignplotscolon - displaylparplotscolon

   - odeplotlparrescomma lsqbtcomma lpary verbarverbar nrpar(t)

  rsqbcomma t equals 0 periodperiod 210rparcommaplotscolon - 

  pointplot((seq((data(icomma 1)comma data(icomma n + 1))comma i 

  equals 1 periodperiod 8)))rparsemi  end dosemi plotscolon - 

  display(p) , Typesetting:-merror("unable to match delimiters")))

 I cannot solve this equation. Please help me

How can I use maple to get the numerical solution of this non-autonomous system

sys_ode := diff(x(t), t) = ax(t)-bx(t)y(t)+2sin(t), diff(y(t), t) = -cy(t)+dx(t)y(t)

where a,b,c,d is parameters, they can take any value.