When using the following sequence of commands, the sixth line doesn't have
the value of omega substituted as can be seen in the output of line 7 which
still displays omega and not the value I want it to have in line 6. Can
anyone tell me what I'm doing wrong and how to fix it?
h(s):=tau[1]*s/(1+(tau[1]+tau[2]+r[1]*c[2])*s+tau[1]*tau[2]*s^2);
s:%i*omega;
abs(h(s))^2;
diff(%,omega);
solve(%,omega);
omega:1/(sqrt(tau[1]*tau[2]));
abs(h(s));