% TD7 - Van der Pol oscillator
% Describing function
%
% solution
clear all
close all

% definition du système
alpha = 1;

G = tf([alpha],[1 -alpha 1]); 

Ag = [0 1;-1 alpha];
Bg = [0 ; alpha];
Cg = [1 0];
Dg = 0;
sysG = ss(Ag,Bg,Cg,Dg);

% à faire :
% tracer sur la même figure
%       - le diagramme nyquist du systeme G
%       - le lieu critique -1/N(A)

A = 0:0.001:10;
w1 = ...;
NA = ...;