To summarize, I have a variable ε = order(1), which maple has assumed is a funtion ε(x,y,z) and so when I differentiate epsilon with respect to x (or y or z) I do not get 0. I get εx (or εy, εz). How do I ensure maple does not assume this?
More detail of my process:
I declare functions,
I have the function I want to transform,
Now I transform the variables to the new co-ordinate system. i.e. from (x,n,q) to (s,Y,z)
Good! - Everything correct so far.
Now I want to linearise so i introduce x=x0+ε*x1; and the same for (Y,z),
As you can see, epsilon has derivatives, which it should not.
How I can avoid this?
Thanks in advance - im well and truely stumped over this.
P.s. if the images do not show, the script can be found here: