Create sigmafun.m