Problem running symmgrp2009.max with recent versions of Maxima



Am Dienstag, den 15.11.2011, 12:51 +0000 schrieb Wiltshire R J (AT):
> I do hope the files I have sent you are helpful. Please let me know if they are not and I shall see if I can provide alternatives.
> 
> Ron Wiltshire
> 
> On 13 Nov 2011, at 15:05,  wrote:

I have extracted out of the informations the following test example. The
files are renamed to fit the standard extensions of Maxima.

load("symmgrp2009.mac")$
load("h_hd.mac")$

symmetry(1,0,0)$
printeqn(lode)$

The following shows a complete session in a Maxima console. I think it
works correctly. There seems to be no problem.

(%i1) load("symmgrp2009.mac")$

Code symmgrp2009.max of Janurary 1, 2010 is being loaded. 
Code symmgrp2009.max of January 1, 2010 was successfully loaded. 
(%i2) load("h_hd.mac")$

(%i3) symmetry(1,0,0)$

  
Code runs under the commercial software Macsyma as well as the 
public domain software Maxima (version 3.16.3 or higher). 
Please report problems to Willy Hereman, whereman at mines.edu. 
 -----------------------------------------------------------------  
Computation started. Execution may be slow for complicated cases! 
 -----------------------------------------------------------------  
/*********************************************************/ 
/*  Welcome to the Macsyma program for the computation   */ 
/*   of Lie-point symmetries of differential equations  */ 
/*     Written by Benoit Champagne and Willy Hereman     */ 
/*  Code adjusted for the public domain software Maxima  */ 
/*         by Benoit Huard and Willy Hereman             */ 
/*       Project supervision: Pavel Winternitz           */ 
/*        Version 3.1 released on January 1, 2010        */ 
/*           Program name: symmgrp2009.max               */ 
/*                Copyright 1991-2009                    */ 
/*********************************************************/ 
 -----------------------------------------------------------------  
Computation started. Execution may be slow for complicated cases! 
 -----------------------------------------------------------------  
*** Number of determining equations before simplifications: 21 . *** 
Number determining eqs. in lode[1] before removing duplicates: 7 
Number determining eqs. in lode[1] after removing duplicates: 6 
Number determining eqs. in lode[2] before removing duplicates: 14 
Number determining eqs. in lode[2] after removing duplicates: 13 
List of factors that were canceled:  [u[1]^2,u[1]^3,u[1]^6] 
Number of determining equations after simplifications: 8 
Number of determining equations before removing duplicates: 8 
Number of determining equations after removing duplicates: 8 
Number of determining equations before resorting lode: 8 
Number of determining equations after resorting lode: 8 
*** Number of determining equations after all simplifications: 8 . *** 
*** These determining equations are stored in lode. *** 
(%i4) printeqn(lode)$

  
Equation 1 : 'diff(eta[2],u[1],1) = 0 
Equation 2 : 'diff(eta[2],x[1],1) = 0 
Equation 3 : 'diff(eta[1],u[1],1) = 0 
Equation 4 : 'diff(phi[1],u[1],2) = 0 
Equation 5 : 
'diff(phi[1],u[1],1,x[1],1)-'diff(eta[1],x[1],2) = 0 
Equation 6 : 
'diff(phi[1],x[2],1)-u[1]^3*'diff(phi[1],x[1],3) = 0 
Equation 7 : 
3*u[1]^3*'diff(phi[1],u[1],1,x[1],2)+'diff(eta[1],x[2],1)
                                    -u[1]^3*'diff(eta[1],x[1],3)
  = 0
  
Equation 8 : 
u[1]*'diff(eta[2],x[2],1)-3*u[1]*'diff(eta[1],x[1],1)+3*phi[1] = 0 

Dieter Kaiser