# Python stdlib from math import sqrt, pi from random import randint # Numpy stuff from Numeric import * from MLab import eye from LinearAlgebra import singular_value_decomposition, determinant from RandomArray import random, normal # Biopython stuff from Bio.PDB import * S=eye(3,3, typecode='d') S[2,2]=-1 print S**3.8