Hi everyone.

I have a 2D function and I wanna after Differentiating from it with respect to tau (at any amount of sigma value) and equaling this derivative to zero solve the infinite system of equations.

P[n](tau)==LegendreP(n - 1/2, cosh(tau)) , Q[n](tau)==LegendreQ(n - 1/2, cosh(tau)) 

are Legendre function.

Thanks in advanced

FUNCTION_f.mw

Download FUNCTION_f.mw

P1 := x^2+y^2-4:
P2 := y^2-2*x+2:

Original question is find CAD of (some y)[P1 <0 and P2 <0]

how to use maple 12 and maple 2015 to find Q1,Q2,Q3 which are projection of P1 and P2

my book show sample points are [-4,-1-sqrt(7),-3,-2,0,1,3/2,-1+sqrt(7),9/5,2,3]
but FindSamples result is not the same with my book, is it my book wrong or FindSamples function wrong?
I find result of my script is the same as book's quantifier position at 7,8,9 though sample points has little different

how to generalize my following script to multiple variables x, y, z, and more ?

and

I compare with maple 2015 result are different from my book solution, is maple 2015 more advanced version CAD? 

with(ListTools):

P1 := x^2+y^2-4:
P2 := y^2-2*x+2:

Q1 := x^2 + 2*x - 6;
Q2 := x^2 - 4;
Q3 := x - 1;

sourcesamples := sort(evalf([solve(Q1), solve(Q2), solve(Q3)]),`<`);

FindSamples:=proc(sourcesamples)
local N, P;
N:=nops(sourcesamples);
P:=proc(a,b)
local a1, b1, m1, n, m;
if a=b then error "Should be a<>b" fi;
a1,b1:=op(convert(sort([a,b],(x,y)->evalf(x)<evalf(y)),rational));
count := 0:
for n from 1 do
m1:=a1*n;
m:=`if`(type(m1,integer),m1+1,ceil(m1));
count := count + 1:
if is(m/n>a1) and is(m/n<b1) then return m/n fi;
od;
print("count=",count);
end proc:
[ceil(sourcesamples[1])-1, seq(op([sourcesamples[i],P(sourcesamples[i],sourcesamples[i+1])]), i=1..N-1),sourcesamples[N],floor(sourcesamples[N])+1];
end proc:

RemoveComplex := proc(yy)
local result, k:
result := []:
for k in yy do
if Im(k) = 0 then
result := [op(result), k]:
end if:
od:
if result = [] then
result := []:
end if:
return result:
end proc:

Joinsolution := proc(param1, param2group)
local result, k:
result := []:
for k in param2group do
result := [op(result), [param1, k]]:
od:
return result:
end proc:

CADsamples := FindSamples(sourcesamples):
CADresult1 := []:
for mm in CADsamples do
#print(mm):
if MakeUnique(RemoveComplex([solve(subs(x=mm, P1)), solve(subs(x=mm, P2))])) = [] then
CADresult1 := [op(CADresult1), op(Joinsolution(mm,[0]))];
else
CADresult1 := [op(CADresult1), op(Joinsolution(mm,FindSamples(sort(evalf(MakeUnique(RemoveComplex([solve(subs(x=mm, P1)), solve(subs(x=mm, P2))]))),`<`))))];
end if:
od:
CADresult1;

for mm in CADresult1 do
if subs(x=mm[1],subs(y=mm[2], P1)) < 0 and subs(x=mm[1],subs(y=mm[2], P2)) < 0 then
print("solution ", mm, SearchAll(mm[1],CADsamples), evalf(mm)):
end if:
od:

Compare with

with(RegularChains):
with(ChainTools):
with(MatrixTools):
with(ConstructibleSetTools):
with(ParametricSystemTools):
with(SemiAlgebraicSetTools):
with(FastArithmeticTools):
R := PolynomialRing([x,y]):
sys := [x^2+y^2-4,y^2-2*x+2]:
N := []:
P := []: 
H := [x]:
dec := RealTriangularize(sys,N,P,H,R):
proj := Projection(sys, N, P, H, 1, R);
Display(dec, R);

P := SamplePoints(sys, R);
Display(P, R);
cad := CylindricalAlgebraicDecompose(sys, R);
 

Hello everyone, my computer installed Maple 2016, 2020 and set Maple 2016 as default opening files .mw, .mws. How to change default to Maple 2020 . I tried searching on Google but not found.

Thank you very much.

How I can prove the following equation in red box.

Also, Pn(v) and qn(v) are the real combinations of half-integer Legendre functions.

For more details please see 

https://math.stackexchange.com/questions/2746660/potential-flow-around-a-torus-laplace-equation-in-toroidal-coordinates/3809487#3809487

Determine how many transitive relations there are on a set with n elements for all positive integers n with n<=7.

Custom Palettes and Palette Entries on windows Pro 7 
Read a tutorial from a user @les 185

https://www.mapleprimes.com/posts/207781-Custom-Palettes-And-Palette-Entries

Is this still valid this tutorial from 2015  in the current version of Maple ?
This seems to be userunfriendly.

Hi,

Thank you all for participating to my questions before.

I was wondering if we can change an equation into quartic form or quadratic form in maple.

I have succesfully got an equation like the following:

`S__2 ` = sin(alpha - phi)*sin(-beta + alpha)*(gamma*H^2*sin(beta - varepsilon)*sin(alpha) - h^2*sin(beta)*sin(alpha - varepsilon)*gamma + h^2*sin(beta)*sin(alpha - varepsilon)*psi)/(2*sin(-beta - delta - phi + alpha)*sin(beta)^2*sin(alpha - varepsilon)*sin(alpha)) - S__1

From the paper I studied, they change it to form an equation like this

`S__2 ` = 1/2*gamma*H^2*sin(beta - alpha)*(M__3(X)^3 + M__2(X)^2 + M__1*X)/(sin(beta)^2*(D__3(X)^3 + D__2(X)^2 + D__1*X + D__0)) - S__1

Where value of M, D and X is sets of long equation. Can someone teach me how to assign maple to change this kind of equation to another form of equation. It's good enough if I can learn how to learn the basic.

Thank you.

Regards

Faiz Farhan

Dear All,

I want to apply the ‘simplify’ command, in parallel, for the simplification of some parameters. Both Grid:-Map and Grid:-Run commands are tested. There is no error in both, whereas no simplification is implemented. It seems that the ‘simplify’ command correctly works on only ‘Master’ node, namely anywhere we are typing.

Can anyone help me to simplify in parallel. I examined two following codes.

1)

with (Grid);

for i from 1 to nops(dummy_UU1) do

freenode:=WaitForFirst():

Run(freenode,simplify,[dummy_UU1[i]],assignto='dummy_UU2'[i]):

end do:

Wait();

2)

dummy_UU2:=Map[tasksize=1](simplify,[seq(dummy_UU1[i],i=1..nops(dummy_UU1))]):

The following code is correctly executed and resulted in the simplification of dummy_UU1 components in serial.

for i from 1 to nops(dummy_UU1) do

dummy_UU2[i]:=simplify(dummy_UU1[i]):

end do:

Could anyone help me out to convert the equation into differential transform method

Hello Everyone, can anyone explain how to import a mathematical equation from maple to word directly?

I want to to solve the system of partial differential equation using maple. I tried it but I am not able to solve it ... please help.

the equations are as follows

Download Finding_transformation_eqn_between_zero_and_harmonic_with_conformal1.mw

Hi,

I am trying to solve a differentiation but I think I am stuck since the solution is not what it should be.

So, I got the equation below

eq4 := `S__2 ` = sin(alpha - phi)*sin(-beta + alpha)*H^2*M/(2*sin(-beta - delta - phi + alpha)*sin(beta)*sin(alpha)) - S__1

And according to the paper I read, to get the maximum value of alpha for maximum value of S_2, I need to make differentiation to first derivative where dS_2/d(alpha) = 0

Then I should substitute back value of alpha to equation above and the paper shows that i should get equation below.

`S__2 ` = 1/2*M*K__a*H^2 - `S__2 `

where K_a is

`K__a`= [(sin(beta+phi))/((sin(beta))/(sqrt(sin(beta+delta))+(sqrt(sin(phi+delta)*(sin(phi-varepsilon)))/(sin(beta-varepsilon))))]

I know its really hard but hope someone can give some idea how to do it.

Thank you very much.

Kind regards

Faiz Farhan

Got a lot of worksheets who are not complete anymore once opened in maple 2020

It can be only opened with a old version of Maple
Can it be imported in Maple 2020?

example 

Dynmod03.mws

Please can you help me in resolving this error?

Here is the codeOptimal_control_model_of_DF_and_LP_2.mw