Greetings,

I have a differential equation:

an the solution:

I want to substitute

ekf := omega = sqrt(c/m);

algsubs(ekf, sol);
but sol does not change.

This works a bit but does not substitute everything in a way that you would expedect

This is the result:

Thanks in advance