add mathematica code, try new c/d vals