Matlab cálculo simbólico, sustituir datos

[draykaler]

Buenos días, estoy intentando sacar la incertidumbre del campo magnético mediante matlab. He conseguido hacerlo bien con un copia pega de la ventana de comandos al script de las derivadas, pero intento hacerlo todo en un script.
Lo que intento, y pregunto, es como hacer para que después de hacer el cálculo de las derivadas parciales las variables simbólicas se sustituyan por su valor correspondiente.
Lo que he intentado pero no funciona es:

syms u i r d

A =diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),i);
B=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),r);
C=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),d);

i=double(50)
Di=2;
r=double(0.03);
Dr=0.001;
d=double(0.01);
Dd=0.001;
u=double(4*pi*10^(-7));

incert=sqrt((A*Di)^2+(B*Dr)^2+(C*Dd)^2)

Muchas gracias

[draykaler]

Re: cálculo simbólico, sustituir datos

Ya encontre la solución, por si a alguien le interesa:

%%%% CALCULO DE INCERTIDUMBRES DE FORMA SENCILLA EN MATLAB
%%%% ejemplo hecho con el campo magnético en el eje de la espira


syms u i r d

A=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),i);
B=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),r);
C=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),d);


Di=2; %incertidumbres sistematicas
Dr=0.001;
Dd=0.001;

%%%%%%%% I r d u
X=subs(A,{i,r,d,u},{50,0.03,0.02,4*pi*10^(-7)});%cuidado con los paréntesis y los corchetes
Y=subs(B,{i,r,d,u},{50,0.03,0.02,4*pi*10^(-7)});
Z=subs(C,{i,r,d,u},{50,0.03,0.02,4*pi*10^(-7)});

incert=sqrt((X*Di)^2+(Y*Dr)^2+(Z*Dd)^2)