L10

2.10.01

Käytännöllinen resonanssi

>

> g:=omega->m^2*(w0^2-omega^2)^2+c^2*omega^2;

g := proc (omega) options operator, arrow; m^2*(w0^...

> dg:=diff(g(omega),omega);

dg := -4*m^2*(w0^2-omega^2)*omega+2*c^2*omega

> wmax:=solve(dg=0,omega);

wmax := 0, 1/2*sqrt(4*m^2*w0^2-2*c^2)/m, -1/2*sqrt(...

> wmax:=%[2];

wmax := 1/2*sqrt(4*m^2*w0^2-2*c^2)/m

> simplify(g(wmax));

1/4*c^2*(-c^2+4*m^2*w0^2)/(m^2)

> Cti:=F0/simplify(sqrt(%),symbolic); # Cti tarkoittaa C* (KRE-merkinnöin)
# simplify(..,symbolic); antaa vastuun käyttäjälle.

Cti := 2*F0*m/(c*sqrt(-c^2+4*m^2*w0^2))

>