--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/src/ltpda_ssmsim/do_a_run.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,21 @@
+function [t,y,x] = do_a_run(Nruns, Nsamples, Nstates, Nstatesout, Ninputs, Noutputs)
+  
+  x = [];
+  
+  lastX     = zeros(Nstates,1);
+  A         = 0.1*rand(Nstates);
+  Coutputs  = rand(Noutputs,Nstates);
+  Cstates   = eye(Nstates);
+  Cstates   = Cstates(1:Nstatesout,1:Nstates);
+  Baos      = rand(Nstates,Ninputs);
+  Daos      = rand(Noutputs, Ninputs);
+  input     = randn(Ninputs, Nsamples);
+  
+  
+  tic
+  for rr=1:Nruns
+    [y,lx] = ltpda_ssmsim(lastX, A.', Coutputs.', Cstates.', Baos.', Daos.', input);
+  end
+  mxtime = toc;
+  t = mxtime/Nruns;
+end