From: katherle@daisy.uwaterloo.ca (Kate Atherley)
Subject: Re: Maple Problem : Error, (in normal/doit) floats not handled
Organization: University of Waterloo
Date: Tue, 15 Sep 1992 14:54:01 GMT

In article <1992Sep13.143621.1@ntuvax.ntu.ac.sg> mc1586217@ntuvax.ntu.ac.sg writes:
>The following maple program was ran under DEC VMS and the last line gave the
>following error :-
>
>Error, (in normal/doit) floats not handled
:
:
>
>Does anyone know what's the problem. Thank you.
>

This is a bug in the current version of Maple.  This has been fixed in
the research version, and you can expect to see the fix in the 
next release of Maple.  I include a session with the result of the
calculation.

Unfortunately, there is no patch for the current version.

    |\^/|
._|\|   |/|_. Symbolic Computation Group, University of Waterloo
 \  MAPLE  /  Research Version --- Mon Sep 14 16:24:13 EDT 1992
 <____ ____>  Type ? for help.
      |
> with(linalg):
Warning: new definition for   norm
Warning: new definition for   trace
> Q:=array(symmetric,1..3,1..3,[(1,1)=20.0,(1,2)=0.7,(1,3)=0.0,(2,2)=2.0,\
> (2,3)=0.0,(3,3)=0.7]);
                                  [ 20.0   .7   0 ]
                                  [               ]
                             Q := [  .7   2.0   0 ]
                                  [               ]
                                  [   0    0   .7 ]

> m:=cos(theta);
                                m := cos(theta)

> n:=sin(theta);
                                n := sin(theta)

> T:=array(1..3,1..3,[(1,1)=m^2,(1,2)=n^2,(1,3)=2*m*n,(2,1)=n^2\
> ,(2,2)=m^2,(2,3)=-(2*m*n),(3,1)=-(m*n),(3,2)=m*n,(3,3)=(m^2-n^2)]);
               [           2            2                            ]
               [ cos(theta)   sin(theta)              2 %1           ]
               [                                                     ]
               [           2            2                            ]
          T := [ sin(theta)   cos(theta)             - 2 %1          ]
               [                                                     ]
               [                                     2             2 ]
               [     - %1          %1      cos(theta)  - sin(theta)  ]

%1 :=                       cos(theta) sin(theta)

> Qbar.k:=evalm(inverse(T) &* Q);
       Qbark :=                2                         2
                 20. cos(theta)  + .7000000000 sin(theta)
                [-----------------------------------------,
                                     %2

                                          2                2
                    .7000000000 cos(theta)  + 2. sin(theta)
                    ----------------------------------------,  - 1.4 %4]
                                       %2

                                2                         2
                  20. sin(theta)  + .7000000000 cos(theta)
                 [-----------------------------------------,
                                      %2

                                           2                2
                     .7000000000 sin(theta)  + 2. cos(theta)
                     ----------------------------------------, 1.4 %4]
                                        %2

                                                     2             2
                                           cos(theta)  - sin(theta)
                   [19.3 %4,  - 1.3 %4, .7 -------------------------]
                                                       %3

                                     2           2
%1 :=                      cos(theta)  sin(theta)

                                4                     4
%2 :=                 cos(theta)  + 2. %1 + sin(theta)

                                4                    4
%3 :=                 cos(theta)  + 2 %1 + sin(theta)

                            cos(theta) sin(theta)
%4 :=                       ---------------------
                                      %3

> quit
bytes used=517204, alloc=458668, time=1.78

Kate Atherley
Technical Support
Waterloo Maple Software
support@maplesoft.on.ca
Newsgroups: sci.math.symbolic