Problem running symmgrp2009.max with recent versions of Maxima
Subject: Problem running symmgrp2009.max with recent versions of Maxima
From: Dieter Kaiser
Date: Wed, 16 Nov 2011 19:47:07 +0100
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