(%i1) e1 : j0 + m1· d = jc ;
(e1) d m1 + j0 = jc
(%i2) se1 : solve( e1, d) ;
(se1) [ d = jc j0 m1 ]
(%i3) d : rhs( se1[ 1]) ;
(d) jc j0 m1
(%i4) j1 : jc + m2·( 1 d) ;
(j1) ( 1 jc j0 m1 ) m2 + jc
(%i5) j1 : ratsimp( j1) ;
(j1) ( m1 jc + j0 ) m2 + jc m1 m1
(%i6) j1 : factor( j1) ;
(j1) m1 m2 jc m2 + j0 m2 + jc m1 m1
(%i7) e2 : j1 = j0 ;
(e2) m1 m2 jc m2 + j0 m2 + jc m1 m1 = j0
(%i8) se2 : solve( e2, j0) ;
(se2) [ j0 = ( m1 jc ) m2 + jc m1 m2 m1 ]
(%i9) J0 : rhs( se2[ 1]) ;
(J0) ( m1 jc ) m2 + jc m1 m2 m1
(%i10) J0 : factor( J0) ;
(J0) m1 m2 jc m2 + jc m1 m2 m1
(%i11) D0 : ev( d, j0 = J0) ;
(D0) m1 m2 jc m2 + jc m1 m2 m1 + jc m1
(%i12) D0 : ratsimp( D0) ;
(D0) m2 m2 m1
(%i13) Jav :( J0 + jc) / 2 ;
(Jav) jc m1 m2 jc m2 + jc m1 m2 m1 2
(%i14) Jav : ratsimp( Jav) ;
(Jav) ( m1 2 jc ) m2 + 2 jc m1 2 m2 2 m1
(%i15) Jav : factor( Jav) ;
(Jav) m1 m2 2 jc m2 + 2 jc m1 2 ( m2 m1 )
(%i16) dJav : Jav jc ;
(dJav) m1 m2 2 jc m2 + 2 jc m1 2 ( m2 m1 ) jc
(%i17) dJav : ratsimp( dJav) ;
(dJav) m1 m2 2 m2 2 m1
(%i18) dJav : factor( dJav) ;
(dJav) m1 m2 2 ( m2 m1 )
(%i19) Jav : jc + dJav ;
(Jav) jc m1 m2 2 ( m2 m1 )
(%i20) djc : diff( Jav, jc) ;
(djc) 1
(%i21) dm1 : diff( Jav, m1) ;
(dm1) m2 2 ( m2 m1 ) m1 m2 2 ( m2 m1 ) 2
(%i22) dm1 : ratsimp( dm1) ;
(dm1) m2 2 2 m2 2 4 m1 m2 + 2 m1 2
(%i23) dm1 : factor( dm1) ;
(dm1) m2 2 2 ( m2 m1 ) 2
(%i24) dm2 : diff( Jav, m2) ;
(dm2) m1 m2 2 ( m2 m1 ) 2 m1 2 ( m2 m1 )
(%i25) dm2 : ratsimp( dm2) ;
(dm2) m1 2 2 m2 2 4 m1 m2 + 2 m1 2
(%i26) dm2 : factor( dm2) ;
(dm2) m1 2 2 ( m2 m1 ) 2
(%i27) D0 ;
(%o27) m2 m2 m1
(%i28) e3 : D = D0 ;
(e3) D = m2 m2 m1
(%i29) se3 : solve( e3, m2) ;
(se3) [ m2 = D m1 D 1 ]
(%i30) M2 : rhs( se3[ 1]) ;
(M2) D m1 D 1
(%i31) dM1 : ev( dm1, m2 = M2) ;
(dM1) D 2 m1 2 2 ( D 1 ) 2 ( D m1 D 1 m1 ) 2
(%i32) dM1 : ratsimp( dM1) ;
(dM1) D 2 2
(%i33) dM2 : ev( dm2, m2 = M2) ;
(dM2) m1 2 2 ( D m1 D 1 m1 ) 2
(%i34) dM2 : ratsimp( dM2) ;
(dM2) D 2 2 D + 1 2
(%i35) dM2 : factor( dM2) ;
(dM2) ( D 1 ) 2 2

Created with wxMaxima.