R: Magic squares prime order

magic.prime {magic} | R Documentation

## Magic squares prime order

### Description

Produces magic squares of prime order using the standard method

### Usage

    magic.prime(n,i=2,j=3)

### Arguments

`n` | The order of the square
`i` | row number of increment
`j` | column number of increment

### Details

Claimed to work for prime order, but I've tried it (with the defaults for `i` and `j`) for many composite integers of the form 6n+1 and 6n-1 and found no exceptions; indeed, they all seem to be panmagic. It is not clear to me when the process works and when it doesn't.

### Author(s)

Robin K. S. Hankin

### References

http://www.magic-squares.de/general/general.html

### Examples

    magic.prime(7)
    f <- function(n){is.magic(magic.prime(n))}
    all(sapply(6*1:30+1,f))
    all(sapply(6*1:30-1,f))

    is.magic(magic.prime(9,i=2,j=4),give.answers=TRUE)
    magic.prime(7,i=2,j=4)

[Package magic version 1.5-6 Index]