I am new to Matlab and would like to solve the following equation:
β*∂T/∂t = q*dx/dt - K*∂T/∂z
I have discretized it in space to make it a system of odes and I have come up with:
β*(∂T_i)/∂t = q*dx/dt - K*1/2h*(T_(i+1)-T_(i-1)
How do I implement this in matlab to solve it?