User-defined function help

function [Temp] = TempProfile(r)
%UNTITLED9 Summary of this function goes here
%   Detailed explanation goes here
Temp = 1500-(((9.61*10e26)*r.^2)/(4*0.05));
end

I have the function above. I am trying to display the Temp for the value of r; r = 0:0.41 How would I display all the values of Temp for all the values of r? Where would I define r as 0:0.41? Thank you so much!

[Walter Roberson]

There are approximately 2^59 representable numbers in the range 0 to 0.41. You do not have have enough memory to store all of the calculations and calculating them all would take many centuries.

I recommend that you reconsider displaying all of those values. Perhaps you do not need more than (say) 1000000000000 of the values in that range.

[Prajith Chilummula]

Consider the number of values you want to use between 0 and 0.41 and modify the increment step likely.Then you can pass the vector as argument and get the output temp array.You can plot the arrays.

Example:

r=0:0.1:0.41;
temp=TempProfile(r);
plot(r,temp);