On 5/15/07, Ryan Krauss <ryanlists at gmail.com> wrote:
>
> cos(beta)*cosh(beta)*L*c*s+cos(beta)*cosh(beta)*L*k
>
> I don't seem to be able to make Maxima factor/simplify this as
> cos(beta)*cosh(beta)*L*(c*s+k)
>
factorsum will take the first expression to the second.
However, factorsum only operates on the numerator. To apply to both the
numerator and the denominator, the simplest way is map(factorsum,test).
You can also control the application of factor and factorsum in more detail
using substpart (qv). Here's a little transcript (using very wide lines so
the structure will be clear)...
test =>
3
L (cos(beta) sinh(beta) c s L - cosh(beta) sin(beta) c s L + cos(beta)
sinh(beta) k L - cosh(beta) sin(beta) k L - 2 beta sin(beta) sinh(beta) EI)
---------------------------------------------------------------------------------------------------------------------------------------------------
3
beta EI (cos(beta) cosh(beta) c s L + c s L + cos(beta) cosh(beta) k
L + k L + beta cos(beta) sinh(beta) EI - beta cosh(beta) sin(beta) EI)
map(factorsum,%) =>
3
L (cos(beta) sinh(beta) c s L - cosh(beta) sin(beta) c s L + cos(beta)
sinh(beta) k L - cosh(beta) sin(beta) k L - 2 beta sin(beta) sinh(beta) EI)
---------------------------------------------------------------------------------------------------------------------------------------------------
3
beta EI ((cos(beta) cosh(beta) + 1) (c s + k) L +
beta (cos(beta) sinh(beta) - cosh(beta) sin(beta)) EI)
substpart(factor(piece),%,1,2,[1,2,3,4]) =>
3
L ((cos(beta) sinh(beta) - cosh(beta) sin(beta)) (c s + k) L - 2
beta sin(beta) sinh(beta) EI)
---------------------------------------------------------------------------------------------------------
3
beta EI ((cos(beta) cosh(beta) + 1) (c s + k) L + beta (cos(beta)
sinh(beta) - cosh(beta) sin(beta)) EI)
Alternatively:
testfs1: substpart(factorsum(piece),testfs,1,2,[1,3,5]);
3
L (sinh(beta) (cos(beta) (c s + k) L - 2 beta sin(beta) EI) -
cosh(beta) sin(beta) c s L - cosh(beta) sin(beta) k L)
---------------------------------------------------------------------------------------------------------------------
3
beta EI ((cos(beta) cosh(beta) + 1) (c s + k) L + beta (cos(beta)
sinh(beta) - cosh(beta) sin(beta)) EI)
substpart(factorsum(piece),testfs1,1,2,[2,3]);
3
L (sinh(beta) (cos(beta) (c s + k) L - 2 beta sin(beta) EI) -
cosh(beta) sin(beta) (c s + k) L)
---------------------------------------------------------------------------------------------------------
3
beta EI ((cos(beta) cosh(beta) + 1) (c s + k) L + beta (cos(beta)
sinh(beta) - cosh(beta) sin(beta)) EI)
Substpart is worth getting comfortable with if you want to do fine-tuning of
the form of an expression.
-s