Hello people in mapleprimes,

I have a question about how Int does.

The following function spy returns 0, of course, with a side effect of listing the value of x one by one to secrets.

secrets := NULL:

spy:=proc(x::{name,numeric})
  global secrets;
  if type(x,name) then
    return 'procname'(args)
  else
  secrets:=x,secrets;
  return 0;
  end if;
end proc;

And, with this function, calculation of the Int, that is, following brings a sequence of numbers:

evalf(Int(spy,0..1));

secrets;

.7506605773, .2493394227, .9118140517, 0.881859483e-1, .9970470440, 0.29529560e-2, 1.0000000000, 2.2449529449*10^(-11), .5000000000

The question I have is why the number of this sequence is not from smaller( or greater) to greater (smaller) in order,

but in random order. And, numerical calculation of Int can be done with only 9 points extracted?

Best wishes.

taro

This post is my attempt to answer the question from here .  

The procedure  ContoursWithLabels  has 2 required parameters: Expr  is an expression in  x  and  y  variables,  Range1  and  Range2  are ranges for  x  and  y . In this case, the output is the list of floats for the contours and 8 black contours (with labels) (the axis of coordinates as a box). 

The optional parameters: Number is positive integer - the number of contours (by default Number=8),  S is a set of real numbers  C  for contours (for which Expr=C) (by default  S={}),  GraphicOptions  is a list of graphic options for plotting (by default  GraphicOptions=[color = black, axes = box]),  Coloring  is an equality  Coloring=list of color options for  plots[dencityplot]  command (by default Coloring=NULL). 

The code of the procedure:

restart;

ContoursWithLabels := proc (Expr, Range1::(range(realcons)), Range2::(range(realcons)), Number::posint := 8, S::(set(realcons)) := {}, GraphicOptions::list := [color = black, axes = box], Coloring::`=` := NULL)

local r1, r2, L, f, L1, h, S1, P, P1, r, M, C, T, p, p1, m, n, A, B, E;

uses plots, plottools;

f := unapply(Expr, x, y);

if S = {} then r1 := rand(convert(Range1, float)); r2 := rand(convert(Range2, float));

L := [seq([r1(), r2()], i = 1 .. 205)];

L1 := convert(sort(select(a->type(a, realcons), [seq(f(op(t)), t = L)]), (a, b) ->is(abs(a) < abs(b))), set);

h := (L1[-6]-L1[1])/Number;

S1 := [seq(L1[1]+(1/2)*h+h*(n-1), n = 1 .. Number)] else

S1 := convert(S, list)  fi;

print(Contours = evalf[2](S1));

r := k->rand(20 .. k-20); M := []; T := [];

for C in S1 do

P := implicitplot(Expr = C, x = Range1, y = Range2, op(GraphicOptions), gridrefine = 3);

P1 := [getdata(P)];

for p in P1 do

p1 := convert(p[3], listlist); n := nops(p1);

if n < 500 then m := `if`(40 < n, (r(n))(), round((1/2)*n)); M := `if`(40 < n, [op(M), p1[1 .. m-11], p1[m+11 .. n]], [op(M), p1]); T := [op(T), [op(p1[m]), evalf[2](C)]] else

if 500 <= n then h := floor((1/2)*n); m := (r(h))(); M := [op(M), p1[1 .. m-11], p1[m+11 .. m+h-11], p1[m+h+11 .. n]]; T := [op(T), [op(p1[m]), evalf[2](C)], [op(p1[m+h]), evalf[2](C)]]

fi; fi; od; od;

A := plot(M, op(GraphicOptions));

B := plots:-textplot(T);

if Coloring = NULL then E := NULL else E := ([plots:-densityplot])(Expr, x = Range1, y = Range2, op(rhs(Coloring)))  fi;

display(E, A, B);

end proc:

Examples of use:

ContoursWithLabels(x^2+y^2, -3 .. 3, -3 .. 3);

ContoursWithLabels(x^2-y^2, -5 .. 5, -5 .. 5, {-20, -15, -10, -5, 0, 5, 10, 15, 20}, [color = black, thickness = 2, axes = box], Coloring = [colorstyle = HUE, colorscheme = ["White", "Red"], style = surface]);

The next example, I took from here:

ContoursWithLabels(sin(1.3*x)*cos(.9*y)+cos(.8*x)*sin(1.9*y)+cos(.2*x*y), -5 .. 0, 2 .. 5, {seq(-2 .. 2, 0.5)}, [color = black, axes = box], Coloring = [colorstyle = HUE, colorscheme = ["Cyan", "Red"], style = surface]);

There are many more examples can be found in the attached file. 

ContoursWithLabels1.mw

Edit. The attached file has been corrected.

Hi,

I recently noticed that Maple 2015 has become irresponsive on mac (macbook pro) with the latest java release (Java 8 update 66). Did anywone experience the same problem ?

I have functions like

g := (y, t)->sin(y(t))

diff(g(y, t), t) # ok

How to determine the derivative below with maple :

diff(g(y, t), y(t))

Thanks

Iam a newbie, just two weeks into my 30-days trial. I have been exploring the symmetry aspect of PDEtools, gone through materials in the help section but still having problem in some of my analysis. The answer to the titled question "Symmetry analysis with parameters" was really helpful but did not work out for me when the parameters are more than one. Attached is a sample question.

sample_question.mw

hi.for linear coupling equation

dsys3 := {-72.49829200*(diff(f1(x), x, x))+0.8377580411e-2*(diff(f2(x), x))-8.873545400*10^9*(diff(f3(x), x))+2.114533515*10^18*f1(x), -878.8477313*(diff(f2(x), x, x))+1.590065471*10^20*f2(x)-7.353421206*10^(-26)*(diff(f3(x), x, x))+4.891459762*10^10*f3(x), 4.027667395*10^(-20)*(diff(f3(x), x, x, x, x))-0.6274394007e-2*(diff(f3(x), x, x))+8.873545401*10^9*(diff(f1(x), x))-7.352113720*10^(-26)*(diff(f2(x), x, x))+4.904509456*10^10*f2(x)+1.208381068*10^19*f3(x)-2.499990383*10^26*omega*(diff(f3(x), x, x))}

by assuming and to expand these functions(f1-f2-f3) in polynomial form(e.g. Chebyshev, power polynomials, Legendre and etc).:

f1(x):=(&sum;)H[i] *(e)^(lambda[i] *x); f2(x):=(&sum;)alpha[i]*H[i] *(e)^(lambda[i] *x);f3(x):=(&sum;)GAMMA[i]*H[i] *(e)^(lambda[i] *x)

how i detemine lambda[i] which are roots of the  characteristic equations?in other word how i can build characteristic relation for coupling equation?

2)how i can gain value for alpha[i] and GAMMA[i]

by using equation

Q1 := subs(x = 0, sum(H[i]*exp(lambda[i]*x), i = 1 .. 8)); Q2 := subs(x = L, sum(H[i]*exp(lambda[i]*x), i = 1 .. 8)); Q3 := subs(x = 0, sum(alpha[i]*H[i]*exp(lambda[i]*x), i = 1 .. 8)); Q4 := subs(x = L, sum(alpha[i]*H[i]*exp(lambda[i]*x), i = 1 .. 8)); Q5 := subs(x = 0, sum(GAMMA[i]*H[i]*exp(lambda[i]*x), i = 1 .. 8)); Q6 := subs(x = L, sum(GAMMA[i]*H[i]*exp(lambda[i]*x), i = 1 .. 8)); M := diff(sum(GAMMA[i]*H[i]*exp(lambda[i]*x), i = 1 .. 8), x); Q7 := subs(x = 0, M); Q8 := subs(x = L, M)

????

thanks...

chebyshev.mw

When executing the following orders,

a:=0.006:

convert(a, symbol);

the result is .6e-2.

How to make Maple return the form 0.006, just as the original?

I am trying to solve 4 nonlinear equations for four variables using fsolve  and the output that i am getting is basically the same equations repeated after some time.  I even tried reducing one of the equations using assumptions from my side but it results in same behaviour..  Quite new to maple, would like some advice as to this behaviour. Thanks

 Here's the file

fsolve_1.mw

PS- using do loop is part of the solving so i cannot remove that

plotpoints.mw

I want to plot points when i_t =1,2,3,..,11,12 instead of a continous line displayed in the worksheet I uploaded. How to modify the function? Thank you for helping:)

Occasionally I have need to import 1D Maple expressions into Mathematica, because I believe mma has superior simplify functionality (sorry Maple). Mathematica has its own syntax which (of course) differs to maple.  eg sin(x) in maple is Sin[x] in mma ; ln(x)=Log[x] ;....... I'm aware of the package MmaTranslator, but I think it only converts mma to maple.

so i need some code pls.

(i have bolded the functions which need to be converted)

eg: i/p:

(1/365)*(365*B21*EP*ln(1+i)*(1+i)^B21*hr*kw*p-365*EP*hr*kw*p*(1+i)^(-1+B21)*i-365*EP*hr*kw*p*(1+i)^(-1+B21)+365*EP*hr*kw*p+SC*i)/(hr*kw*(EP*ln(1+i)*(1+i)^B21*p-FIT*i*p+cos(omega)))

o/p:

(1/365)*(365*B21*EP*Log[1+i]*(1+i)^B21*hr*kw*p-365*EP*hr*kw*p*(1+i)^(-1+B21)*i-365*EP*hr*kw*p*(1+i)^(-1+B21)+365*EP*hr*kw*p+SC*i)/(hr*kw*(EP*Log[1+i]*(1+i)^B21*p-FIT*i*p+Cos[omega]))

Hi everyone

I've come across a RAM-issue.

My program creates a lot of plots using millions of variables-values and it eats up my computers memory since I use the plots to make an animation.

The plots are simple 3d plots created using the surfdata function.

How do I save plots into a single file as the program creates them concurrently such that after the program has terminated I can use the plots to create an animation?

and how do I efficiently delete variables I dont not need anymore?

Dear Friends

            Hope everything going fine with you. I want the numerical solution of nonlinear system of ordinary differential equations using RK method. The system of ODEs and their required results are present in attached file. I am waiting your quick response.

Help.docx

With my best regards and sincerely.

Mob #: 0086-13001903838

Does [Length of output exceeds limit of 1000000] mean that the function couldn't be solved and I have to modify the mathematical model?
If the results are just too long to be displayed, how can I see the results?

how to make physical quantities graphs from RK45 ode bvp convective boundary layer flow rather than using lines

i use lines command copy one by one yje data from the solution

are they any simplest way to sketch graphs for physical quantities like f'', thetha prime versus 3 parameters?

I try to find the exact (symbolic) value of

(-2*sqrt(7)-4)*EllipticK((1/8)*sqrt(2)*(-3+sqrt(7)))^2+4*EllipticE(-(1/8)*sqrt(2)*(-3+sqrt(7)))*sqrt(7)*EllipticK((1/8)*sqrt(2)*(-3+sqrt(7)))

I tried 'simplify' with different options and 'convert'. It would be pi=3.141... as numerical approximation suggests.

Many thanks.