Dear Octave Community Together with the package maintainer, I'm working on a function called margin.m for the Octave Control Package, version 1.0.14:
http://octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/main/control/inst/margin.m?revision=6163&view=markup At the end of the file, there's a test. Unfortunately, we need very high tolerances (about 3000*eps) to get a pass, even on my own machine. Does somebody know what's going on here? I calculated the expected values with the following m-file: clear all; close all; clc; format long; num = [24]; den = [1, 6, 11, 6]; sysc = tf (num, den); sysd = c2d (sysc, 0.3); [gamma_c, phi_c, w_gamma_c, w_phi_c] = margin (sysc) [gamma_d, phi_d, w_gamma_d, w_phi_d] = margin (sysd) BTW: I'm using Octave.app version 3.2.2 on Mac OS X 10.5.8 Intel. Thanks for any help Lukas ------------------------------------------------------------------------------ Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day trial. Simplify your report design, integration and deployment - and focus on what you do best, core application coding. Discover what's new with Crystal Reports now. http://p.sf.net/sfu/bobj-july _______________________________________________ Octave-dev mailing list Octave-dev@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/octave-dev