vv

14257 Reputation

20 Badges

11 years, 12 days

MaplePrimes Activity


These are replies submitted by vv

@minhthien2016 

The real option in limit(f(x), x=a, real)  makes sense only when a=infinity and means that limit(f(x), x=infinity) and limit(f(x), x=-infinity) both exist and are equal.

For example,
limit(sqrt(1-x^2)/abs(x), x=infinity, real);  # the limit is complex!
        I

 

@Rouben Rostamian  

Nice problem.
The exact answers are:
Q1:

P:= n -> add( binomial(n-5*k,k)/6^(6*k)*(-1)^(k-1), k=1..floor(n/6))

Q2:
N = 32341

(P(32340)<1/2, P(32341)>1/2)

 

@Rouben Rostamian  

It's only a (~good) approximation for small n.

E.g. P(12) = 7/6^6 - 1/6^12.

Probably the answer for Q2 is almost correct.

 

@Joe Riel 

It would be interesting to see the Maple answer in the form
Y := n -> [1,-1,0][n mod 3 +1]

@JAMET 

I don't see any problem.

 

> 

FC:= proc(n)
uses plots,plottools;
local a,b,u;
{seq(seq(a/b, a = 0 .. b), b = 1 .. n)};
display(seq(circle([u,1/2/denom(u)^2],1/2/denom(u)^2), u=%)
        ,scaling = constrained, color = blue, size = [800, 800])
end proc:

> 

FC(6);

 
> 

 


Download FC.mw

@JAMET 

u runs in the set {0, 1/n, ..., a/b, ..., 1}  constructed just above it.

Please explain for this simpler expression.

f:=diff(v(x, y), x, y) * diff(u(x, y), y, y)^2;

What are you expecting for int(f, x=0..1, y=0..1); ?

 

@Carl Love 

Are you saying that the global variables are saved in a .m file without being explicitely present in the arguments of save?

I have not seen this behavior. Probably you mean something else.

@acer 

Thank you, good to know. Probably I should have guessed it. Actually, if B is not definite, some eigenvalues could be undefined (even if symmetric and commuting A,B implies real eigenvalues -- but their number could be < n).

@Carl Love 

You are right, I have assumed (without checking) that A,B commute (A.B=B.A).

Edit. Even when A,B commute and have hfloats, the eigenvalues still have imaginary 0s.
So, the algorithm in this case ignores the symmetry.

@Carl Love 

Yes, M is symmetric, but it seems that this is ignored by Eigenvalues, because 0 imaginary parts still appear. Probably the symmetry is taken into account only for hfloats.

@tomleslie

So, you think that these are OK?

 

> 

restart;

> 

J:= Matrix(2, shape=identity);

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

(1)
> 

K:= Matrix(2, shape=triangular[lower,unit]);

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

(2)
> 

2*J-K, 2*K-J, 3*J-2*K;
J-K, K-J;

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

 

Matrix(%id = 18446744074412918894), Matrix(%id = 18446744074412919494)

(3)
> 

 

Download JK.mw

@tomleslie 

The OP does not want to reassign anything; he wants some computations with U.

@mmcdara 

You missed 13 cycles.

{[1,3,5],[1,3,6],[1,4,5],[1,4,6],[2,3,5],[2,3,6],[2,4,5],[2,4,6],[1,3,2,4],[1,3,2,5],[1,3,2,6],[1,3,5,4],[1,3,6,4],[1,4,2,5],[1,4,2,6],[1,5,2,6],[1,5,3,6],[1,5,4,6],[2,3,5,4],[2,3,6,4],[2,5,3,6],[2,5,4,6],[3,5,4,6],[1,3,2,4,5],[1,3,2,4,6],[1,3,2,5,4],[1,3,2,6,4],[1,3,5,2,4],[1,3,5,2,6],[1,3,5,4,6],[1,3,6,2,4],[1,3,6,2,5],[1,3,6,4,5],[1,4,2,3,5],[1,4,2,3,6],[1,4,5,2,6],[1,4,5,3,6],[1,4,6,2,5],[1,4,6,3,5],[1,5,2,3,6],[1,5,2,4,6],[1,5,3,2,6],[1,5,4,2,6],[2,3,5,4,6],[2,3,6,4,5],[2,4,5,3,6],[2,4,6,3,5],[1,3,2,5,4,6],[1,3,2,6,4,5],[1,3,5,2,4,6],[1,3,5,2,6,4],[1,3,5,4,2,6],[1,3,6,2,4,5],[1,3,6,2,5,4],[1,3,6,4,2,5],[1,4,2,5,3,6],[1,4,2,6,3,5],[1,4,5,2,3,6],[1,4,5,3,2,6],[1,4,6,2,3,5],[1,4,6,3,2,5],[1,5,3,2,4,6],[1,5,4,2,3,6]}

@weidade37211 

When working with such "extreme" numbers a good idea is to transform the equation symbolically.
 

> 

c = p_l^k*(1-p_l)^(n-k)*theta/(p_l^k*(1-p_l)^(n-k)*theta+p^k*(1-p)^n*(1-theta));

c = p_l^k*(1-p_l)^(n-k)*theta/(p_l^k*(1-p_l)^(n-k)*theta+p^k*(1-p)^n*(1-theta))

(1)
> 

eq:=eval(%,[k=1, theta=9/10]);

c = (9/10)*p_l*(1-p_l)^(n-1)/((9/10)*p_l*(1-p_l)^(n-1)+(1/10)*p*(1-p)^n)

(2)
> 

isolate(eq, p*(1-p)^n);

p*(1-p)^n = 9*p_l*(1-p_l)^(n-1)/c-9*p_l*(1-p_l)^(n-1)

(3)
> 

factor(map(ln,%));

ln(p*(1-p)^n) = ln(-9*p_l*(1-p_l)^(n-1)*(c-1)/c)

(4)
> 

expand(%) assuming p<1,p_l<1,p_l>0,n>0;

n*ln(1-p)+ln(p) = 2*ln(3)+ln(p_l)+(n-1)*ln(1-p_l)+ln((-c+1)/c)

(5)
> 

EQ:=eval(%, [c=95/100, p_l=10^(-15), n=10^10, c=95/100]);

10000000000*ln(1-p)+ln(p) = 2*ln(3)-15*ln(10)+9999999999*ln(999999999999999/1000000000000000)-ln(19)

(6)
> 

fsolve(EQ);

0.1496592586e-8

(7)
> 

Digits:=30;

30

(8)
> 

fsolve(EQ);

0.149659258635919070563992798396e-8

(9)
First 69 70 71 72 73 74 75 Last Page 71 of 178