contrast.m control.emf diophantine.p getG.p getGij.p GPC_basic.m GPC_basic1.m GPC_basic2.m output.emf rankn.p tf2AB.p