Advanced Search

zy

30 Reputation

4 Badges

11 years, 210 days
Contact zy

MaplePrimes Activity


These are replies submitted by zy

the attachment...

zy 30
January 13 2014
0 0



X[1] := expand(subs(r = m, binomial(d, r)*binomial(v-d, r)));

binomial(v-d, d)

 (1)

 

a[1] := expand(binomial(d-1, r-1)*binomial(v-d-1, r)+binomial(d-1, r)*binomial(v-d-1, r-1)):

a[2] := expand(binomial(d-1, r)*binomial(v-d-1, r)):

a[3] := expand(binomial(d-1, r-1)*binomial(v-d-1, r-1)):

 (2)

b[1] := expand(binomial(d-1, r)*binomial(v-d-2, r-2)+binomial(d-1, r-1)*binomial(v-d-2, r)):

b[2] := expand(binomial(d-1, r-1)*binomial(v-d-2, r-1)):

b[3] := expand(binomial(d-1, r)*binomial(v-d-2, r-1)):NULL

b[4] := expand(binomial(d-1, r-1)*binomial(v-d-2, r-2)):

 (3)

b[5] := expand(binomial(d-1, r)*binomial(v-d-2, r)):

e[1] := expand(binomial(d-2, r)*binomial(v-d-1, r-1)+binomial(d-2, r-2)*binomial(v-d-1, r)):
NULL

e[2] := expand(binomial(d-2, r-1)*binomial(v-d-1, r-1)):

 (4)

e[3] := expand(binomial(d-2, r-1)*binomial(v-d-1, r)):NULL

e[4] := expand(binomial(d-2, r)*binomial(v-d-1, r)):NULL

e[5] := expand(binomial(d-2, r-2)*binomial(v-d-1, r-1)):

f[1] := expand(binomial(d-2, r)*binomial(v-d-2, r-2)+binomial(d-2, r-2)*binomial(v-d-2, r)+2*binomial(d-2, r-1)*binomial(v-d-2, r-1)):

f[2] := expand(binomial(d-2, r-1)*binomial(v-d-2, r-2)+binomial(d-2, r-2)*binomial(v-d-2, r-1)):

f[3] := expand(binomial(d-2, r)*binomial(v-d-2, r-1)+binomial(d-2, r-1)*binomial(v-d-2, r)):
NULL

f[4] := expand(binomial(d-2, r-1)*binomial(v-d-2, r-1)):NULL

f[5] := expand(binomial(d-2, r)*binomial(v-d-2, r)):
NULL

Q := proc (j) options operator, arrow; (d*(v-d)-v*j)*(v-1)/(d*(v-d)) end proc;

proc (j) options operator, arrow; (d*(v-d)-v*j)*(v-1)/(d*(v-d)) end proc

 (5)

d[0] := W[1]+W[2];

W[1]+W[2]

 (6)

c[0] := simplify(expand(subs(r = m, W[1]*(a[1]*Q(r)^2+a[2]*Q(r+1)^2+a[3]*Q(r-1)^2)/X[1])+subs(r = d, W[2]*(a[1]*Q(r)^2+a[2]*Q(r+1)^2+a[3]*Q(r-1)^2)/X[2]))):

c[11] := simplify(expand(subs(r = m, W[1]*(b[1]*Q(r)^2+2*b[2]*Q(r)*Q(r-1)+2*b[3]*Q(r)*Q(r+1)+b[4]*Q(r-1)^2+b[5]*Q(r+1)^2)/X[1])+subs(r = d, W[2]*(b[1]*Q(r)^2+2*b[2]*Q(r)*Q(r-1)+2*b[3]*Q(r)*Q(r+1)+b[4]*Q(r-1)^2+b[5]*Q(r+1)^2)/X[2]))):

c[12] := simplify(expand(subs(r = m, W[1]*(e[1]*Q(r)^2+2*e[2]*Q(r)*Q(r-1)+2*e[3]*Q(r)*Q(r+1)+e[4]*Q(r+1)^2+e[5]*Q(r-1)^2)/X[1])+subs(r = d, W[2]*(e[1]*Q(r)^2+2*e[2]*Q(r)*Q(r-1)+2*e[3]*Q(r)*Q(r+1)+e[4]*Q(r+1)^2+e[5]*Q(r-1)^2)/X[2]))):

NULL

NULL

 (7)

c[2] := simplify(expand(subs(r = m, W[1]*(f[1]*Q(r)^2+2*f[2]*Q(r)*Q(r-1)+2*f[3]*Q(r)*Q(r+1)+2*f[4]*Q(r+1)*Q(r-1)+f[5]*Q(r+1)^2+f[6]*Q(r-1)^2)/X[1])+subs(r = d, W[2]*(f[1]*Q(r)^2+2*f[2]*Q(r)*Q(r-1)+2*f[3]*Q(r)*Q(r+1)+2*f[4]*Q(r+1)*Q(r-1)+f[5]*Q(r+1)^2+f[6]*Q(r-1)^2)/X[2]))):

NULL

c[15] := simplify(expand(subs(r = m, W[1]*(a[1]*Q(r)+a[2]*Q(r+1)+a[3]*Q(r-1))/X[1])+subs(r = d, W[2]*(a[1]*Q(r)+a[2]*Q(r+1)+a[3]*Q(r-1))/X[2]))):

NULL

NULL

 (8)

with(LinearAlgebra):
NULL

B := Matrix([[c[0]+(v-d-2)*c[11], (d-1)*c[2], c[11], c[15]], [(v-d-1)*c[2], c[0]+(d-2)*c[12], c[12], c[15]], [(v-d-1)*c[11], (d-1)*c[12], c[0], c[15]], [(v-d-1)*c[15], (d-1)*c[15], c[15], d[0]]]):

C := Matrix([[c[0]+(v-d-2)*c[11], c[2]], [(v-d-1)*c[2], c[12]]]):

F[1] := -c[2]*(c[0]-c[12]):

h[1] := simplify(expand(subs(r = m, ((d-r)*G[2]+(c[0]-c[12])*G[1]+(r-1)*F[1])*Q(r-1)+((v-d-r)*F[1]-(d-r)*G[2])*Q(r)))):

h[2] := simplify(expand(subs(r = m, ((d-1)*G[2]+(c[0]-c[12])*G[1]+(r-1)*F[2]+(r-1)*F[3])*Q(r-1)+((v-d-r)*F[2]+(d-r)*F[3])*Q(r)))):

g[1] := simplify(expand(subs(r = d, ((d-r)*G[2]+(c[0]-c[12])*G[1]+(r-1)*F[1])*Q(r-1)+((v-d-r)*F[1]-(d-r)*G[2])*Q(r)))):

g[2] := simplify(expand(subs(r = d, ((d-1)*G[2]+(c[0]-c[12])*G[1]+(r-1)*F[2]+(r-1)*F[3])*Q(r-1)+((v-d-r)*F[2]+(d-r)*F[3])*Q(r)))):

R[7] := factor(simplify(expand(subs(a = 0, r[1] = m, r[2] = d, (r[1]-1)*h[1]*g[1]/((c[0]-c[12])*((d-r[1])*G[2]+(c[0]-c[12])*G[1])*((d-1)*G[2]+(c[0]-c[12])*G[1]))))));

R[8] := factor(simplify(expand(subs(a = 0, r[1] = m, r[2] = d, h[2]*g[2]/(G[3]*((d-1)*G[2]+(c[0]-c[12])*G[1]))))))

Error, Did not recognize this structure's syntax

R[8] := factor(simplify(expand(subs(a = 0, r[1] = m, r[2] = d, h[2]*g[2]/(G[3]*((d-1)*G[2]+(c[0]-c[12])*G[1]))))))

  

R[9] := factor(simplify(expand(subs(a = 0, r[1] = m, r[2] = d, h[3]*g[3]/(G[3]*G[4])))))

NULL



Download error.mw

@Preben Alsholm 

thanks...

zy 30
October 30 2013
0 0

@Carl Love

The problem is solved and thanks for your help. 

Page 1 of 1