Advanced Search
All Posts : MaplePrimes Commons General Technical Discussions
hi , can anyone tell me what's wrong with my codes ?? THANK YOU > restart: > with(LinearAlgebra): > A:=Matrix(2,2,[[100,600],[200,400]]); [100 600] A := [ ] [200 400] > list2:='list2': > list1:=list1: > n1[0]:=200: n2[0]:=200: list1[1]:=[0,n1[0]]: > for k from 1 to 1000 do > alpha[k]:=n1[k-1]*(100*n1[k-1]+600*n2[k-1]); > beta[k]:=n1[k-1]*(n1[k-1]+n2[k-1]); > gammas[k]:=(n2[k-1]*(200*n1[k-1]+400*n2[k-1]))/(n1[k-1]+n2[k-1]); > five[k]:=n2[k-1]*(n1[k-1]+n2[k-1]); > alphahat[k]:=(alpha[k])/(alpha[k]+beta[k]+gammas[k]+five[k]); > betahat[k]:=beta[k]/(alpha[k]+beta[k]+gammas[k]+five[k]); > gammahat[k]:=gammas[k]/(alpha[k]+beta[k]+gammas[k]+five[k]); > fivehat[k]:=five[k]/(alpha[k]+beta[k]+gammas[k]+five[k]); > r:=rand()*10.0^(-12);evalf(r); > if r

Please Wait...