Advanced Search

vv

12950 Reputation

20 Badges

8 years, 362 days
Contact vv

MaplePrimes Activity


These are answers submitted by vv

Just use ReducedRowEchelon...

vv 12950
July 19 2023
3 2

Just use ReducedRowEchelonForm and append the identity matrix to your square matrix.

 

> 

restart;
> 

with(LinearAlgebra):
> 

n:=3:
> 

A:=RandomMatrix(n);

Matrix(3, 3, {(1, 1) = 27, (1, 2) = 99, (1, 3) = 92, (2, 1) = 8, (2, 2) = 29, (2, 3) = -31, (3, 1) = 69, (3, 2) = 44, (3, 3) = 67})

 (1)
> 

LinearAlgebra:-ReducedRowEchelonForm(<A|IdentityMatrix(n)>);

Matrix(3, 6, {(1, 1) = 1, (1, 2) = 0, (1, 3) = 0, (1, 4) = -3307/327244, (1, 5) = 2585/327244, (1, 6) = 5737/327244, (2, 1) = 0, (2, 2) = 1, (2, 3) = 0, (2, 4) = 2675/327244, (2, 5) = 4539/327244, (2, 6) = -1573/327244, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1, (3, 4) = 1649/327244, (3, 5) = -5643/327244, (3, 6) = 9/327244})

 (2)
> 

B:=%[..,n+1..];

Matrix(3, 3, {(1, 1) = -3307/327244, (1, 2) = 2585/327244, (1, 3) = 5737/327244, (2, 1) = 2675/327244, (2, 2) = 4539/327244, (2, 3) = -1573/327244, (3, 1) = 1649/327244, (3, 2) = -5643/327244, (3, 3) = 9/327244})

 (3)
> 

A.B; #check
 

Matrix(3, 3, {(1, 1) = 1, (1, 2) = 0, (1, 3) = 0, (2, 1) = 0, (2, 2) = 1, (2, 3) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1})

 (4)
> 

 

 

You just need a pencil for this.c=5, a=1...

vv 12950
July 11 2023
2 2

You just need a pencil for this.

c=5, a=12/2=6, b^2 = a^2 - c^2 = 11 

So, x^2/36 + y^2/11 = 1.

  I see the QuantifierElimi...

vv 12950
July 01 2023
1 1

 

I see the QuantifierElimination package (new in Maple 2023) as experimental.
Here is a much simpler problem and the solution obtained by  QuantifierElimination.

 

> 

restart;
> 

with( QuantifierElimination );

[CylindricalAlgebraicDecompose, DeleteFormula, InsertFormula, PartialCylindricalAlgebraicDecompose, QuantifierEliminate, QuantifierTools]

 (1)
> 

f:=(x-1)/(x^2+x+1);

(x-1)/(x^2+x+1)

 (2)
> 

fmax:=evala(maximize(f));

-1+(2/3)*3^(1/2)

 (3)
> 

forall(x, f <= K);

forall(x, (x-1)/(x^2+x+1) <= K)

 (4)
> 

sol:=QuantifierEliminate( % ); ### ?

And(Or(And(-K <= 0, Or(K <> 0, And(2*K <= 1, Or(2*K-1 <> 0, And(-K <= 0, Or(K <> 0, And(K <= 0, Or(K <> 0, -K < 1)))))))), And(K = 0, K-1 = 0, K+1 = 0)), Or(K <> 0, K <= 1, -K^5-2*K^4-4*K^3-6*K^2-3*K < 0, And(K^3+2*K^2+K = 0, K^2+3 <> 0)), Or(K <> 0, -K <= -1, -K^5-2*K^4-4*K^3-6*K^2-3*K < 0, And(K^3+2*K^2+K = 0, K^2+3 <> 0)), Or(K <> 0, K-1 = 0, And(-K^5-2*K^4-4*K^3-6*K^2-3*K <= 0, Or(K^5+2*K^4+4*K^3+6*K^2+3*K <> 0, And(2*K^5+7*K^4+6*K^2-18*K <= -3, Or(2*K^5+7*K^4+6*K^2-18*K+3 <> 0, And(-K^3-3*K^2-5*K <= -1, Or(K^3+3*K^2+5*K-1 <> 0, And(2*K^3+5*K^2-8*K <= -1, Or(2*K^3+5*K^2-8*K+1 <> 0, -K < 0)))))))), And(K = 0, K-1 = 0, K+1 = 0)), Or(-K <= 0, -3*K^2-6*K <= -1, -3*K^3+K^2 < 0, K^3 <> 0), Or(K <= 0, -3*K^2-6*K <= -1, -3*K^3+K^2 < 0, K^3 <> 0), Or(K <= 0, 3*K^2+6*K-1 <> 0, 3*K^2-K <> 0), Or(K <> 0, -K <= -1, And(-K^5-2*K^4-4*K^3-6*K^2-3*K <= 0, Or(K^5+2*K^4+4*K^3+6*K^2+3*K <> 0, And(2*K^5+7*K^4+6*K^2-18*K <= -3, Or(2*K^5+7*K^4+6*K^2-18*K+3 <> 0, And(-K^3-3*K^2-5*K <= -1, Or(K^3+3*K^2+5*K-1 <> 0, And(2*K^3+5*K^2-8*K <= -1, Or(2*K^3+5*K^2-8*K+1 <> 0, -K < 0)))))))), And(K = 0, K-1 = 0, K+1 = 0)), Or(K <= 0, -3*K^2-6*K <= -1, 3*K^3-K^2 < 0, K^3 <> 0), Or(-K <= 0, -3*K^2-6*K <= -1, 3*K^3-K^2 < 0, K^3 <> 0), Or(-K <= 0, 3*K^2+6*K-1 <> 0, 3*K^2-K <> 0), Or(K <= 0, -3*K^2-6*K <= -1, And(Or(3*K^3+6*K^2-K <= 0, 3*K^5+6*K^4-K^3 <= 0), Or(-3*K^2+K < 0, And(3*K^3+6*K^2-K <= 0, -3*K^5-6*K^4+K^3 <= 0)), Or(9*K^5+15*K^4-9*K^3+K^2 < 0, 3*K^5+6*K^4-K^3 <> 0, And(Or(2*K^2+5*K <= 1, K^3+5*K^2+6*K <= 1), 2*K^2+5*K <= 1, -K^3-5*K^2-6*K <= -1, Or(-2*K^2-5*K < -1, K^3+5*K^2+6*K-1 <> 0, And(Or(-K <= 0, 4*K^3+8*K^2-K <= 0), Or(K < 0, And(-K <= 0, -4*K^3-8*K^2+K <= 0)), Or(K^2 < 0, 4*K^3+8*K^2-K <> 0, -K < 0)))))), And(K = 0, K-1 = 0, K+1 = 0)), Or(-K <= 0, -3*K^2-6*K <= -1, And(Or(-3*K^3-6*K^2+K <= 0, -3*K^5-6*K^4+K^3 <= 0), Or(3*K^2-K < 0, And(-3*K^3-6*K^2+K <= 0, 3*K^5+6*K^4-K^3 <= 0)), Or(9*K^5+15*K^4-9*K^3+K^2 < 0, 3*K^5+6*K^4-K^3 <> 0, And(Or(-2*K^2-5*K <= -1, -K^3-5*K^2-6*K <= -1), Or(-2*K^2-5*K < -1, K^3+5*K^2+6*K-1 <> 0, And(Or(K <= 0, -4*K^3-8*K^2+K <= 0), Or(-K < 0, And(K <= 0, 4*K^3+8*K^2-K <= 0)), Or(K^2 < 0, 4*K^3+8*K^2-K <> 0, -K < 0)))))), And(K = 0, K-1 = 0, K+1 = 0)), Or(-K <= 0, -3*K^2-6*K <= -1, And(Or(-3*K^3-6*K^2+K <= 0, -3*K^5-6*K^4+K^3 <= 0), Or(-3*K^2+K < 0, And(-3*K^3-6*K^2+K <= 0, 3*K^5+6*K^4-K^3 <= 0)), Or(-9*K^5-15*K^4+9*K^3-K^2 < 0, 3*K^5+6*K^4-K^3 <> 0, And(Or(-2*K^2-5*K <= -1, -K^3-5*K^2-6*K <= -1), -2*K^2-5*K <= -1, K^3+5*K^2+6*K <= 1, Or(2*K^2+5*K < 1, K^3+5*K^2+6*K-1 <> 0, And(Or(K <= 0, -4*K^3-8*K^2+K <= 0), Or(K < 0, And(K <= 0, 4*K^3+8*K^2-K <= 0)), Or(-K^2 < 0, 4*K^3+8*K^2-K <> 0, -K < 0)))))), And(K = 0, K-1 = 0, K+1 = 0)), Or(K = 0, -3*K^2-6*K < -1, And(Or(3*K^4+6*K^3-K^2 <= 0, 3*K^6+6*K^5-K^4 <= 0), Or(-3*K^3+K^2 < 0, And(3*K^4+6*K^3-K^2 <= 0, -3*K^6-6*K^5+K^4 <= 0)), Or(9*K^5+15*K^4-9*K^3+K^2 < 0, 3*K^5+6*K^4-K^3 <> 0, And(Or(2*K^3+5*K^2-K <= 0, K^4+5*K^3+6*K^2-K <= 0), Or(K < 0, And(2*K^3+5*K^2-K <= 0, -K^4-5*K^3-6*K^2+K <= 0)), Or(-2*K^2-5*K < -1, K^3+5*K^2+6*K-1 <> 0, And(Or(-K^2 <= 0, 4*K^4+8*K^3-K^2 <= 0), Or(K^2 < 0, And(-K^2 <= 0, -4*K^4-8*K^3+K^2 <= 0)), Or(K^2 < 0, 4*K^3+8*K^2-K <> 0, -K < 0)))))), And(K = 0, K-1 = 0, K+1 = 0)), Or(K <= 0, -3*K^2-6*K <= -1, And(Or(3*K^3+6*K^2-K <= 0, 3*K^5+6*K^4-K^3 <= 0), Or(3*K^2-K < 0, And(3*K^3+6*K^2-K <= 0, -3*K^5-6*K^4+K^3 <= 0)), Or(-9*K^5-15*K^4+9*K^3-K^2 < 0, 3*K^5+6*K^4-K^3 <> 0, And(Or(2*K^2+5*K <= 1, K^3+5*K^2+6*K <= 1), Or(2*K^2+5*K < 1, K^3+5*K^2+6*K-1 <> 0, And(Or(-K <= 0, 4*K^3+8*K^2-K <= 0), Or(-K < 0, And(-K <= 0, -4*K^3-8*K^2+K <= 0)), Or(-K^2 < 0, 4*K^3+8*K^2-K <> 0, -K < 0)))))), And(K = 0, K-1 = 0, K+1 = 0)))

 (5)
> 

should_be = (K >= fmax);   # !!

should_be = (-1+(2/3)*3^(1/2) <= K)

 (6)
> 

is(eval(sol, K=fmax));            # ckeck
> 

is(eval(sol, K=fmax-10^(-20)));   #

true

  

false

 (7)

 

N.B. Why do you use such an unusual UserProfile?  The output of   2*x + 30  is:

 

3 half-lines...

vv 12950
June 26 2023
1 1

The system is not very simple, but I agree that SolveTools:-SemiAlgebraic is very slow in general, and I had not the patience for a result.

Anyway, the result (for eq1) seems to be (and is confirmed by MMA):

a := RootOf(_Z^3 - 4*_Z^2 - 27, 5.054 .. 5.063):
sol := {z=0, y = x, a < x}, {y=0, z = x, a < x}, {x=0, y = z, a < z};

plots:-spacecurve([[t,t,0], [t,0,t], [0,t,t]], t=a..3*a, color=[red,blue,green], thickness=5);

Answer...

vv 12950
June 12 2023
3 0

Can the equation (x-a)^4 + (x-b)^4 = c have four different integer solutions,where  a, b, c  are integers?

 

The answer is NO.  Actually not even four real solutions!

Suppose the polynomial  f := (x-a)^4 + (x-b)^4 - c  has 4 real solutions. After a shift, we may suppose that one of the solutions is 0, so c = a^4 + b^4,  a, b nonzero.

 

> 

restart;
> 

f := (x-a)^4 + (x-b)^4 - a^4 - b^4;

(x-a)^4+(x-b)^4-a^4-b^4

 (1)
> 

factor(f);

-2*x*(a+b-x)*(2*a^2-2*a*b-a*x+2*b^2-b*x+x^2)

 (2)

 

> 

f1:=op(-1,%);

2*a^2-2*a*b-a*x+2*b^2-b*x+x^2

 (3)

f1  must have two real solutions

> 

d1:=discrim(f1,x); # should be >0. But

-7*a^2+10*a*b-7*b^2

 (4)
> 

d2:=discrim(%,a);  # should be >= 0 (otherwise d2<0 implies d1<0).  ==>

-96*b^2

 (5)
> 

b:=0;  # contradiction.

0

 (6)
> 

solve(f1,x);  # Actually:

(1/2+((1/2)*I)*7^(1/2))*a, (1/2-((1/2)*I)*7^(1/2))*a

 (7)
> 

# Contradiction. Q.E.D.

 

ALG...

vv 12950
June 11 2023
1 0

The procedure ALG works only for finite sets! Actually, the support for infinite sets in Maple is very limited, see ?SetOf.

[0,2] is a list, not an interval. [0,1) is a nonsense in Maple (for 1D input).

Note also that ALG computes the algebra generated by C, not the sigma-algebra; of course, X, C must be finite, so this is the same thing.

Pure maths...

vv 12950
May 31 2023
2 0

You cannot solve such theoretical problems using Maple. You cannot even formulate them within Maple, because a CAS does other things!

But your problem is very simple if you know basic measure theory:

If A is a Borel subset in R^n and  g : A --> R  and  h : R^n \ A  --> R are continuous
then the function f : R^n --> R,  f(x) = g(x), if x in A, and f(x) = h(x), if x in R^n \ A
is a Borel function.

In your example, A = {(0,0)}; in your previous version A = {0} x R.

utf-8...

vv 12950
May 30 2023
4 1

Use:

LEN:=proc(s::string) Python:-EvalFunction("len",s) end proc:

It works for utf-8 characters, used by Maple, see https://mapleprimes.com/posts/214347-Multibyte-Characters

For example

LEN("România"); #  7

Note that the solution proposed by mmcdara works only for characters coded with at most 2 bytes. 

workaround...

vv 12950
May 29 2023
2 3

Use instead:

LC := proc(expr)
local i;
  if type(expr, {numeric,name,string}) then return 1  # or maybe atomic
  else add(i, i = map(thisproc,[op(expr)]))  +  1
  end if
end proc;


 

OK...

vv 12950
May 22 2023
1 1

Everything is correct, even for a without assumptions (i.e. for any complex a).

answer...

vv 12950
May 19 2023
0 0 
a:=sqrt(3):
plot3d(y^2,  x=y/a .. 1-y/a, y=0...a/2, orientation=[-75,60]);

du*dx etc...

vv 12950
May 09 2023
0 0

I suspect that you actually want du*dx instead of dudx etc. If so, try:

ss:=convert(aa,string):
ss:=StringTools:-SubstituteAll(ss,"dudx","du*dx"):
ss:=StringTools:-SubstituteAll(ss,"dudy","du*dy"):
ss:=StringTools:-SubstituteAll(ss,"dvdx","dv*dx"):
ss:=StringTools:-SubstituteAll(ss,"dvdy","dv*dy"):
ddaa:=parse(ss):
collect(ddaa,[du,dv,dx,dy],distributed);

 

Telescoping sums...

vv 12950
May 07 2023
2 1

Why do you complicate things?

The identity  

holds for k in Z, x in C  (actually for k,x in C). 
Summing for k in 0..n-1  we obtain S1 and summing for -k in 1..n  we obtain S6.

Workaround...

vv 12950
May 06 2023
0 2
> 

restart;
> 

expr := [
exp(y*LambertW(ln(y))),
(ln(y)/LambertW(ln(y)))^y,
eval(x^(x^x), x = exp(LambertW(ln(y)))),
eval(x^(x^x), x = ln(y)/LambertW(ln(y)))
]:
> 

Matrix(4, (i,j)->'is'(expr[i]=expr[j])):
simplify(%) assuming y>1;

Matrix(4, 4, {(1, 1) = true, (1, 2) = FAIL, (1, 3) = FAIL, (1, 4) = FAIL, (2, 1) = FAIL, (2, 2) = true, (2, 3) = FAIL, (2, 4) = FAIL, (3, 1) = FAIL, (3, 2) = FAIL, (3, 3) = true, (3, 4) = true, (4, 1) = FAIL, (4, 2) = FAIL, (4, 3) = true, (4, 4) = true})

 (1)

 

Workaround

 

> 

Y:=solve(LambertW(ln(y))=t, y);

exp(t*exp(t))

 (2)
> 

mat:=simplify(subs(LambertW(ln(y))=t, y=Y, expr)) assuming t>0;

[exp(exp(t*exp(t))*t), exp(exp(t*exp(t))*t), exp(exp(t*exp(t))*t), exp(exp(t*exp(t))*t)]

 (3)
> 

Matrix(4, (i,j)->is(mat[i]=mat[j]) );

Matrix(4, 4, {(1, 1) = true, (1, 2) = true, (1, 3) = true, (1, 4) = true, (2, 1) = true, (2, 2) = true, (2, 3) = true, (2, 4) = true, (3, 1) = true, (3, 2) = true, (3, 3) = true, (3, 4) = true, (4, 1) = true, (4, 2) = true, (4, 3) = true, (4, 4) = true})

 (4)

 

direct...

vv 12950
May 05 2023
2 0

We don't need Maple for this.

ΔABC is similar to ΔEBC and ΔFDC. So, D is the midpoint of BC. DF=AE=EB=3, DE=AF=FC=4, DB=DC=5 (Pythagoras).

H being the centroid of ΔABC we have AD=BD=DC=5, and finally DH = AD/3 = 5/3.

 

5 6 7 8 9 10 11 Last Page 7 of 115