#ifndef BANCROFT_H
#define BANCROFT_H

#include <newmatap.h> 

void bancroft(const Matrix& BBpass, ColumnVector& pos);

inline double lorentz(const ColumnVector& aa, const ColumnVector& bb);

#endif