How to fix this error

I need some help fixing this error, I assume it's in the way I call the function. Any help would be much obliged.

Answer 1

[James Tursa]

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