 # 2010-10-7

 clear all
 clc
 close all

 W = 0:0001:1000;
 R = 1; C = 1;

 x=(R.^2).*(C.^2).*(W.^2)./(1+(R.^2).*(C.^2).*(W.^2));
 y=R.*C.*W./(1+(R.^2).*(C.^2).*(W.^2));

 Zcua = x.^2 + y.^2;

loglog(W,Zcua)


# control, octave-control - control functions for Octave from Octave-Forge

# transfer function
s = tf('s');
H = R*C*S/(1+R*C*S);
bode(H);

