Hi,

If I use a package  that I don't wnt to load, is these two is the best way to proceedcontains

• alias (P= package): 
• macro (P= package):
   

Thanks in advance

restart;
with(plots); with(LinearAlgebra);
_EnvHorizontalName := 'x';

_EnvVerticalName := 'y';

x1,y1,x2,y2,x3,y3:=0,-3,3,1,5,-2:   
A := [x1, y1]: B := [x2, y2]: C := [x3, y3]:

Barycentre := proc (A, B, t) description "Barycentre de 2 points A(1) et B(t) dans le rapport t";
return [(1-t)*A[1]+t*B[1], (1-t)*A[2]+t*B[2]] end proc;
ellip := proc (r1, r2) local a, b, c, d, e, f, D, E, F, eq1, eq2, eq3, eq4, eq5, eq6, x0, y0, EE, r3, sol, Ff, Tg;
global A, B, C;
r3 := -1/(r2*r1);
D := Barycentre(C, B, 1/(1-r1)); E := Barycentre(A, C, 1/(1-r2)); F := Barycentre(B, A, 1/(1-r3));
Ff := proc (x, y) options operator, arrow; a*x^2+2*b*x*y+c*y^2+2*d*x+2*e*y+f end proc;
Tg := proc (x0, y0, x, y) options operator, arrow; a*x*x0+b*(x*y0+y*x0)+c*y*y0+d*(x+x0)+e*(y+y0)+f end proc;
eq1 := Ff(D[1], D[2]);
eq2 := Ff(E[1], E[2]);
eq3 := Ff(F[1], F[2]);
eq4 := Tg(F[1], F[2], x1, y1);
eq5 := Tg(D[1], D[2], x2, y2);
eq6 := Tg(E[1], E[2], x3, y3);
sol := op(solve([eq1, eq2, eq3, eq4, eq5, eq6], [a, b, c, d, e]));
assign(sol);
EE := subs(f = 1, Ff(x, y) = 0) end proc;

ellip(-1, -7); tri := plot([A, B, C, A], color = blue):
 
po := plot([A, B, C], style = point, symbolsize = 15, symbol = solidcircle, color = red);
tp := textplot([[A[], "A"], [B[], "B"], [C[], "C"]], 'align' = {'above', 'left'});
x := 'x'; y := 'y';
ELL := seq(implicitplot(ellip(-7/11-(1/11)*j, -1/17-3*j*(1/17)), x = 0 .. 5, y = -3 .. 1, color = ColorTools:-Color([rand()/10^12, rand()/10^12, rand()/10^12])), j = 1 .. 17);
display([tri, ELL, po, tp], view = [-.5 .. 5.5, -4 .. 1.5], axes = none, scaling = constrained, size = [500, 500]);
Explore(implicitplot(ellip(r1, r2), x = 0 .. 5, y = -3 .. 1), parameters = [r1 = -2.18 .. -.7, r2 = -3 .. -.23]);
Can you tell me why this last instruction does't work ? Thank you.
 

Hi guys ,

I have the equation dt=(L/r)*(1 - a^3/r^3)^(-1/2) dr which want to integrate on both side and then solve in terms of r.(L and a are constants). I know the answer (r=a (cosh(3t/2*L)^2/3),but it seems maple can not produce it!

Mathematical almost compute it corectly! altough i think  tanh-1, should be cosh^-1(x/a)^3/2

I would appreciate if someone can help me

Thanks so muchproblem.mw

Why does Maple write

restart;
eq:=x-infinity=0;

as 

Hi sir, I hope everyone are good in Covid situation.

I wish to obtain the dual branch solution as given in the Article.

Please help me to obtain the dual branch solution in Maple.

The article is here. akbar2014_(1).pdf

Thank you..

Hi

How can i get y value from a x value to show in a graph plot?. Like i can with the point probe, but i would like to mark it with a line.

Something like this (this is from a solution that isn't mine):

Where i am supposed to find y values to specific time values..

I rewrite my code with Grid library using local vCPU of amazon 

discover no license of distributed HPC when setup command show need go to acresso.com

can Grid still be used and function with local 96 number of vCPU?

then when I test code, I can not pkill mserver in ps -aux in LInux 

can the performance really improved ? Because I suppose 3 to 5 minutes mserver will end and disappear from ps -aux but grid node of mserver still running

originally I can run 100 batches every day., but I have to monitor decrease of memory in order to determine whether can kill mserver for next batch 

I expect run 400 batches per day. But Can not kill mserver when using grid

In my code I had using time limit(30, ...) 

when using Grid seq of function , can lprint work normally to get results into text ?

i notice Grid up to 57, do I need to recalculate and revise code to fit 96 vCPU for grid number 96?

Hello;I want to solve the following system.the result is not what I need.

with(linalg):

f1:=b*R + b*S - d*S - bcd*I1*S - bcd*I2*S - (e*g*S)/K - (g*I1*S)/K - ( g*I2*S)/K - (g*R*S)/K - (g*S^2)/K=0;

f2:=-d*e - (e^2*g)/K - (e*g*I1)/K - (e*g*I2)/K - (e*g*R)/K + bcd*I1*S +  bcd*I2*S - (e*g*S)/K - e*scd1=0;

f3:=-d*I1 - (e*g*I1)/K - (g*I1^2)/K - (g*I1*I2)/K - (g*I1*R)/K - ( g*I1*S)/K + e*scd1 - e*f*scd1 - I1*scd2=0;

f4:=-acd*I2 - d*I2 - (e*g*I2)/K - (g*I1*I2)/K - (g*I2^2)/K - ( g*I2*R)/K - (g*I2*S)/K + e*f*scd1=0;

f5:=-d*R - (e*g*R)/K - (g*I1*R)/K - (g*I2*R)/K - (g*R^2)/K - (g*R*S)/K +  I1*scd2=0;

solve({f1,f2,f3,f4,f5},{S,e,I1,I2,R});

What Maple function answers the following:

Produce (a minimal set of) sets that generate, by union, every intersection of sets given.

Examples:
F({1, 2, 3}, {1, 4, 5}, {2, 3, 4});
{}, {1}, {4}, {5}, {2, 3};

F({1, 2, 3, 6}, {1, 4, 5, 6}, {2, 3, 4, 6});
{1}, {4}, {5}, {6}, {2, 3};

Correction:
Produce (a minimal set of) disjoint sets that generate, by union, every intersection of sets given.

Intersection and union here are n-ary operators and every combination of sets is considered for intersection.

I'd like to delete my mapleprime account and purge all associated data. How to achieve this?

I have attached a maplework please check it. I want to solution graph by using RK 4th order with shooting technique or RK Fehlberg method. help.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/help.mw .
 

Download help.mw

Mal_simulation.mw

Sketch the 3d graph of solid bounded inside a sphere of x²+y²+z²=4z and below a cone z=√(x²+y²) on xy-plane.

Any one has smart way to help verify this Maple solution to this ODE? I can't figure it out

restart;
ode:=diff(y(x),x)-y(x) = x*y(x)^(1/2):
ic:= y(0)=4:
maple_sol:=dsolve([ode,ic],'implicit'):
the_diff_1:=odetest( (lhs-rhs)(maple_sol)=0,[ode,ic]);


maple_sol:=dsolve([ode,ic],'explicit'):
the_diff_2:=odetest( (lhs-rhs)(maple_sol)=0,[ode,ic]);

Tried few simplifications but not everything. I am assuming Maple solution is correct.

Maple 2021.1

Hello there, 

While I was trying to calculate the complex (i.e., P and Q) power equations, I got this question. 

Thus, I wonder if I can have a chance to ask the question here:

a. When the voltage and current are defined as functions, taking the t (time) as a input variable, I was not able to get the proper outcome from the operation. Those attempts are labeled as '#1' and '#2' in the attached worksheet. 

b. When the voltage and current are defined as simple complex expression, I was able to see what I wanted. The expression is named as 'ideal' in the attached worksheet. 

c. Therefore, my question is to see if there is a way to make the operaion using the functions in the same way as done by the simple expressions. 

d. Any insight would be appreciated. 

Thank you. 

Q20210617.mw