Recibí unafor i=1:15. Dentro genero una variabled=1:0.01:10, que es el eje x'x y en base a esto, creo una función continua F (d) que tiene 2 variables únicas pitch y yaw. Luego trazo esto usando diferentes colores en cada recursión usandocmap = hsv(15);. Entonces entonces es:

d=1:0.01:10;
cmap = hsv(15);

for i=1:15
    pitch = unidrnd(10);
    yaw   = unidrnd(10);

    for j=1:length(d)
        F(j) = d(j)*3*pitch*yaw; %// some long calculation here
    end

    p1 = plot(d,F,'Linewidth', 1.0);
    title ('blah blah')
    set(p1, 'Color', cmap(i,:));
    hold on;
    legend (['pitch,yaw:', num2str(pitch) num2str(yaw)]) 
end 
hold off;

Este código actualiza los valores únicos de inclinación, guiñada en cada recursión (sin espacio entre ellos, por lo que es un poco irritante) pero no logra:

Aplique el color apropiado, visible en la figura.Mantenga el color de la iteración anterior y los valores depitch,yaw.