Solution

> plot( [ f, f_avg ], w=0..4, 0..0.1, title=`Function and average value` );

>

> sol := [ solve( f=f_avg, w ) ];

> w_avg := sol[1];

> evalf( w_avg );

>