--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/lisovfit/test_lisovfit.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,91 @@
+mc
+
+%%
+% make ao of desired response
+apl = plist('filename', 'mdc2_1_tf11_np.dat', 'type', 'fsdata')
+a = ao(apl);
+
+% now make LISO fit file
+
+% start with a guess pzmodel
+
+z0(1) = pz(1e-6);
+zl(1) = pz(1e-07);
+zu(1) = pz(1e-05);
+
+p0(1) = pz(400e-9);
+pl(1) = pz(1e-07);
+pu(1) = pz(1e-05);
+
+p0(2) = pz(1e-6);
+pl(2) = pz(1e-07);
+pu(2) = pz(1e-05);
+
+p0(3) = pz(30e-6);
+pl(3) = pz(1e-06);
+pu(3) = pz(1e-4);
+
+z0(2) = pz(50e-6);
+zl(2) = pz(1e-06);
+zu(2) = pz(1e-4);
+
+z0(3) = pz(50e-6);
+zl(3) = pz(1e-06);
+zu(3) = pz(1e-4);
+
+z0(4) = pz(50e-6);
+zl(4) = pz(1e-06);
+zu(4) = pz(1e-4);
+
+z0(5) = pz(50e-6);
+zl(5) = pz(1e-06);
+zu(5) = pz(1e-4);
+
+p0(4) = pz(7e-3);
+pl(4) = pz(1e-4);
+pu(4) = pz(1e-1);
+
+p0(5) = pz(7e-3);
+pl(5) = pz(1e-4);
+pu(5) = pz(1e-1);
+
+p0(6) = pz(40e-3);
+pl(6) = pz(1e-4);
+pu(6) = pz(1);
+
+p0(7) = pz(40e-3);
+pl(7) = pz(1e-4);
+pu(7) = pz(1);
+
+z0(6) = pz(0.5);
+zl(6) = pz(0.01);
+zu(6) = pz(1);
+
+p0(8) = pz(0.1);
+pl(8) = pz(0.01);
+pu(8) = pz(10);
+
+z0(6) = pz(0.5);
+zl(6) = pz(0.01);
+zu(6) = pz(1);
+
+
+pzm0 = pzmodel(1e-10, p0, z0);
+pzml = pzmodel(1e-12, pl, zl);
+pzmu = pzmodel(1e-8, pu, zu);
+
+resp(pzm0,pzml,pzmu)
+
+
+%%
+  fpl = plist('PZM0', pzm0, 'PZML', pzml, 'PZMU', pzmu, 'delay', [0 1 10], 'f1', 1e-4, 'f2', 1, 'nf', 100);
+
+pzmfit = lisovfit(a, fpl);
+
+%%
+
+f0   = resp(pzm0);
+rfit = resp(pzmfit);
+
+iplot(f0, rfit)
+