Symbolic Math Toolbox may be used to implement the above equation in MATLAB. Here are a few suggestions which may be helpful in the implementation
- Declare “n” as a symbolic variable and “N” as a symbolic function of two variables “x” and “t” using the following syntax
- To implement the summation term, use the following code
F = exp(-((2*n+1)*pi)^2*D*t/L^2)*sin((2*n+1)*pi*x/L)/((2*n+1)*pi^3);
Fsum = symsum(F,n,0,inf);
Here “n” is replaced by “2n+1” so that only odd terms are considered in the summation.
The function “symsum()” sets the limits of summation from 0 to inf.
Remaining terms can be written as normal equation in MATLAB
Please use the following documentation links on how to create symbolic functions and “symsum()” to get a better understanding