Matlab - cómo reemplazar una variable por otra ?

%Vamos a llamar a la matriz de coeficiente estequiométrico como A%

A=[1,0,-2,0,0; 2, 0, 0, -1,0; 0 , 2 , -3, -2,0]

A(2,:)= -2*A(1,:)+A(2,:)

A([2 3],:)=A([3 2],:)

%Armamos las ecuaciones%

syms a b c d t

Ec1= 0==a.*A(1,1) + b.*A(1,2)+c.*A(1,3)+d.*A(1,4)

Ec2= 0==a.*A(2,1) + b.*A(2,2)+c.*A(2,3)+d.*A(2,4)

Ec3= 0==a.*A(3,1) + b.*A(3,2)+c.*A(3,3)+d.*A(3,4)

c= solve(Ec3,c)

b= solve(Ec2,b)

a= solve(Ec1,a)

%entonces, obtengo las variables a,b, y c; pero me gustaría dejar todo en funcion de la variable d, es decir, tengo que a=2c y que c=d/4. Como puedo hacer que a sea igual a d/2

A=[1,0,-2,0,0; 2, 0, 0, -1,0; 0 , 2 , -3, -2,0];

A(2,:)= -2*A(1,:)+A(2,:)

A([2 3],:)=A([3 2],:)

%Armamos las ecuaciones%

syms a b c d t

Ec1= 0==a.*A(1,1) + b.*A(1,2)+c.*A(1,3)+d.*A(1,4);

Ec2= 0==a.*A(2,1) + b.*A(2,2)+c.*A(2,3)+d.*A(2,4);

Ec3= 0==a.*A(3,1) + b.*A(3,2)+c.*A(3,3)+d.*A(3,4);

 [valorc,valorb,valora]=solve([Ec1; Ec2; Ec2],[c; b; a])