Create calculateAB.m

Function to calculate a and b for the stiffness matrix