clear all;
part_init;

u0 = 0; 
v0 = 0;
w0 = Wg;
x0 = Ro;
y0 = 0;
z0 = 0;

v = [u0; v0; w0; x0; y0; z0];

[f, pf_pw] = part_pf(v, p);

[VV,DD] = eig(pf_pw);

lambda = sort(abs(diag(DD)/(Wg/L)));

la_min = lambda(2);
la_max = lambda(6);

econd = la_max/la_min;

tcond = la_max;

fprintf('Eigenvalue condition number: %f\n',econd);
fprintf('T/tau_min  condition number: %f\n',tcond);



