tauD=((MvTotD*16)/(pi*slutDiameterD.^3)+ (4*TmaxD/(pi*slutDiameterD.^2)));

h3=figure('name','Tau-plots');
plot(Xiy, tauD)

sigmaD = ((4*Nmax)/(pi*D.^2) + (32*MtotD)/(pi*D.^3))

h4=figure('name','sigma-plot');
plot(Xiy, sigmaD)

mises=sqrt(3*(tauD.^2) + sigmaD.^2)
h5=figure('name','mises');
plot(mises)