Sure, it has to be possible.
But you cannot expect a routine which is written for a scalar T to work for a vector T, as you did.
The easies way in my opinion would be to write an additional loop inside your program that cycles through all values of T. You would have to put the caclculations of all variables which depend on T inside this program loop.
That way you could return a 1560x1000 matrix, but if you are anly after the maximum I would rewrite the program to just remember the max after each iteration and return the max only.
Is there a reason you use square root of uuu^2 rather than the absolute value of uuu?
