volumenes parcialas atravez de matlap

Para la determinación de los volúmenes molares parciales se uso el método de la derivada
tangente con Matlab en donde tenia la siguiente programación.
clc
fprintf('CALCULO DE LOS VOLUMENES MOLARESPARCIALES DE UNA MEZCLA BINARIA
A\n')
fprintf('PARTIR DE DATOS EXPERIMENTALES.\n')
fprintf('\n')
fprintf('METODO DE LA PENDIENTE DE LA TANGENTE.\n')
fprintf('SISTEMA POR DEFAULT ESAGUA(1)/METANOL(2)\n')
fprintf('V2 es el vol del más volátil (metanol) y V1 el menos volátil
(agua)\n')
fprintf('W1 es la masa del picnometro vacío y seco y W2 la del picnometro
con la mezcla\n')
fprintf('y Vp el vol delpic y dsol son las densidades de las
soluciones\n')
fprintf('y n es el grado del polinomio\n')
datos = inputdlg({'volmasvolátil (V2) en ml','volmenosvolátil (V1)','masa
pic (W1) en gr','masapic+mezcla (W2)','vol pic (Vp)','masa molar met
(Mm)','Grado del polinomio de ajuste'},'Ingrese los
datos',1,{'[]','[]','','[]','','',''});
V2 = str2num(datos{1});
V1 = str2num(datos{2});
W1 =str2num(datos{3});
W2 = str2num(datos{4});
Vp = str2num(datos{5});
Mm = str2num(datos{6});
n = str2num(datos{7});
fprintf('La densidad de la soluciones dsol=(W2-W1)/Vp son:')
dsol=(W2-W1)./Vp
fprintf('Observe quedsol(9) es la densidad del más volátil\n')
fprintf('y que dsol(1) es la densidad del menos volátil\n')
md=((V2.*dsol(9))./(V1.*dsol(1)))*(1000/Mm);
%Donde Mm es la masa molar del metanol
%Calculando masasol Wsol con 1 kg de solvente
Wsol=md.*Mm+1000;
% Vs es el volumen de la solucion con 1 kg de solvente
Vs=Wsol./dsol;
Vms=Vs./(55.5556+md);
%Calculo de fraccion molar
Xm=md./(md+55.5556);
%EliminandoNaN
VsExp=[Vs(2:8)];
fprintf('LA MOLALIDAD DE LAS SOLUCIONES CON RESPECTO AL MÁS VOLÁTIL
SON:\n')
md=[md(2:8)]
P=polyfit(md,VsExp,n);
VsCalc=polyval(P,md);
fprintf('OBSERVE EN LA GRAFICA SI ELPOLINOMIO SE AJUSTA A LOS DATOS
EXPERIMENTALES.\n')
fprintf('EN CASO CONTRARIO CAMBIE EL GRADO DEL POLINOMIO.\n')
%Calculando la derivadas de VsCalc
fprintf('\n')
fprintf('ESTOS SON LOS VOLUMENES MOLARES...