Voor matlab moet je ook alles uitkauwen, hij kan toch zelf ook wel een beetje z'n best doen
T =
1/8*l^2*sin(-gam+alp)^2*alpd^2*m+1/8*l^2*cos(-gam+alp)^2*alpd^2*m+1/24*alpd^2*m*l^2+1/2*l^2*sin(-gam+alp)^2*alpd^2*M+1/2*l^2*cos(-gam+alp)^2*alpd^2*M+(-1/2*l*sin(-gam+alp)*alpd+1/4*l*sin(gam+bet)*betd)*m*(-l*sin(-gam+alp)*alpd+1/2*l*sin(gam+bet)*betd)+(-1/2*l*cos(-gam+alp)*alpd-1/4*l*cos(gam+bet)*betd)*m*(-l*cos(-gam+alp)*alpd-1/2*l*cos(gam+bet)*betd)+1/24*betd^2*m*l^2
>> simlpe(T)
??? Undefined function or method 'simlpe' for input arguments of type 'sym'.
>> simple(T)
simplify:
1/6*l^2*(4*m*alpd^2+3*alpd^2*M-3*m*sin(-gam+alp)*alpd*sin(gam+bet)*betd+betd^2*m+3*m*cos(-gam+alp)*alpd*cos(gam+bet)*betd)
radsimp:
5/8*l^2*sin(-gam+alp)^2*alpd^2*m+5/8*l^2*cos(-gam+alp)^2*alpd^2*m+1/24*m*l^2*alpd^2+1/2*l^2*sin(-gam+alp)^2*alpd^2*M+1/2*l^2*cos(-gam+alp)^2*alpd^2*M-1/2*m*l^2*sin(-gam+alp)*alpd*sin(gam+bet)*betd+1/8*m*l^2*sin(gam+bet)^2*betd^2+1/2*m*l^2*cos(-gam+alp)*alpd*cos(gam+bet)*betd+1/8*m*l^2*cos(gam+bet)^2*betd^2+1/24*betd^2*m*l^2
combine(trig):
2/3*m*l^2*alpd^2+1/2*l^2*alpd^2*M+1/2*m*l^2*alpd*betd*cos(alp+bet)+1/6*betd^2*m*l^2
factor:
1/24*l^2*(15*sin(-gam+alp)^2*alpd^2*m+15*cos(-gam+alp)^2*alpd^2*m+m*alpd^2+12*sin(-gam+alp)^2*alpd^2*M+12*cos(-gam+alp)^2*alpd^2*M-12*m*sin(-gam+alp)*alpd*sin(gam+bet)*betd+3*m*sin(gam+bet)^2*betd^2+12*m*cos(-gam+alp)*alpd*cos(gam+bet)*betd+3*m*cos(gam+bet)^2*betd^2+betd^2*m)
expand:
1/24*m*l^2*alpd^2+5/8*m*l^2*alpd^2*cos(gam)^2*sin(alp)^2+5/8*m*l^2*alpd^2*sin(gam)^2*cos(alp)^2+5/8*m*l^2*alpd^2*cos(gam)^2*cos(alp)^2+5/8*m*l^2*alpd^2*sin(gam)^2*sin(alp)^2+1/2*l^2*alpd^2*M*sin(gam)^2*cos(alp)^2+1/2*l^2*alpd^2*M*cos(gam)^2*sin(alp)^2+1/2*l^2*alpd^2*M*cos(gam)^2*cos(alp)^2+1/2*l^2*alpd^2*M*sin(gam)^2*sin(alp)^2+1/2*m*l^2*alpd*sin(gam)^2*cos(alp)*betd*cos(bet)-1/2*m*l^2*alpd*cos(gam)^2*sin(alp)*betd*sin(bet)+1/8*m*l^2*betd^2*sin(gam)^2*cos(bet)^2+1/8*m*l^2*betd^2*cos(gam)^2*sin(bet)^2+1/2*m*l^2*alpd*cos(gam)^2*cos(alp)*betd*cos(bet)-1/2*m*l^2*alpd*sin(gam)^2*sin(alp)*betd*sin(bet)+1/8*m*l^2*betd^2*cos(gam)^2*cos(bet)^2+1/8*m*l^2*betd^2*sin(gam)^2*sin(bet)^2+1/24*betd^2*m*l^2
combine:
2/3*m*l^2*alpd^2+1/2*l^2*alpd^2*M+1/2*m*l^2*alpd*betd*cos(alp+bet)+1/6*betd^2*m*l^2
convert(exp):
-1/32*l^2*(exp(i*(-gam+alp))-1/exp(i*(-gam+alp)))^2*alpd^2*m+1/8*l^2*(1/2*exp(i*(-gam+alp))+1/2/exp(i*(-gam+alp)))^2*alpd^2*m+1/24*m*l^2*alpd^2-1/8*l^2*(exp(i*(-gam+alp))-1/exp(i*(-gam+alp)))^2*alpd^2*M+1/2*l^2*(1/2*exp(i*(-gam+alp))+1/2/exp(i*(-gam+alp)))^2*alpd^2*M+(1/4*i*l*(exp(i*(-gam+alp))-1/exp(i*(-gam+alp)))*alpd-1/8*i*l*(exp(i*(gam+bet))-1/exp(i*(gam+bet)))*betd)*m*(1/2*i*l*(exp(i*(-gam+alp))-1/exp(i*(-gam+alp)))*alpd-1/4*i*l*(exp(i*(gam+bet))-1/exp(i*(gam+bet)))*betd)+(-1/2*l*(1/2*exp(i*(-gam+alp))+1/2/exp(i*(-gam+alp)))*alpd-1/4*l*(1/2*exp(i*(gam+bet))+1/2/exp(i*(gam+bet)))*betd)*m*(-l*(1/2*exp(i*(-gam+alp))+1/2/exp(i*(-gam+alp)))*alpd-1/2*l*(1/2*exp(i*(gam+bet))+1/2/exp(i*(gam+bet)))*betd)+1/24*betd^2*m*l^2
convert(sincos):
1/8*l^2*sin(-gam+alp)^2*alpd^2*m+1/8*l^2*cos(-gam+alp)^2*alpd^2*m+1/24*m*l^2*alpd^2+1/2*l^2*sin(-gam+alp)^2*alpd^2*M+1/2*l^2*cos(-gam+alp)^2*alpd^2*M+(-1/2*l*sin(-gam+alp)*alpd+1/4*l*sin(gam+bet)*betd)*m*(-l*sin(-gam+alp)*alpd+1/2*l*sin(gam+bet)*betd)+(-1/2*l*cos(-gam+alp)*alpd-1/4*l*cos(gam+bet)*betd)*m*(-l*cos(-gam+alp)*alpd-1/2*l*cos(gam+bet)*betd)+1/24*betd^2*m*l^2
convert(tan):
1/2*l^2*tan(-1/2*gam+1/2*alp)^2/(1+tan(-1/2*gam+1/2*alp)^2)^2*alpd^2*m+1/8*l^2*(1-tan(-1/2*gam+1/2*alp)^2)^2/(1+tan(-1/2*gam+1/2*alp)^2)^2*alpd^2*m+1/24*m*l^2*alpd^2+2*l^2*tan(-1/2*gam+1/2*alp)^2/(1+tan(-1/2*gam+1/2*alp)^2)^2*alpd^2*M+1/2*l^2*(1-tan(-1/2*gam+1/2*alp)^2)^2/(1+tan(-1/2*gam+1/2*alp)^2)^2*alpd^2*M+(-l*tan(-1/2*gam+1/2*alp)/(1+tan(-1/2*gam+1/2*alp)^2)*alpd+1/2*l*tan(1/2*gam+1/2*bet)/(1+tan(1/2*gam+1/2*bet)^2)*betd)*m*(-2*l*tan(-1/2*gam+1/2*alp)/(1+tan(-1/2*gam+1/2*alp)^2)*alpd+l*tan(1/2*gam+1/2*bet)/(1+tan(1/2*gam+1/2*bet)^2)*betd)+(-1/2*l*(1-tan(-1/2*gam+1/2*alp)^2)/(1+tan(-1/2*gam+1/2*alp)^2)*alpd-1/4*l*(1-tan(1/2*gam+1/2*bet)^2)/(1+tan(1/2*gam+1/2*bet)^2)*betd)*m*(-l*(1-tan(-1/2*gam+1/2*alp)^2)/(1+tan(-1/2*gam+1/2*alp)^2)*alpd-1/2*l*(1-tan(1/2*gam+1/2*bet)^2)/(1+tan(1/2*gam+1/2*bet)^2)*betd)+1/24*betd^2*m*l^2
collect(m):
(1/8*l^2*cos(-gam+alp)^2*alpd^2+1/24*alpd^2*l^2+1/8*l^2*sin(-gam+alp)^2*alpd^2+(-1/2*l*sin(-gam+alp)*alpd+1/4*l*sin(gam+bet)*betd)*(-l*sin(-gam+alp)*alpd+1/2*l*sin(gam+bet)*betd)+(-1/2*l*cos(-gam+alp)*alpd-1/4*l*cos(gam+bet)*betd)*(-l*cos(-gam+alp)*alpd-1/2*l*cos(gam+bet)*betd)+1/24*betd^2*l^2)*m+1/2*l^2*cos(-gam+alp)^2*alpd^2*M+1/2*l^2*sin(-gam+alp)^2*alpd^2*M
mwcos2sin:
1/8*l^2*sin(-gam+alp)^2*alpd^2*m+1/8*l^2*(1-sin(-gam+alp)^2)*alpd^2*m+1/24*m*l^2*alpd^2+1/2*l^2*sin(-gam+alp)^2*alpd^2*M+1/2*l^2*(1-sin(-gam+alp)^2)*alpd^2*M+(-1/2*l*sin(-gam+alp)*alpd+1/4*l*sin(gam+bet)*betd)*m*(-l*sin(-gam+alp)*alpd+1/2*l*sin(gam+bet)*betd)+(-1/2*l*cos(-gam+alp)*alpd-1/4*l*cos(gam+bet)*betd)*m*(-l*cos(-gam+alp)*alpd-1/2*l*cos(gam+bet)*betd)+1/24*betd^2*m*l^2
ans =
2/3*m*l^2*alpd^2+1/2*l^2*alpd^2*M+1/2*m*l^2*alpd*betd*cos(alp+bet)+1/6*betd^2*m*l^2