I am getting an empty plot.

x1 = 0.01;
x2 = 0.2;
d0min = fminbnd(@Sold0, x1, x2);

d0min
disp(d0min);

d0 = 0.01:0.2;
diff = Sold0(d0);
plot (d0,diff)
xlim([0.01 0.08])
ylim([0 inf])
%xlabel('Tube outside diameter (d0)');
%ylabel('Assumed Flux(Q) - Obtained Flux(qr)');

>> Probably because, the y axis range is too big. Y varies from 10^6 till 10^14. How do I plot it here? PS: Attached files for reference.

Answer 1

[KSSV]

Replace plot command:

plot (d0,diffD)

with

plot (d0,diffD,'.')

[Devdatt Thengdi]

Okay. So, same code I should get Sold0(d0) over the entire range of d0 (0.01:0.2) but, I am getting only a single point at d0 = 0.01 and Sold0 = 3.132e^6. Why?

[KSSV]

Take d0 as a vector:

d0 = linspace(0.01,0.2);