T0 =[383.49500 383.36684 383.20479 383.00086 382.74576 382.42899 382.03927 381.56517 380.99632 380.32493 379.54776 378.66813 377.69746 376.65586 375.57116 374.47628 373.40536 372.38966 371.45409 370.61525 369.88114 369.55450 369.15979 368.68736 368.12836 367.47581 366.72599 365.87989 364.94444 363.93329 362.86662 361.76992 360.67171 359.60070 358.58279 357.63881 356.78307 356.02320 355.36073 354.79237 354.31141 ]';
T1 =[383.16783 382.90758 382.58409 382.18577 381.70096 381.11914 380.43259 379.63843 378.74067 377.75177 376.69315 375.59395 374.48821 373.41077 372.39301 371.45944 370.62588 369.89936 369.27931 368.75957 368.33044 367.73681 367.05478 366.28404 365.42916 364.50038 363.51372 362.49029 361.45463 360.43249 359.44829 358.52292 357.67210 356.90577 356.22826 355.63902 355.13384 354.70601 354.34744 354.04954 353.80382 ]';
T2 =[383.60000 383.51509 383.40675 383.26894 383.09435 382.87428 382.59865 382.25618 381.83488 381.32289 380.70982 379.98862 379.15767 378.22291 377.19927 376.11076 374.98886 373.86917 372.78700 371.77306 370.85032 370.67096 370.44694 370.16866 369.82532 369.40520 368.89636 368.28760 367.56992 366.73838 365.79403 364.74568 363.61084 362.41541 361.19177 359.97549 358.80132 357.69929 356.69200 355.79336 355.00880 ]';
stage =[1.00000 2.00000 3.00000 4.00000 5.00000 6.00000 7.00000 8.00000 9.00000 10.00000 11.00000 12.00000 13.00000 14.00000 15.00000 16.00000 17.00000 18.00000 19.00000 20.00000 21.00000 22.00000 23.00000 24.00000 25.00000 26.00000 27.00000 28.00000 29.00000 30.00000 31.00000 32.00000 33.00000 34.00000 35.00000 36.00000 37.00000 38.00000 39.00000 40.00000 41.00000 ]';
hold on
xlabel('Stage number')
ylabel('Temperature')
axis([0 41 354.82 383.56])
plot(stage,T0,'-k','Linewidth',1.1)
plot(stage,T1,'--b','Linewidth',1)
plot(stage,T2,'--r','Linewidth',1)
legend('Nominal','\Delta LT = +0.01','\Delta VB = +0.01')
set(gca,'box','on','ticklength',[0.02 0.05])
hold off