function f=Au(t,u)
%
global A
%
size(A)
size(u)
f=A*u;