Thanks, but it work only for text and not for variables and command (mod math) in expressions or equations

In Classic worksheet ok it work 

Could you give K2(x) ?

Thanks for the two answers

selectremove(has,f,i) do the trick

Thanks for the two answers

selectremove(has,f,i) do the trick

optimization.mws

optimization.mws

This is the same in mathematica, matlab, fortran etc...

This is the same in mathematica, matlab, fortran etc...

Yes, it is a challenge for Maplesoft

Thanks for your answer

Thank you for the answers

I will chose one way

Thank you for the answers

I will chose one way

In the case of your s1 and s2, you have a system of equations in f1 and f2 (if f[0](x) is known)

In this case, you can not solve it with a for command

you should write : dsolve({s[1](x),s[2](x)},{f[1](x),f[2](x)})

and add boundary conditions as (D@@2)(f[1]) because you have a third derivative

In the case of your s1 and s2, you have a system of equations in f1 and f2 (if f[0](x) is known)

In this case, you can not solve it with a for command

you should write : dsolve({s[1](x),s[2](x)},{f[1](x),f[2](x)})

and add boundary conditions as (D@@2)(f[1]) because you have a third derivative

@Markiyan Hirnyk 

particle.zip

Upload don't accept .nb file. I ziped it

I executed your maple worksheet and mathematica notebook.

In maple take a lot of time but in mathematica few seconds

I am interested also to know how to modify your commands to get more and more less time with maple ?