function [xdot,Av] = fL63(x,v,param);

%sigma = param.sigma;
%beta = param.beta;
%r = param.r;
sigma = 10;
rho = 28;
beta = 8/3;

xdot = [sigma*(x(2)-x(1));
    	x(1)*(rho - x(3)) - x(2);
    	x(1)*x(2) - beta*x(3)];
    
A = [-sigma, sigma, 0;
    rho-x(3), -1, -x(1);
    x(2), x(1), -beta];
Av = A*v;
  