Advanced Search

vv

14122 Reputation

20 Badges

10 years, 191 days
Contact vv

MaplePrimes Activity


These are answers submitted by vv

Solution...

vv 14122
March 28 2026
1 0

    Note first that if the triangle SBC is not supposed to be isosceles and BC is not parallel to the plane π then the minimum is 0, S being the intersection between BC and the plane π (so SBC is of course degenerated).
    
    If the triangle is isosceles (SB=SC) and BC is not perpendicular to the plane π, then S must belong to the mediator plane π′ of the segment BC (i.e. the plane perpendicular on BC through the midpoint M of the segment BC).
    So, S∈π∩π′. For a minimal area, S must be the projection of the point M onto the straight line π∩π′.

(For a Maple solution one has to find the equation of π∩π′ and minimize the quadratic expression MS^2 for S∈π∩π′).

Modifications...

vv 14122
March 04 2026
1 0 
A1 := 
plot([L12, L22], Ce = 0 .. 0.8, 
color = ["Red", "Blue"], 
linestyle = [solid, longdash], 
labels = [C__e, pi__r], labeldirections = ["horizontal", "vertical"], 
legend = ['one'=0.11, 'two'=0.12], legendstyle=['location=right'], 
axis[2] = [color = "#600000"]);

Similar for A2, A3.

x_...

vv 14122
February 04 2026
1 0 
F:=x-> (x^4+3*x^3-3*x^2-2*x-24)/(x^4-4*x^3-13*x^2+62*x-56);
u:=x->piecewise(x=-4, limit(F(x_), x_=x), true, F(x));   # or x_=-4
u(-4);

convert/sin,cos...

vv 14122
January 23 2026
3 0

Local solution...

vv 14122
December 13 2025
3 2

The solution has not an unexpected behavior.

The rhs of the ode is obviously C^1, so locally Lipschitz, so the solution exists and is unique in a nbd of IC.
The solution found by dsolve is valid in  `#mrow(mi("a"),mo(","),mo("∞"))`

where a=-1.3...   (see below).

 

A sufficient condition for F(x,y)  to be Lipschitz wrt  y near (x0,y0)  is that D[2](F) be without singularities near (x0,y0).   

 

> 

F := (x,y) -> x^2 - y^2;
ode := diff(y(x), x) = x^2 - y(x)^2:
IC := y(0) = 1:
sol := dsolve([ode, IC]);

F := proc (x, y) options operator, arrow; x^2-y^2 end proc

  

y(x) = piecewise(x < 0, x*((1/2)*BesselI(-3/4, (1/2)*x^2)*Pi*(GAMMA(3/4)^2*2^(1/2)-Pi)/GAMMA(3/4)^2-BesselK(3/4, (1/2)*x^2))/((1/2)*Pi*(GAMMA(3/4)^2*2^(1/2)-Pi)*BesselI(1/4, (1/2)*x^2)/GAMMA(3/4)^2+BesselK(1/4, (1/2)*x^2)), x = 0, 1, 0 < x, x*((1/2)*BesselI(-3/4, (1/2)*x^2)*Pi*(GAMMA(3/4)^2*2^(1/2)+Pi)/GAMMA(3/4)^2-BesselK(3/4, (1/2)*x^2))/((1/2)*Pi*(GAMMA(3/4)^2*2^(1/2)+Pi)*BesselI(1/4, (1/2)*x^2)/GAMMA(3/4)^2+BesselK(1/4, (1/2)*x^2)))

 (1)
> 

Y1:=rhs(sol) assuming x<0: Y2:=rhs(sol) assuming x>0:
> 

plot([Y1, Y2], x=-3..3);

 
> 

plot(rhs(sol), x=-3..3, discont=[usefdiscont], title=`ode solution` );

 
> 

fdiscont(Y1, x=-2..-1);

[HFloat(-1.0378337700000781) .. HFloat(-1.0368198762094487)]

 (2)
> 

discont(D[2](F)(x,y), x), discont(D[2](F)(x,y), y);

{}, {}

 (3)
> 

limit(D[2](F)(x,y), [x=0,y=1]);

-2

 (4)


Download LipODE-vv.mw

Even in the plane...

vv 14122
November 23 2025
4 0

Such triangles exist even in the plane. (AE is the angle bisector in  ΔABC.)

restart;

IntegerTriangle := proc(iter)
local p,q,r,s, A:=[0,0], B,C,E,AB,AC,BC,AE,  den,
      nor:=V -> sqrt(V[1]^2+V[2]^2):
for p to iter do
for q to p-1  do   if igcd(p,q)>1 or irem(p,2)=irem(q,2) then next fi;
for r from p to iter do
for s to r-1  do   if igcd(r,s)>1 or irem(p,2)=irem(q,2) then next fi;
  B:=[p^2-q^2, 2*p*q]: AB:=p^2+q^2:
  C:=[r^2-s^2, 2*r*s]: AC:=r^2+s^2:
  BC:=nor(B-C):
  if not type(BC, rational) then next fi:
  E:=(AC*~B+AB*~C)/~(AC+AB):
  AE:=nor(E):
  if not type(AE, rational) then next fi:
  if (AB=AC)or(AB=BC)or(AC=BC) then next fi:
  if (AB+AC+BC)<=2*max(AB,AC,BC) then next fi:
  if (AB^2=AC^2+BC^2)or(AC^2=AB^2+BC^2)or(BC^2=AB^2+AC^2) then next fi:
  den:=ilcm(denom~([E[],AE])[]): B:=B*den: C:=C*den: E:=E*den;
  print(ABCE=[A,B,C,E]);
od od od od:
print("OK")
end proc:
IntegerTriangle(100);

 ABCE = [[0, 0], [4797, 1804], [4797, 280604], [4797, 6804]]
 ABCE = [[0, 0], [697, 14280], [18377, 14280], [7425, 14280]]
ABCE = [[0, 0], [1989, 116348], [55029, 45628], [34845, 72540]]
ABCE = [[0, 0], [7605, 444860], [226395, 153140], [143136, 264152] ]
"OK"

CompleteSquare...

vv 14122
November 17 2025
3 4

 

restart;
eq := x^2 + y^2 + z^2 - 2*x + 8*y - 6*z - 30 = 0:
Student[Precalculus][CompleteSquare](eq);

                            (z-3)^2 + (y+4)^2 + (x-1)^2 - 56  =  0

Unevaluate 'expr'...

vv 14122
November 15 2025
2 0

You must unevaluate the argument of latex using right single quotes 'expr'.

latex('diff(phi(t), t) = (`&Dscr;`[1, 4](phi(t)) + `&Dscr;`[4, 1](phi(t)))*diff(phi(t), t, x) + (`&Dscr;`[2, 4](phi(t)) + `&Dscr;`[4, 2](phi(t)))*diff(phi(t), t, y) + (`&Dscr;`[3, 4](phi(t)) + `&Dscr;`[4, 3](phi(t)))*diff(phi(t), t, z) + (`&Dscr;`[1, 2](phi(t)) + `&Dscr;`[2, 1](phi(t)))*diff(phi(t), x, y) + (`&Dscr;`[1, 3](phi(t)) + `&Dscr;`[3, 1](phi(t)))*diff(phi(t), x, z) + (`&Dscr;`[2, 3](phi(t)) + `&Dscr;`[3, 2](phi(t)))*diff(phi(t), y, z) + `&Dscr;`[4, 4](phi(t))*diff(phi(t), t, t) + `&Dscr;`[1, 1](phi(t))*diff(phi(t), x, x) + `&Dscr;`[2, 2](phi(t))*diff(phi(t), y, y) + `&Dscr;`[3, 3](phi(t))*diff(phi(t), z, z) + (diff(phi(t), t) + diff(phi(t), x) + diff(phi(t), y) + diff(phi(t), z))*(D(`&Dscr;`[1, j])(phi(t))*diff(phi(t), x) + D(`&Dscr;`[2, j])(phi(t))*diff(phi(t), y) + D(`&Dscr;`[3, j])(phi(t))*diff(phi(t), z) + D(`&Dscr;`[4, j])(phi(t))*diff(phi(t), t))');

Otherwise your expression is evaluated, and for example diff(phi(t), t, x) evaluates to 0.

A way...

vv 14122
November 02 2025
1 0

A way to obtain the plot file is the following:

In the worksheet "f.mw" subject to RunWorksheet, don't use plotsetup but use:

...
p := plot(...):
...
return p;

(so, return the PLOT structure).
In the calling worksheet use:

ret := DocumentTools :-RunWorksheet("f.mw"):
plottools:-exportplot( "d:/temp/myplot.png", ret);

Obvious general solution...

vv 14122
November 02 2025
1 1

General solution: p, q, r, n arbitrary integers (n>0) and s = k^3 - p*A - q*B - r*C, where k is an arbitrary integer.
(Not very interesting).

Primfield...

vv 14122
November 01 2025
4 1
> 

restart;
> 

u := convert(3^(1/3),RootOf);
v := convert(4^(1/4), RootOf);
uv1 := u + v + 1^(1/1);

RootOf(_Z^3-3, index = 1)

  

RootOf(_Z^2-2, index = 1)

  

RootOf(_Z^3-3, index = 1)+RootOf(_Z^2-2, index = 1)+1

 (1)

"I would like to represent  u  and  v  as elements of the field generated by a = u+v+1" .

> 

pf:=evala(Primfield({u,v}));

[[RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1) = RootOf(_Z^2-2, index = 1)+RootOf(_Z^3-3, index = 1)], [RootOf(_Z^2-2, index = 1) = (48/755)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)^5+(27/755)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)^4-(64/151)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)^3-(468/755)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)^2+(879/755)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)-1092/755, RootOf(_Z^3-3, index = 1) = -(124/755)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)-(48/755)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)^5-(27/755)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)^4+(64/151)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)^3+(468/755)*RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)^2+1092/755]]

 (2)
> 

op([1,1,1],pf);

RootOf(_Z^6-6*_Z^4-6*_Z^3+12*_Z^2-36*_Z+1)

 (3)
> 

alias(alpha=op([1,1,1],pf));

alpha

 (4)
> 

evala(eval(uv1, pf[2]));

alpha+1

 (5)
> 

simplify(eval(pf[2], alpha=a-1));

[2^(1/2) = (48/755)*a^5-(213/755)*a^4+(52/755)*a^3+(174/755)*a^2+(987/755)*a-428/151, 3^(1/3) = -(232/755)*a+277/151-(48/755)*a^5+(213/755)*a^4-(52/755)*a^3-(174/755)*a^2]

 (6)
> 

simplify(eval(%, a=u+v+1)); # check

[2^(1/2) = 2^(1/2), 3^(1/3) = 3^(1/3)]

 (7)


Download Primfield-vv.mw

sort, key=evalf...

vv 14122
October 30 2025
4 1

Just replace
S := sort([sqrt(x2), sqrt(y2), sqrt(z2)]);
with

S := sort([sqrt(x2), sqrt(y2), sqrt(z2)], key=evalf);

This is necessary because we need a numeric rather than symbolic ordering (see ?sort).

Solution...

vv 14122
October 26 2025
2 4 
Q5:=n -> irem(add(convert(n,base,10)),5):
for n to 100000 do
if (Q5(n)+Q5(n+1)=0) then print(nmin=n); break fi od:

                            nmin = 49999

answer...

vv 14122
October 25 2025
2 0
> 

 

During a birthday party, the birthday child realizes:

In 1968, I was the same age as the sum of the digits of my birth year. How old will I be now at the end of 2025?

> 

restart;
> 

# y = birth year, 1800<y<1968
> 

s:=y -> add(convert(y, base,10))

proc (y) options operator, arrow; add(convert(y, base, 10)) end proc

 (1)
> 

eq := 1968-y = 's'(y);

1968-y = s(y)

 (2)
> 

for y from 1800 to 1968 do
if eq then Y:=y; fi od;
> 

Y, 2025-Y;

1947, 78

 (3)
 

 

Download birth-vv.mw

Solution...

vv 14122
October 25 2025
5 0
> 

restart

Denote by u the number with first digit (7) deleted and by n the number of its digits.

Then:

> 

(10*u+7)*3 = 7*10^n+u;

30*u+21 = 7*10^n+u

 (1)
> 

u:=solve(%,u);

(7/29)*10^n-21/29

 (2)
> 

for n do
if type(u,integer) then print(n); break fi od:

27

 (3)
> 

k=u+7*10^n;

k = 7241379310344827586206896551

 (4)
 

 

Download sol-vv.mw

1 2 3 4 5 6 7 Last Page 1 of 121