In Maple V the display of a lengthy expression often contained "auxiliary variables" named %1, %2, ... to make it more readable ( (a personal point of view, see for instance Introduction to Maple, A. Heck, p 97).
When did this feature disappear? Has it been replaced by something else?

TIA

 

how this integration can be computed and plotted? tnx in advance

 

Download code.mw

I have the following equations and plot in Maple, but I would like for Maple to only consider positive values of k, but I cannot find a way to do this as k is dependent on x. Is there any way to do this either in the code for k, or in the code for the plot?

For an expression like , I want to extract the subscripts in an another set/list [{1,1}, {2,3}, {1,1,2}]. Any help is much appreciated. Thanks

Dear maple users,

A fine day wishes to all.

Here, we have computed the fN(x,t) value by using pdsolve.

We have to integrate the computed value and need to find the values with the sequence of x.

A1:=int(fN(x,1.12),x) 

A2:=seq(A1,x=0..1,0.1)

How to integrate the computed values 

And, How to find the values for the sequence of x.

JVB1.mw
 

Download JVB1.mw

I would like to show : in a quadrilateral circumscribed to an ellipse, the line passing through the middle of the diagonals passes through the centre of the ellipse.
My code is :

restart; with(geometry): with(plots): `local`(O):
_EnvHorizontalName := x: _EnvVerticalName := y:

alias(coor = coordinates):
ell := x^2/a^2+y^2/b^2 = 1:
point(P1,a*cos(omega), b*sin(omega)):
point(P2,a*cos(omega-(1/2)*Pi), b*sin(omega-(1/2)*Pi)): 
point(P3,a*cos(omega+(8/7)*Pi), b*sin(omega+(8/7)*Pi)):
point(P4,a*cos(omega+5*Pi*(1/2)), b*sin(omega+5*Pi*(1/2))):
a := 5: b := 3: omega := (1/5)*Pi:
Ell := implicitplot(ell, x = -a .. a, y = -b .. b, color = red): 
dr := draw([seq(P || k, k = 1 .. 4)], axes = normal, printtext = true):

for i from 1 to 4 do tgP||i := x*coor(P||i)[1]/a^2+y*coor(P||i)[2]/b^2 = 1 od:
poly := Matrix([coor(P1), coor(P2), coor(P3), coor(P4)]):
Quadri := polygonplot(poly, axes = normal, color = "DarkGreen", transparency = .8):

with(combinat): with(ListTools):
L := [1, 2, 3, 4]:
for i from 1 to 4 do Rotate(L, i)[1] od:
for i to 4 do solve({(tgP || Rotate)(L, i)[1], tgP || i}, {x, y}); point(S || i, subs(%, x), subs(%, y)); coor(S || i) end do;
Error, invalid input: subs received 1, which is not valid for its 1st argument
#otherwise
solve({tgP1, tgP2}, {x, y}): point(S1, subs(%, x), subs(%, y)); coor(S1):
                               S1
solve({tgP2, tgP3}, {x, y}): point(S2, subs(%, x), subs(%, y)); coor(S2):
                               S2
solve({tgP3, tgP4}, {x, y}): point(S3, subs(%, x), subs(%, y)); coor(S3):
                               S3
solve({tgP1, tgP4}, {x, y}): point(S4, subs(%, x), subs(%, y)); coor(S4):
                               S4

poly := Matrix([coor(S1), coor(S2), coor(S3), coor(S4)]):
Quadri2 := polygonplot(poly, axes = normal, color = "DarkGreen", transparency = .9):
#dr2:=draw(seq(S||k,k =1..4), axes = normal, printtext = true):
line(diag13, [S1, S3]): line(diag24, [S2, S4]): midpoint(M1, S1, S3): midpoint(M2, S4, S2): 
line(Lm, [M1, M2]):
dr2 := draw([S1, S2, S3, S4, M1, M2, Lm(color = black), diag13, diag24], axes = normal, printtext = true):
for i from 1 to 4 do
TgP||i := implicitplot(tgP||i, x = -a-5 .. a+5, y = -b-5 .. b+5, color = blue) od:
display([Ell, seq(TgP||i,i=1..4), Quadri, Quadri2,dr,dr2], view = [-a-5 .. a+3, -b-2 .. b+2], 
scaling = constrained, size = [700, 700]); Thank you for your answere.

How create a rtable with a global indexing variable outside the procedures and use it inside the procedure's the same rttable with index as global

RA := rtable(1 .. 5, 1 .. 5)

ik := 1

s := proc() global ik; RA[ik, 1] := 1; ik := ik + 1; RA[ik, 2] := 2; ik := ik + 1; end proc;

t := proc() global ik; RA[ik, 3] := 3; ik := ik + 1; RA[ik, 4] := 4; ik := ik + 1; end proc;

I dont want pass the RA to the procedure as a parameter I want declare the table RA as global and it to useable inside every procedure how to did if I have lot procedure and procedure within a procedure calls

B(4) := x^4 - 4*x^3/(2 + p) - 6*(p - 1)*x^2/((3 + p)*(2 + p)) - 4*p*(p - 5)*x/((4 + p)*(3 + p)*(2 + p)) - (p - 1)*(p^2 - 15*p - 4)/((5 + p)*(4 + p)*(3 + p)*(2 + p));
for p from -1 to 5000 do
    A(p) := fsolve(B(4), x, ':-complex');
end do;
Error, (in fsolve) p is in the equation, and is not solved for
ptlist := [0];
for j from -1 to 5000 do
    ptlist := [op(ptlist), A(j)];
end do;
with(plots);
complexplot(ptlist, x = -1 .. 1.5, y = -0.5 .. 0.5, style = point);

I get an error with this code, but when I replace B(4) with its assignment in fsovle, it works. Why is this, and is there a way I can write B(4) instead of the longer expression x^4 - 4*x^3/(2 + p) - 6*(p - 1)*x^2/((3 + p)*(2 + p)) - 4*p*(p - 5)*x/((4 + p)*(3 + p)*(2 + p)) - (p - 1)*(p^2 - 15*p - 4)/((5 + p)*(4 + p)*(3 + p)*(2 + p)) in the fsolve argument?

ishak.mw
 

Download ishak.mw

Good day every one;

please im soliciting for a help on how to plot my Nusselt number and Skin friction.

The attached is plotting against the dependent variable (eta) but i want skin friction (f  ' ' ) against N or Pr not against eta 

thank you for your help in aticipation

I have the following polynomials, and I want to multiply each term of the polynomial by the least common denomiator, so that the denominators are cleared. However, I want to keep the products and powers of terms together. So for example, B(2) should become (2+p)(3+p)x^2-2(3+p)x-p-1.

How can I accomplish this? Alternatively, given the following polynomials, I want to factor only the expressions of p in the brackets, so we have (2+p), (2+p)(3+p), etc.

Hi, 

I'm trying to solve an ODE. I don't understand why Maple returns nothing when I run this code. 

with(ODETools):
dsolve(diff(y(x), x) = (6*y(x)^5 - 3*y(x)*x^2 - 20*y(x)^3*x)/(-4*x^3 + 30*y(x)^2*x^2 - 30*y(x)^4 + 7*y(x)^6), y(x));

I don't have any error messages, it just doesn't return anything.

Thank you very much for your time and your help. 

Hello. Plotting a graph in the polar coordinate system using polar plot (example). Is it possible to make it display 90 degrees instead of Pi/2, 45 degrees instead of Pi/4, and so on. Thanks

Hello,

How to solve an equation like and show all values ?

Thanks

Whenever i press F5, Maple will switch to nonexecutable math so I have to press it once more to go to math. This is very annoying. I would very much like to have it just switch between math and text when i press F5 like it used to do. 
How do I make it so?

Please help.

I didn't realize until now that ExcelTools:-Import will not return a value for a cell that has an equation inside. Am I doing something wrong?