Create alp_beta.m

function that finds the shape parameter beta and variance parameter alpha using ML estimate.