2304 Reputation

9 Badges

1 years, 240 days

MaplePrimes Activity

These are replies submitted by vv


The ode could be in an interval (0,a), a>0, and f'' unbounded.


You are right, I missed the "e-2" at the end.
For me it is more annoying that LSSolve is not able to manage Digits>15.


@Carl Love 

Yes, and it works also for X(n) for n>4  (see a previous comment).
But it does not work for arctan because of combine. E.g.

x1 := arctan(-sqrt(2)-1+sqrt(2-sqrt(2))+sqrt(2+sqrt(2)))

(The tan(k*x1) approach works).



Of course. I wonder if such a systematic approach will work for


or n>4. (X(4) = x1).


@Rouben Rostamian  

rij hat  is a versor (= rij/rij) , so the bold formula is OK

@Joe Riel 

Unfortunately Iterator will be anyway limited by the huge dimension of the discrete structure.
It will be not possible to use Permute for n=15 because 15! > 10^12;  20! > 10^18.
We will have to wait for the quantum computers.

@Joe Riel 

Thank you for the answer.
It is also strange for me that the procedure is also very slow for the very simple graph defined by

A:=Matrix(n, (i,j)->`if`(i=j,0,i+j)):
and here the triangle inequality holds)

but becomes very fast for e.g.

A:=Matrix(n, (i,j) -> evalf(abs(sin(i)-sin(j))));
(at least for n=10).


Very nice, vote up.
As optimizations I'd suggest:

1. Remove  HC(S):= P[k]  and recompute it at the end of MinPart (being a waste of memory).
2. Probably  combinat:-permute  should be called separately (at start) only r times for MinPart([m1,...,mr])
and then only translate the indices.

@Carl Love 


@Carl Love 

Yes, I know and I wanted to mention this. The problem is that for each set of each partition we must find  a minimal Hamiltonian path so, probably the running time will be prohibitive (even if parallel computations are used).

(But many sets will be eliminated because of the capacity constraint, so after all the attempt could be succesful).

@Carl Love 

An execution group may have several prompts (why?). For example, if a .mpl file is opened (via File > Open) all its contents is placed into a single execution group and each line has a prompt. The only way I know to delete the prompts is to copy & paste into Notepad (to remove the format information) and then paste back. Do you know another method? Thank you.

@Markiyan Hirnyk 

Simply solving the first 3 inequalities wrt y3 and using
y3 >= a, y3 >= b, y3 >= c  <==>  y3 >= max(a,b,c).

@Markiyan Hirnyk 

But compare with a "by hand" solution:

sol:={ y1 <= 0, y2 >= 0, y3 >= max( -(4/3)*y1-(2/3)*y2-5/3, -(3/2)*y1-(5/2)*y2-1, -y1-2*y2-1 ) };


@Carl Love 

Yes, but not "much more willing"; e.g.  for qgamma(0.8, 0.999);  # should be 1.164159836


@Joe Riel 

R := rtable(antisymmetric,1..3, 1..3, (i,j)->i-j);




I agree, but I think that the command should have been programmed more carefully.

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