You can’t.
Use meshgrid or ndgrid to calculate the ‘x’ and ‘t’ matrices. (You will have to have them as matrices for your function to work.)
You don’t need to do an infinite sum. Instead, use a while loop to stop when the value of ‘T(x,t)’ is no longer changing between iterations.
Example (sine series):
x = pi/9;
n = 0;
S = 0;
Sp = Inf;
err = 1;
while err > 1E-8
S = S + ((-1)^n)/factorial(2*n+1) * x^(2*n+1)
err = abs(S - Sp);
Sp = S;
n = n + 1;
end
That produces a decent approximation, and converges in a few iterations (depending on the threshold set for ‘err’). Since your function will return a matrix, you will have to write your loop to recognise that, and calculate ‘err’ (in this example, call it what you like in your code) correctly.
