Is it  possible to draw an ellipse inscribled in a given triangle and we know its orthoptist circle ? Thank you very much.

Etude d'un cas particulier a := 5: b := 7: k := 9: A := [a, 0]: B := [0, b]: #A et B fixes P := [t, 0]: Q := [0, k/t]:#P et Q 2 points mobiles cir := -a*x-b*y+x^2+y^2 = 0: sol := solve(subs(y = 5, cir), x): cen := [solve(diff(cir, x)), solve(diff(cir, y))]: x0 := sol[1]: y0 := 5: M := [x0, y0]: R := sqrt(cen[1]^2+cen[2]^2): beta := arctan(diff(solve(EQ(M, cen), y), x)): Recherche des valeurs de t pour que les 2 droites soient perpendiculaires eq := t^2*(y0-b)+t*(a*b-a*y0+b*x0-k)-x0*(a*b-k) = 0; sol := solve(eq, t); t := sol[1]; tp := sol[2]; P1 := [t, 0]; Q1 := [0, k/t]; PQ1 := simplify(x*(-a*b+b*t+k)+y*t*(t-a)-t*(-a*b+b*t+k)) = 0:#1ere tangente PQ2 := simplify(x*(-a*b+b*tp+k)+y*tp*(tp-a)-tp*(-a*b+b*tp+k)) = 0:#2ième tangente P2 := [tp, 0]; Q2 := [0, k/tp]; CIR := implicitplot(cir, x = -4 .. 8, y = -4 .. 12, color = red); Fig := proc (alpha) local Dr1, DR1, Dr2, DR2, N, u0, v0, Po, t, tp, sol; global a, b, k, cen, R; u0 := cen[1]+R*cos(alpha); v0 := cen[2]+R*sin(alpha); N := [u0, v0]; sol := solve(t^2*(v0-b)+t*(b*u0-a*v0+a*b-k)-u0*(a*b-k) = 0, t); t := sol[1]; tp := sol[2]; Dr1 := simplify(x*(-a*b+b*t+k)+y*t*(t-a)-t*(-a*b+b*t+k)) = 0; DR1 := implicitplot(Dr1, x = -4 .. 8, y = -4 .. 12, color = brown); Dr2 := simplify(x*(-a*b+b*tp+k)+y*tp*(tp-a)-tp*(-a*b+b*tp+k)) = 0; DR2 := implicitplot(Dr2, x = -4 .. 8, y = -4 .. 12, color = pink); Po := pointplot([N[]], symbol = solidcircle, color = [black], symbolsize = 8); display([Po, DR1, DR2]) end proc; DrPQ1 := implicitplot(PQ1, x = -4 .. 22, y = -4 .. 12, color = blue); DrPQ2 := implicitplot(PQ2, x = -4 .. 22, y = -4 .. 12, color = blue); Points := pointplot([A[], B[], M[], P1[], P2[], Q1[], Q2[], cen[]], symbol = solidcircle, color = [green], symbolsize = 10); T := plots:-textplot([[A[], "A"], [B[], "B"], [M[], "M"], [P1[], "P1"], [P2[], "P2"], [Q1[], "Q1"], [Q2[], "Q2"], [cen[], "cen"]], font = [times, 10], align = {below, left}); n := 19; display([seq(Fig(2*i*Pi/n), i = 0 .. n), Fig(beta), CIR, DrPQ1, DrPQ2, Points, T], scaling = constrained, size = [500, 500]); I would find out the focus of the ellipse. Thank you.

Please, anyone with the idea on how to plot D1 against alpha[n]. The link is attached below

Help!!!.mw

 

I am trying to use task programming but getting incorrect results and taking more time when using a single core.  Please help

TASK_PARALLEL3.mw

Dear friends, please I would like to ask for your help with the following problem: 

I need to invoke the number of elements of an Array working with parallel programming in the task programming model. I've tried to used the command rtable_num_elems as it is contained in the thread safe functions lists. However, Maple does not recognize it as I obtain the error "Bad index into array". Using the same code, I've substituted the array for a list and rtable_num_elems for nops and the code works perfectly. What could I be doing wrong? I need to use arrays given the extension of the data I'm handling. 

Many thanks for your kind help. 

Anyone have any thoughts on how I can combine these two terms? (see screen shot)  The error message implies the units are somehow not really the same but there is no help page for this error.  Any insight would be appreciated.  Thanks.

Here is the file:  units_issue.mw

 

I am trying a multistart local search over a multimodal function. I need to locate around 700 local minima where the dimension of the function is very high. It will be really helpful if the computational time can be reduced. Thank you.

 

opt.mw

example let P be a permutation: P: = Matrix ([[1,2,3], [4,5,6]])
calculate P ^ 4567 knowing that the order of P = 20

how to determine the matrix of a permutation by knowing the orbit in Maple?

Why doesn't Maple show a solution to the following odesys of second order and not show an error?
 

Differential Equations Trebuchet, Phase II_2020-06-05   (1) (2) (3) (4) (5) (6)   (7) (8) (9) (10) (11) (12) Error, (in ODEtools/info) Required a specification of the indeterminate function

 

Download Differential_Equations_Trebuchet_Phase_II_2020-06-05.mw

 

Hi,

I need taylor expasion of

A(w):=M/sqrt(M^2+2*w)-mup*sqrt(((mu*M^2+3*sigmap-2*w)-sqrt((mu*M^2+3*sigmap-2*w)^2-12*mu*sigmap*M^2))/(6*sigmap))

I used taylor(A,w,4), but I had a problem!

If I know that A(w) is a real function, how can I remove ‘csgn’ from the  ‘taylor(A,w,4)’ command:

Thanks.

12.mws

restart; with(plots):with(LinearAlgebra):unprotect(O); alias(conj = conjugate); conj z = lambda*v+a; La droite D est représentée par son équation complexe Appelons H l'affixe h le pied de la perpendicukaire abaissée de O sur (D) Les vecteurs OH et V sont orthogonaux donc z = lambda v + a h*conj(v)+conj(h)*v = 0; Le point H appartient à la droite (D) donc : h = lambda*v+a; conj(h) = conj(a)+lambda*conj(v); (lambda*v+a)*conj(v)+(conj(a)+lambda*conj(v))*v = 0; solve(%, lambda); h := simplify(subs(lambda = %, lambda*v+a)); a := 3-I*4; v := -2/3+4*I; evalc(h); H := [Re(h), Im(h)]; Représentation graphique d'un cas particulier : f := proc (x) options operator, arrow; -3*x+5 end proc: a := 3: A := [a, f(a)]:O:=[0,0]: zo := [8/3+I*f(8/3)]; ze := [2+I^(eval(diff(f(x), x), x = 2))]; Zo := [8/3, f(8/3)]; Ze := [2, f(2)]; ex := -3*x+5; V := `

Good day sirs, I write a system of DAE but giving me this code "(The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)". The code is attached below.

Thanking you in anticipating for your help.

Help!!!!.mw

 

f := x -> exp(-x)*sin(x); intvx:= 0..3;
f := proc (x) options operator, arrow; exp(-x)*sin(x) end proc
intvx := 0 .. 3
 

And this one below ( i prefer this one ) , but got in worksheet now the one above
Probably a option issue ?

 

f := x -> exp(-x)*sin(x); intvx:= 0..3;

restart;

ii := 50;

seq(ii, ii = 0 .. 5);

evalf(int(ii, ii = 0 .. 5))

 

For sequence the syntax is correct i.e the output is 0,1,2,3,4,5. is ii inside sequence command takes the value 5 assigned it then changes when ii ranges from 0 to 5?

restart;

ii:=50:
seq(ii)

output: 50

restart;

ii:=50:

seq(seq(ii),ii=1..5)

expected answer: 50,50,50,50,50

obtained: 1,2,3,4,5

please explain how seq command works