# chex32.r                    J F Monahan August 2007, June 2010
 #
 # Example 3.2 Solve a general system of linear equations
 #
  A <- matrix( c(1,.5,0,1, 2,1,2,-1, -1,0,-.5,1.5, 0,1,1.5,0), 4, 4)
 #  write out matrix
  A
  b <- c( .5, 1, 1.5, 2)
 #  write out right hand side
  b
 # solve equations
  x <- solve(A,b)
 #  write out solution
  x 
 # clean up
  rm(list=ls())            
  q()