You are now following this question
- You will see updates in your followed content feed.
- You may receive emails, depending on your communication preferences.
how can i substitute some variables in a function with except of others
1 view (last 30 days)
Show older comments
I have an equation with 137 parts, each part has its variables. I need to replace some variables in these parts and except others.
In the following is the example of equation that i want to process.
y = a0*b0*c0*k0 + a1*b0*c0*d*k0 + a0*b0*c0*k4*cos((pi*m)/3) + a0*b0*c0*k2*cos((pi*m)/6) + a0*b0*c4*k0*cos((pi*t)/6) + a0*b0*c2*k0*cos((pi*t)/12) + a0*b4*c0*k0*cos((pi*theta)/90) + a0*b2*c0*k0*cos((pi*theta)/180) + a0*b0*c0*k3*sin((pi*m)/3) + a0*b0*c0*k1*sin((pi*m)/6) + a0*b0*c3*k0*sin((pi*t)/6) ...
I need to replace (a0b0c0k0 with p1) in the first part, (a1b0c0k0 with p2 and leave d) in the second part, (a0b0c0k4 with p3 and leave (cos(pi*m/3))) for the third part and so on.
Please I you can help me with regard.
2 Comments
Walter Roberson
on 17 Aug 2019
Are there any combinations of values that you want to replace in one part but not replace in the other? For example if you did
subs(y, a0*b0*c0*k0, p1)
then that would replace a0*b0*c0*k0 in all parts.
Answers (2)
Mohammed Ismail
on 17 Aug 2019
but the combinatin is not the same in all part. it has different combination and need to be replace and i have 137 part in my equation.
6 Comments
Walter Roberson
on 17 Aug 2019
The question is whether there is one place that you want to substitute p1 for a0*b0*c0*k0 but another place that you do not want to do that substituation for a0*b0*c0*k0 . If for example later there were a a0*b0*c0*k0*sin((pi*m)/17) then should the a0*b0*c0*k0 also be converted to p1 ?
The code is significantly easier to write if it is permitted to make the substitution anywhere it occurs. For example,
subs(y, {a0*b0*c0*k0, a1*b0*c0*k0, a0*b0*c0*k4}, {p1, p2, p3})
Mohammed Ismail
on 17 Aug 2019
Edited: Walter Roberson
on 18 Aug 2019
Thank you for your answer, but really the equation of mine is :
y = a0*b0*c0*k0 + a1*b0*c0*d*k0 + a0*b0*c0*k4*cos((pi*m)/3) + ...
a0*b0*c0*k2*cos((pi*m)/6) + a0*b0*c4*k0*cos((pi*t)/6) + ...
a0*b0*c2*k0*cos((pi*t)/12) + a0*b4*c0*k0*cos((pi*theta)/90) + ...
a0*b2*c0*k0*cos((pi*theta)/180) + a0*b0*c0*k3*sin((pi*m)/3) + ...
a0*b0*c0*k1*sin((pi*m)/6) + a0*b0*c3*k0*sin((pi*t)/6) + ...
a0*b0*c1*k0*sin((pi*t)/12) + a0*b3*c0*k0*sin((pi*theta)/90) + ...
a0*b1*c0*k0*sin((pi*theta)/180) + a2*b0*c0*d^2*k0 + ...
a0*b0*c3*k4*cos((pi*m)/3)*sin((pi*t)/6) + ...
a0*b0*c4*k3*cos((pi*t)/6)*sin((pi*m)/3) + ...
a0*b0*c3*k2*cos((pi*m)/6)*sin((pi*t)/6) + ...
a0*b0*c4*k1*cos((pi*t)/6)*sin((pi*m)/6) + ...
a0*b0*c1*k4*cos((pi*m)/3)*sin((pi*t)/12) + ...
a0*b0*c2*k3*cos((pi*t)/12)*sin((pi*m)/3) + ...
a0*b0*c1*k2*cos((pi*m)/6)*sin((pi*t)/12) + ...
a0*b0*c2*k1*cos((pi*t)/12)*sin((pi*m)/6) + ...
a0*b3*c0*k4*cos((pi*m)/3)*sin((pi*theta)/90) + ...
a0*b4*c0*k3*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a0*b3*c0*k2*cos((pi*m)/6)*sin((pi*theta)/90) + ...
a0*b4*c0*k1*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a0*b1*c0*k4*cos((pi*m)/3)*sin((pi*theta)/180) + ...
a0*b2*c0*k3*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a0*b1*c0*k2*cos((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b2*c0*k1*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a0*b3*c4*k0*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b4*c3*k0*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a0*b3*c2*k0*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b4*c1*k0*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a0*b1*c4*k0*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b2*c3*k0*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a0*b1*c2*k0*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b2*c1*k0*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a0*b0*c3*k3*sin((pi*m)/3)*sin((pi*t)/6) + ...
a0*b0*c3*k1*sin((pi*m)/6)*sin((pi*t)/6) + ...
a0*b0*c1*k3*sin((pi*m)/3)*sin((pi*t)/12) + ...
a0*b0*c1*k1*sin((pi*m)/6)*sin((pi*t)/12) + ...
a0*b3*c0*k3*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a0*b3*c0*k1*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a0*b1*c0*k3*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a0*b1*c0*k1*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b3*c3*k0*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c1*k0*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b1*c3*k0*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c1*k0*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b0*c0*d^2*k4*cos((pi*m)/3) + a2*b0*c0*d^2*k2*cos((pi*m)/6) + ...
a2*b0*c4*d^2*k0*cos((pi*t)/6) + a2*b0*c2*d^2*k0*cos((pi*t)/12) + ...
a2*b4*c0*d^2*k0*cos((pi*theta)/90) + a2*b2*c0*d^2*k0*cos((pi*theta)/180) ...
+ a2*b0*c0*d^2*k3*sin((pi*m)/3) + a2*b0*c0*d^2*k1*sin((pi*m)/6) + ...
a2*b0*c3*d^2*k0*sin((pi*t)/6) + a2*b0*c1*d^2*k0*sin((pi*t)/12) + ...
a2*b3*c0*d^2*k0*sin((pi*theta)/90) + a2*b1*c0*d^2*k0*sin((pi*theta)/180) ...
+ a1*b0*c0*d*k4*cos((pi*m)/3) + a1*b0*c0*d*k2*cos((pi*m)/6) + ...
a1*b0*c4*d*k0*cos((pi*t)/6) + a1*b0*c2*d*k0*cos((pi*t)/12) + ...
a1*b4*c0*d*k0*cos((pi*theta)/90) + a1*b2*c0*d*k0*cos((pi*theta)/180) + ...
a1*b0*c0*d*k3*sin((pi*m)/3) + a1*b0*c0*d*k1*sin((pi*m)/6) + ...
a1*b0*c3*d*k0*sin((pi*t)/6) + a1*b0*c1*d*k0*sin((pi*t)/12) + ...
a1*b3*c0*d*k0*sin((pi*theta)/90) + a1*b1*c0*d*k0*sin((pi*theta)/180) + ...
a0*b0*c4*k4*cos((pi*m)/3)*cos((pi*t)/6) + ...
a0*b0*c4*k2*cos((pi*m)/6)*cos((pi*t)/6) + ...
a0*b0*c2*k4*cos((pi*m)/3)*cos((pi*t)/12) + ...
a0*b0*c2*k2*cos((pi*m)/6)*cos((pi*t)/12) + ...
a0*b4*c0*k4*cos((pi*m)/3)*cos((pi*theta)/90) + ...
a0*b4*c0*k2*cos((pi*m)/6)*cos((pi*theta)/90) + ...
a0*b2*c0*k4*cos((pi*m)/3)*cos((pi*theta)/180) + ...
a0*b2*c0*k2*cos((pi*m)/6)*cos((pi*theta)/180) + ...
a0*b4*c4*k0*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a0*b4*c2*k0*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a0*b2*c4*k0*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a0*b2*c2*k0*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b0*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6) + ...
a1*b0*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6) + ...
a1*b0*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12) + ...
a1*b0*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12) + ...
a1*b4*c0*d*k4*cos((pi*m)/3)*cos((pi*theta)/90) + ...
a1*b4*c0*d*k2*cos((pi*m)/6)*cos((pi*theta)/90) + ...
a1*b2*c0*d*k4*cos((pi*m)/3)*cos((pi*theta)/180) + ...
a1*b2*c0*d*k2*cos((pi*m)/6)*cos((pi*theta)/180) + ...
a1*b4*c4*d*k0*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a1*b4*c2*d*k0*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a1*b2*c4*d*k0*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a1*b2*c2*d*k0*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b0*c3*d*k4*cos((pi*m)/3)*sin((pi*t)/6) + ...
a1*b0*c4*d*k3*cos((pi*t)/6)*sin((pi*m)/3) + ...
a1*b0*c3*d*k2*cos((pi*m)/6)*sin((pi*t)/6) + ...
a1*b0*c4*d*k1*cos((pi*t)/6)*sin((pi*m)/6) + ...
a1*b0*c1*d*k4*cos((pi*m)/3)*sin((pi*t)/12) + ...
a1*b0*c2*d*k3*cos((pi*t)/12)*sin((pi*m)/3) + ...
a1*b0*c1*d*k2*cos((pi*m)/6)*sin((pi*t)/12) + ...
a1*b0*c2*d*k1*cos((pi*t)/12)*sin((pi*m)/6) + ...
a1*b3*c0*d*k4*cos((pi*m)/3)*sin((pi*theta)/90) + ...
a1*b4*c0*d*k3*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a1*b3*c0*d*k2*cos((pi*m)/6)*sin((pi*theta)/90) + ...
a1*b4*c0*d*k1*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a1*b1*c0*d*k4*cos((pi*m)/3)*sin((pi*theta)/180) + ...
a1*b2*c0*d*k3*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a1*b1*c0*d*k2*cos((pi*m)/6)*sin((pi*theta)/180) + ...
a1*b2*c0*d*k1*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a1*b3*c4*d*k0*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k0*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a1*b3*c2*d*k0*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k0*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a1*b1*c4*d*k0*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k0*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a1*b1*c2*d*k0*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b2*c1*d*k0*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a1*b0*c3*d*k3*sin((pi*m)/3)*sin((pi*t)/6) + ...
a1*b0*c3*d*k1*sin((pi*m)/6)*sin((pi*t)/6) + ...
a1*b0*c1*d*k3*sin((pi*m)/3)*sin((pi*t)/12) + ...
a1*b0*c1*d*k1*sin((pi*m)/6)*sin((pi*t)/12) + ...
a1*b3*c0*d*k3*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a1*b3*c0*d*k1*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a1*b1*c0*d*k3*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a1*b1*c0*d*k1*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b4*c4*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a0*b4*c4*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a0*b4*c2*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a0*b4*c2*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a0*b2*c4*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a0*b2*c4*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a0*b2*c2*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a0*b2*c2*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b3*c3*d*k0*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c1*d*k0*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b1*c3*d*k0*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c1*d*k0*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b3*c4*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b4*c3*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a0*b4*c4*k3*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a0*b3*c4*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b4*c3*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a0*b4*c4*k1*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a0*b3*c2*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b4*c1*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a0*b4*c2*k3*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a0*b3*c2*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b4*c1*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a0*b4*c2*k1*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a0*b1*c4*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b2*c3*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a0*b2*c4*k3*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a0*b1*c4*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b2*c3*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a0*b2*c4*k1*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a0*b1*c2*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b2*c1*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a0*b2*c2*k3*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a0*b1*c2*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b2*c1*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a0*b2*c2*k1*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a0*b3*c3*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c4*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a0*b4*c3*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a0*b3*c3*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c4*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a0*b4*c3*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a0*b3*c1*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b3*c2*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a0*b4*c1*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a0*b3*c1*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b3*c2*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a0*b4*c1*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a0*b1*c3*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c4*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a0*b2*c3*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a0*b1*c3*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c4*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b2*c3*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a0*b1*c1*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b1*c2*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a0*b2*c1*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a0*b1*c1*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b1*c2*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b2*c1*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a0*b3*c3*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c3*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c1*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b3*c1*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b1*c3*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c3*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c1*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b1*c1*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b0*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6) + ...
a2*b0*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6) + ...
a2*b0*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12) + ...
a2*b0*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12) + ...
a2*b4*c0*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/90) + ...
a2*b4*c0*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/90) + ...
a2*b2*c0*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/180) + ...
a2*b2*c0*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/180) + ...
a2*b4*c4*d^2*k0*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a2*b4*c2*d^2*k0*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a2*b2*c4*d^2*k0*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a2*b2*c2*d^2*k0*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a2*b0*c3*d^2*k4*cos((pi*m)/3)*sin((pi*t)/6) + ...
a2*b0*c4*d^2*k3*cos((pi*t)/6)*sin((pi*m)/3) + ...
a2*b0*c3*d^2*k2*cos((pi*m)/6)*sin((pi*t)/6) + ...
a2*b0*c4*d^2*k1*cos((pi*t)/6)*sin((pi*m)/6) + ...
a2*b0*c1*d^2*k4*cos((pi*m)/3)*sin((pi*t)/12) + ...
a2*b0*c2*d^2*k3*cos((pi*t)/12)*sin((pi*m)/3) + ...
a2*b0*c1*d^2*k2*cos((pi*m)/6)*sin((pi*t)/12) + ...
a2*b0*c2*d^2*k1*cos((pi*t)/12)*sin((pi*m)/6) + ...
a2*b3*c0*d^2*k4*cos((pi*m)/3)*sin((pi*theta)/90) + ...
a2*b4*c0*d^2*k3*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a2*b3*c0*d^2*k2*cos((pi*m)/6)*sin((pi*theta)/90) + ...
a2*b4*c0*d^2*k1*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a2*b1*c0*d^2*k4*cos((pi*m)/3)*sin((pi*theta)/180) + ...
a2*b2*c0*d^2*k3*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a2*b1*c0*d^2*k2*cos((pi*m)/6)*sin((pi*theta)/180) + ...
a2*b2*c0*d^2*k1*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a2*b3*c4*d^2*k0*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k0*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a2*b3*c2*d^2*k0*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k0*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a2*b1*c4*d^2*k0*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k0*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a2*b1*c2*d^2*k0*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k0*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a2*b0*c3*d^2*k3*sin((pi*m)/3)*sin((pi*t)/6) + ...
a2*b0*c3*d^2*k1*sin((pi*m)/6)*sin((pi*t)/6) + ...
a2*b0*c1*d^2*k3*sin((pi*m)/3)*sin((pi*t)/12) + ...
a2*b0*c1*d^2*k1*sin((pi*m)/6)*sin((pi*t)/12) + ...
a2*b3*c0*d^2*k3*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a2*b3*c0*d^2*k1*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a2*b1*c0*d^2*k3*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a2*b1*c0*d^2*k1*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a2*b3*c3*d^2*k0*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c1*d^2*k0*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b1*c3*d^2*k0*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c1*d^2*k0*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b3*c3*d*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c4*d*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a1*b3*c3*d*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c4*d*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a1*b3*c1*d*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b3*c2*d*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a1*b3*c1*d*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b3*c2*d*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a1*b1*c3*d*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c4*d*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a1*b1*c3*d*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c4*d*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a1*b1*c1*d*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b1*c2*d*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a1*b2*c1*d*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a1*b1*c1*d*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b1*c2*d*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a1*b2*c1*d*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a1*b3*c3*d*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c3*d*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c1*d*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b3*c1*d*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b1*c3*d*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c3*d*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c1*d*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b1*c1*d*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b4*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a2*b4*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a2*b4*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a2*b4*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a2*b2*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a2*b2*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a2*b2*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a2*b2*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a2*b3*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a2*b4*c4*d^2*k3*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a2*b3*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a2*b4*c4*d^2*k1*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a2*b3*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a2*b4*c2*d^2*k3*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a2*b3*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a2*b4*c2*d^2*k1*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a2*b1*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a2*b2*c4*d^2*k3*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a2*b1*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a2*b2*c4*d^2*k1*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a2*b1*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a2*b2*c2*d^2*k3*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a2*b1*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a2*b2*c2*d^2*k1*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a2*b3*c3*d^2*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c4*d^2*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a2*b3*c3*d^2*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c4*d^2*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a2*b3*c1*d^2*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b3*c2*d^2*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a2*b3*c1*d^2*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b3*c2*d^2*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a2*b1*c3*d^2*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c4*d^2*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a2*b1*c3*d^2*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c4*d^2*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a2*b1*c1*d^2*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b1*c2*d^2*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a2*b1*c1*d^2*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b1*c2*d^2*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a2*b3*c3*d^2*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c3*d^2*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c1*d^2*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b3*c1*d^2*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b1*c3*d^2*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c3*d^2*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c1*d^2*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b1*c1*d^2*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b4*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a1*b4*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a1*b4*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a1*b4*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a1*b2*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a1*b2*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a1*b2*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b2*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b3*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a1*b4*c4*d*k3*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a1*b3*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a1*b4*c4*d*k1*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a1*b3*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a1*b4*c2*d*k3*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a1*b3*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a1*b4*c2*d*k1*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a1*b1*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a1*b2*c4*d*k3*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a1*b1*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a1*b2*c4*d*k1*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a1*b1*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b2*c1*d*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a1*b2*c2*d*k3*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a1*b1*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b2*c1*d*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a1*b2*c2*d*k1*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/6)
This is the problem.
Do you have any idea of loop or abbreviation?
Walter Roberson
on 18 Aug 2019
syms a0 a1 a2
syms b0 b1 b2 b3 b4
syms c0 c1 c2 c3 c4
syms k0 k1 k2 k3 k4
syms d m t theta
Y(1) = a0*b0*c0*k0 + a1*b0*c0*d*k0 + a0*b0*c0*k4*cos((pi*m)/3) + ...
a0*b0*c0*k2*cos((pi*m)/6) + a0*b0*c4*k0*cos((pi*t)/6) + ...
a0*b0*c2*k0*cos((pi*t)/12) + a0*b4*c0*k0*cos((pi*theta)/90) + ...
a0*b2*c0*k0*cos((pi*theta)/180) + a0*b0*c0*k3*sin((pi*m)/3) + ...
a0*b0*c0*k1*sin((pi*m)/6) + a0*b0*c3*k0*sin((pi*t)/6) + ...
a0*b0*c1*k0*sin((pi*t)/12) + a0*b3*c0*k0*sin((pi*theta)/90) + ...
a0*b1*c0*k0*sin((pi*theta)/180) + a2*b0*c0*d^2*k0 + ...
a0*b0*c3*k4*cos((pi*m)/3)*sin((pi*t)/6) + ...
a0*b0*c4*k3*cos((pi*t)/6)*sin((pi*m)/3) + ...
a0*b0*c3*k2*cos((pi*m)/6)*sin((pi*t)/6) + ...
a0*b0*c4*k1*cos((pi*t)/6)*sin((pi*m)/6) + ...
a0*b0*c1*k4*cos((pi*m)/3)*sin((pi*t)/12) + ...
a0*b0*c2*k3*cos((pi*t)/12)*sin((pi*m)/3) + ...
a0*b0*c1*k2*cos((pi*m)/6)*sin((pi*t)/12) + ...
a0*b0*c2*k1*cos((pi*t)/12)*sin((pi*m)/6) + ...
a0*b3*c0*k4*cos((pi*m)/3)*sin((pi*theta)/90) + ...
a0*b4*c0*k3*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a0*b3*c0*k2*cos((pi*m)/6)*sin((pi*theta)/90) + ...
a0*b4*c0*k1*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a0*b1*c0*k4*cos((pi*m)/3)*sin((pi*theta)/180) + ...
a0*b2*c0*k3*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a0*b1*c0*k2*cos((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b2*c0*k1*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a0*b3*c4*k0*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b4*c3*k0*cos((pi*theta)/90)*sin((pi*t)/6);
Y(2) = a0*b3*c2*k0*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b4*c1*k0*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a0*b1*c4*k0*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b2*c3*k0*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a0*b1*c2*k0*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b2*c1*k0*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a0*b0*c3*k3*sin((pi*m)/3)*sin((pi*t)/6) + ...
a0*b0*c3*k1*sin((pi*m)/6)*sin((pi*t)/6) + ...
a0*b0*c1*k3*sin((pi*m)/3)*sin((pi*t)/12) + ...
a0*b0*c1*k1*sin((pi*m)/6)*sin((pi*t)/12) + ...
a0*b3*c0*k3*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a0*b3*c0*k1*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a0*b1*c0*k3*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a0*b1*c0*k1*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b3*c3*k0*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c1*k0*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b1*c3*k0*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c1*k0*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b0*c0*d^2*k4*cos((pi*m)/3) + a2*b0*c0*d^2*k2*cos((pi*m)/6) + ...
a2*b0*c4*d^2*k0*cos((pi*t)/6) + a2*b0*c2*d^2*k0*cos((pi*t)/12) + ...
a2*b4*c0*d^2*k0*cos((pi*theta)/90) + a2*b2*c0*d^2*k0*cos((pi*theta)/180) ...
+ a2*b0*c0*d^2*k3*sin((pi*m)/3) + a2*b0*c0*d^2*k1*sin((pi*m)/6) + ...
a2*b0*c3*d^2*k0*sin((pi*t)/6) + a2*b0*c1*d^2*k0*sin((pi*t)/12) + ...
a2*b3*c0*d^2*k0*sin((pi*theta)/90) + a2*b1*c0*d^2*k0*sin((pi*theta)/180) ...
+ a1*b0*c0*d*k4*cos((pi*m)/3) + a1*b0*c0*d*k2*cos((pi*m)/6);
Y(3) = a1*b0*c4*d*k0*cos((pi*t)/6) + a1*b0*c2*d*k0*cos((pi*t)/12) + ...
a1*b4*c0*d*k0*cos((pi*theta)/90) + a1*b2*c0*d*k0*cos((pi*theta)/180) + ...
a1*b0*c0*d*k3*sin((pi*m)/3) + a1*b0*c0*d*k1*sin((pi*m)/6) + ...
a1*b0*c3*d*k0*sin((pi*t)/6) + a1*b0*c1*d*k0*sin((pi*t)/12) + ...
a1*b3*c0*d*k0*sin((pi*theta)/90) + a1*b1*c0*d*k0*sin((pi*theta)/180) + ...
a0*b0*c4*k4*cos((pi*m)/3)*cos((pi*t)/6) + ...
a0*b0*c4*k2*cos((pi*m)/6)*cos((pi*t)/6) + ...
a0*b0*c2*k4*cos((pi*m)/3)*cos((pi*t)/12) + ...
a0*b0*c2*k2*cos((pi*m)/6)*cos((pi*t)/12) + ...
a0*b4*c0*k4*cos((pi*m)/3)*cos((pi*theta)/90) + ...
a0*b4*c0*k2*cos((pi*m)/6)*cos((pi*theta)/90) + ...
a0*b2*c0*k4*cos((pi*m)/3)*cos((pi*theta)/180) + ...
a0*b2*c0*k2*cos((pi*m)/6)*cos((pi*theta)/180) + ...
a0*b4*c4*k0*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a0*b4*c2*k0*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a0*b2*c4*k0*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a0*b2*c2*k0*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b0*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6) + ...
a1*b0*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6) + ...
a1*b0*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12) + ...
a1*b0*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12) + ...
a1*b4*c0*d*k4*cos((pi*m)/3)*cos((pi*theta)/90) + ...
a1*b4*c0*d*k2*cos((pi*m)/6)*cos((pi*theta)/90) + ...
a1*b2*c0*d*k4*cos((pi*m)/3)*cos((pi*theta)/180) + ...
a1*b2*c0*d*k2*cos((pi*m)/6)*cos((pi*theta)/180);
Y(4) = a1*b4*c4*d*k0*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a1*b4*c2*d*k0*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a1*b2*c4*d*k0*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a1*b2*c2*d*k0*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b0*c3*d*k4*cos((pi*m)/3)*sin((pi*t)/6) + ...
a1*b0*c4*d*k3*cos((pi*t)/6)*sin((pi*m)/3) + ...
a1*b0*c3*d*k2*cos((pi*m)/6)*sin((pi*t)/6) + ...
a1*b0*c4*d*k1*cos((pi*t)/6)*sin((pi*m)/6) + ...
a1*b0*c1*d*k4*cos((pi*m)/3)*sin((pi*t)/12) + ...
a1*b0*c2*d*k3*cos((pi*t)/12)*sin((pi*m)/3) + ...
a1*b0*c1*d*k2*cos((pi*m)/6)*sin((pi*t)/12) + ...
a1*b0*c2*d*k1*cos((pi*t)/12)*sin((pi*m)/6) + ...
a1*b3*c0*d*k4*cos((pi*m)/3)*sin((pi*theta)/90) + ...
a1*b4*c0*d*k3*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a1*b3*c0*d*k2*cos((pi*m)/6)*sin((pi*theta)/90) + ...
a1*b4*c0*d*k1*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a1*b1*c0*d*k4*cos((pi*m)/3)*sin((pi*theta)/180) + ...
a1*b2*c0*d*k3*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a1*b1*c0*d*k2*cos((pi*m)/6)*sin((pi*theta)/180) + ...
a1*b2*c0*d*k1*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a1*b3*c4*d*k0*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k0*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a1*b3*c2*d*k0*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k0*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a1*b1*c4*d*k0*cos((pi*t)/6)*sin((pi*theta)/180);
Y(5) = a1*b2*c3*d*k0*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a1*b1*c2*d*k0*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b2*c1*d*k0*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a1*b0*c3*d*k3*sin((pi*m)/3)*sin((pi*t)/6) + ...
a1*b0*c3*d*k1*sin((pi*m)/6)*sin((pi*t)/6) + ...
a1*b0*c1*d*k3*sin((pi*m)/3)*sin((pi*t)/12) + ...
a1*b0*c1*d*k1*sin((pi*m)/6)*sin((pi*t)/12) + ...
a1*b3*c0*d*k3*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a1*b3*c0*d*k1*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a1*b1*c0*d*k3*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a1*b1*c0*d*k1*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b4*c4*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a0*b4*c4*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a0*b4*c2*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a0*b4*c2*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a0*b2*c4*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a0*b2*c4*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a0*b2*c2*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a0*b2*c2*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b3*c3*d*k0*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c1*d*k0*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b1*c3*d*k0*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c1*d*k0*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b3*c4*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b4*c3*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/6);
Y(6) = a0*b4*c4*k3*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a0*b3*c4*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b4*c3*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a0*b4*c4*k1*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a0*b3*c2*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b4*c1*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a0*b4*c2*k3*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a0*b3*c2*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b4*c1*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a0*b4*c2*k1*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a0*b1*c4*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b2*c3*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a0*b2*c4*k3*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a0*b1*c4*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b2*c3*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a0*b2*c4*k1*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a0*b1*c2*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b2*c1*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a0*b2*c2*k3*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a0*b1*c2*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b2*c1*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a0*b2*c2*k1*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a0*b3*c3*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c4*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a0*b4*c3*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/6);
Y(7) = a0*b3*c3*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c4*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a0*b4*c3*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a0*b3*c1*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b3*c2*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a0*b4*c1*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a0*b3*c1*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b3*c2*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a0*b4*c1*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a0*b1*c3*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c4*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a0*b2*c3*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a0*b1*c3*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c4*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b2*c3*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a0*b1*c1*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b1*c2*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a0*b2*c1*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a0*b1*c1*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b1*c2*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a0*b2*c1*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a0*b3*c3*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c3*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a0*b3*c1*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a0*b3*c1*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90);
Y(8) = a0*b1*c3*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c3*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a0*b1*c1*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a0*b1*c1*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b0*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6) + ...
a2*b0*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6) + ...
a2*b0*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12) + ...
a2*b0*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12) + ...
a2*b4*c0*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/90) + ...
a2*b4*c0*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/90) + ...
a2*b2*c0*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/180) + ...
a2*b2*c0*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/180) + ...
a2*b4*c4*d^2*k0*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a2*b4*c2*d^2*k0*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a2*b2*c4*d^2*k0*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a2*b2*c2*d^2*k0*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a2*b0*c3*d^2*k4*cos((pi*m)/3)*sin((pi*t)/6) + ...
a2*b0*c4*d^2*k3*cos((pi*t)/6)*sin((pi*m)/3) + ...
a2*b0*c3*d^2*k2*cos((pi*m)/6)*sin((pi*t)/6) + ...
a2*b0*c4*d^2*k1*cos((pi*t)/6)*sin((pi*m)/6) + ...
a2*b0*c1*d^2*k4*cos((pi*m)/3)*sin((pi*t)/12) + ...
a2*b0*c2*d^2*k3*cos((pi*t)/12)*sin((pi*m)/3) + ...
a2*b0*c1*d^2*k2*cos((pi*m)/6)*sin((pi*t)/12) + ...
a2*b0*c2*d^2*k1*cos((pi*t)/12)*sin((pi*m)/6) + ...
a2*b3*c0*d^2*k4*cos((pi*m)/3)*sin((pi*theta)/90);
Y(9) = a2*b4*c0*d^2*k3*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a2*b3*c0*d^2*k2*cos((pi*m)/6)*sin((pi*theta)/90) + ...
a2*b4*c0*d^2*k1*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a2*b1*c0*d^2*k4*cos((pi*m)/3)*sin((pi*theta)/180) + ...
a2*b2*c0*d^2*k3*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a2*b1*c0*d^2*k2*cos((pi*m)/6)*sin((pi*theta)/180) + ...
a2*b2*c0*d^2*k1*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a2*b3*c4*d^2*k0*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k0*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a2*b3*c2*d^2*k0*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k0*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a2*b1*c4*d^2*k0*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k0*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a2*b1*c2*d^2*k0*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k0*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a2*b0*c3*d^2*k3*sin((pi*m)/3)*sin((pi*t)/6) + ...
a2*b0*c3*d^2*k1*sin((pi*m)/6)*sin((pi*t)/6) + ...
a2*b0*c1*d^2*k3*sin((pi*m)/3)*sin((pi*t)/12) + ...
a2*b0*c1*d^2*k1*sin((pi*m)/6)*sin((pi*t)/12) + ...
a2*b3*c0*d^2*k3*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a2*b3*c0*d^2*k1*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a2*b1*c0*d^2*k3*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a2*b1*c0*d^2*k1*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a2*b3*c3*d^2*k0*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c1*d^2*k0*sin((pi*t)/12)*sin((pi*theta)/90);
Y(10) = a2*b1*c3*d^2*k0*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c1*d^2*k0*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b3*c3*d*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c4*d*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a1*b3*c3*d*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c4*d*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a1*b3*c1*d*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b3*c2*d*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a1*b3*c1*d*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b3*c2*d*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a1*b1*c3*d*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c4*d*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a1*b1*c3*d*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c4*d*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a1*b1*c1*d*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b1*c2*d*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a1*b2*c1*d*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a1*b1*c1*d*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b1*c2*d*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/180);
Y(11) = a1*b2*c1*d*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a1*b3*c3*d*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c3*d*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b3*c1*d*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b3*c1*d*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b1*c3*d*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c3*d*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b1*c1*d*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b1*c1*d*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b4*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a2*b4*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a2*b4*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a2*b4*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a2*b2*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a2*b2*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a2*b2*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a2*b2*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a2*b3*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a2*b4*c4*d^2*k3*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a2*b3*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a2*b4*c4*d^2*k1*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a2*b3*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/12);
Y(12) = a2*b4*c2*d^2*k3*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a2*b3*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a2*b4*c2*d^2*k1*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a2*b1*c4*d^2*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a2*b2*c4*d^2*k3*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a2*b1*c4*d^2*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a2*b2*c4*d^2*k1*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a2*b1*c2*d^2*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a2*b2*c2*d^2*k3*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a2*b1*c2*d^2*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a2*b2*c2*d^2*k1*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a2*b3*c3*d^2*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c4*d^2*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a2*b3*c3*d^2*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c4*d^2*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a2*b4*c3*d^2*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a2*b3*c1*d^2*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b3*c2*d^2*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k3*cos((pi*theta)/90)*sin((pi*m)/3)*sin((pi*t)/12);
Y(13) = a2*b3*c1*d^2*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b3*c2*d^2*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/90) + ...
a2*b4*c1*d^2*k1*cos((pi*theta)/90)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a2*b1*c3*d^2*k4*cos((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c4*d^2*k3*cos((pi*t)/6)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/6) + ...
a2*b1*c3*d^2*k2*cos((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c4*d^2*k1*cos((pi*t)/6)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a2*b2*c3*d^2*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/6) + ...
a2*b1*c1*d^2*k4*cos((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b1*c2*d^2*k3*cos((pi*t)/12)*sin((pi*m)/3)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k3*cos((pi*theta)/180)*sin((pi*m)/3)*sin((pi*t)/12) + ...
a2*b1*c1*d^2*k2*cos((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b1*c2*d^2*k1*cos((pi*t)/12)*sin((pi*m)/6)*sin((pi*theta)/180) + ...
a2*b2*c1*d^2*k1*cos((pi*theta)/180)*sin((pi*m)/6)*sin((pi*t)/12) + ...
a2*b3*c3*d^2*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c3*d^2*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/90) + ...
a2*b3*c1*d^2*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b3*c1*d^2*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/90) + ...
a2*b1*c3*d^2*k3*sin((pi*m)/3)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c3*d^2*k1*sin((pi*m)/6)*sin((pi*t)/6)*sin((pi*theta)/180) + ...
a2*b1*c1*d^2*k3*sin((pi*m)/3)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a2*b1*c1*d^2*k1*sin((pi*m)/6)*sin((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b4*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/90) + ...
a1*b4*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/90);
Y(14) = a1*b4*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a1*b4*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/90) + ...
a1*b2*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a1*b2*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6)*cos((pi*theta)/180) + ...
a1*b2*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b2*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12)*cos((pi*theta)/180) + ...
a1*b3*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a1*b4*c4*d*k3*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a1*b3*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/90) + ...
a1*b4*c3*d*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/6) + ...
a1*b4*c4*d*k1*cos((pi*t)/6)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a1*b3*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k4*cos((pi*m)/3)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a1*b4*c2*d*k3*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/3) + ...
a1*b3*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/90) + ...
a1*b4*c1*d*k2*cos((pi*m)/6)*cos((pi*theta)/90)*sin((pi*t)/12) + ...
a1*b4*c2*d*k1*cos((pi*t)/12)*cos((pi*theta)/90)*sin((pi*m)/6) + ...
a1*b1*c4*d*k4*cos((pi*m)/3)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a1*b2*c4*d*k3*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a1*b1*c4*d*k2*cos((pi*m)/6)*cos((pi*t)/6)*sin((pi*theta)/180) + ...
a1*b2*c3*d*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/6) + ...
a1*b2*c4*d*k1*cos((pi*t)/6)*cos((pi*theta)/180)*sin((pi*m)/6) + ...
a1*b1*c2*d*k4*cos((pi*m)/3)*cos((pi*t)/12)*sin((pi*theta)/180);
Y(15) = a1*b2*c1*d*k4*cos((pi*m)/3)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a1*b2*c2*d*k3*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/3) + ...
a1*b1*c2*d*k2*cos((pi*m)/6)*cos((pi*t)/12)*sin((pi*theta)/180) + ...
a1*b2*c1*d*k2*cos((pi*m)/6)*cos((pi*theta)/180)*sin((pi*t)/12) + ...
a1*b2*c2*d*k1*cos((pi*t)/12)*cos((pi*theta)/180)*sin((pi*m)/6);
y = sum(Y);
Walter Roberson
on 18 Aug 2019
You have 375 terms. It is not clear what you want substituted.
If you want to convert all of the a b c d k combinations into individual variables, then toss the above code into a script a0b0c0k0.m and do
sin = @(expr) 1; cos = @(expr) 1;
a0b0c0k0;
ch = children(y);
p = sym('p', [1, length(ch)]);
clear sin cos
a0b0c0k0;
transformed = collect(subs(y, ch, p),p);
The result will start out
p1 + sin((m*pi)/6)*p2 + sin((t*pi)/12)*p3 + sin((theta*pi)/180)*p4 + cos((m*pi)/6)*p5
Mohammed Ismail
on 19 Aug 2019
thank u for y kinds. yes, i want to convert all the combination of a b c k and leave the sin cos and d.
Mohammed Ismail
on 21 Aug 2019
thank you for your assistant. Can i ask you some questions.
First, the benifit of using this (sin = @(expr) 1).
Second, when i put in the (sin = @(expr) 1) an (theta*pi/3) the stop and not take it.
1 Comment
Walter Roberson
on 21 Aug 2019
You have a lot of terms involving multiplying several variables together and then multiply by a sin or cos term. One way to extract out the variables without the sin and cos terms is to temporarily redefine sin and cos as returning 1 so that
a0*b2*c1*d*k2*cos(theta*pi/3)
would evaluate as
a0*b2*c1*d*k2*1
Which is
a0*b2*c1*d*k2
It is a trick for stripping off the sin and cos terms leaving just the variables.
See Also
Categories
Find more on Function Creation in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!An Error Occurred
Unable to complete the action because of changes made to the page. Reload the page to see its updated state.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
You can also select a web site from the following list
How to Get Best Site Performance
Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Asia Pacific
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)