(* Import LinearAlgebra package *) <1];]; (* Obtain timing data for Ones matrices *) tms={50,100,150,200}; For[i=1,i<5,i++, b=Table[Floor[Random[]*50],{j,tms[[i]]}]; A=Table[Table[If[Random[Integer]==0,-1,1],{j,tms[[i]]}],{k,tms[[i]]}]; Print[First[Timing[LinearSolve[A,b]]]/.Second->1];]; (* Obtain timing data for Pascal matrices *) tms={10,15,20,25,30}; For[i=1,i<6,i++, b=Table[Floor[Random[]*50],{j,tms[[i]]}]; A=Table[Table[Binomial[j+k,k],{j,tms[[i]]}],{k,tms[[i]]}]; Print[First[Timing[LinearSolve[A,b]]]/.Second->1];];