syms a b c d
[a,b,c,d]=solve('0.005083-0.2985*a-0.00025*b+0.0725*c+0.0875*d=0','0.009083+0.00025*a+0.306*b+0.2025*c-0.14*d=0',...
'0.064167+0.0725*a-0.2025*b-0.10575*c-0.0225*d=0','0.006667+0.0875*a+0.14*b-0.0225*c+0.23325*d=0')
a = -1.264915611178603343104260187311
b =2.0229725158204830541709095591987
c =-3.8908749382444142140830410100597
d =-1.1436129575444336757949538890652