You need to save the output of fun into specific elements of T, P, and D. E.g.,
h=0:1:105;
for k=1:numel(h)
[T(k),P(k),D(k)]=fun(h(k));
end
And account for an h input of 0 in your fun.m code:
if(0<=h && h<=11) % <-- changed 0<h to 0<=h
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!