require(gtools)

data<-as.data.frame(permutations(9,5,1:9))

data<-data[!(data$V1==4|data$V1==6|data$V1==8),]

data$LHS<-data$V1*as.numeric(paste(as.character(data$V2),as.character(data$V2),sep=""))*as.numeric(paste(as.character(data$V2),as.character(data$V3),sep=""))^2
data$RHS<-as.numeric(paste(as.character(data$V4),as.character(data$V5),sep=""))^2-data$V2
data$DIFF<-data$LHS-data$RHS
data<-data[data$DIFF==0,]

data