#include "libscl.h"
using namespace scl;  
using namespace std;

REAL f(realmat& b) {b(1,1)=5.0; return b(1,1);}

int main(int argc, char** argp, char** envp)
{
  realmat a(5,5,0.0);
  REAL r = f(a);
  cout << "a(1,1) = " << a(1,1) << '\n';
  cout << "r = " << r  << '\n';
  return 0;
}