how can I plot this function?

My problem is how to plot this function

[KALYAN ACHARJYA]

Pls refer here, helpful link

https://in.mathworks.com/help/symbolic/sym.symsum.html

[Niloufar]

The power in denominator is two

[David Hill]

Is it a function in f(x,r,t) or f(x,r) or f(r,t) or f(x,t) or f(x), f(r ) or f(t). What are the bounds or inputs for x, r, and t?

[Niloufar]

f(t) the bound for t is -5 to 5

[Niloufar]

close all;clear;clc; t = linspace(-5,5,1000); %ylim([-0.01 0.01]); N = 1000;

syms n t; an = (100*n*pi*(-1).^(n+1).*(1000-121*n.*n*pi*pi))/((1000-121*n.*n*pi*pi)^2 + (1320*n*pi)^2)*(50*n*pi); bn = (100*n*pi*(-1)^(n+1)*(1320*n*pi))/((1000-121*n.*n*pi*pi)^2 + (1320*n*pi)^2)*(50*n*pi); f = symsum(an.*cos(100*n*pi*t) + bn.*sin(100*n*pi*t),n,1,Inf); fplot(f);

I have done it with Symsum but it doesn't work.what is the problem

Answer 1

[David Hill]

t = linspace(-5,5,1000); 
N = 1000;n=1:N;
an = (100*n*pi.*(-1).^(n+1).*(1000-121*n.^2*pi^2))/(((1000-121*n.^2*pi^2).^2 + (1320*n*pi).^2).*(50*n*pi));
bn = (100*n*pi.*(-1).^(n+1).*(1320*n*pi))./(((1000-121*n.^2*pi^2).^2 + (1320*n*pi).^2).*(50*n*pi));
f = @(t)sum(an.*cos(100*n*pi*t) + bn.*sin(100*n*pi*t)); 
for k=1:length(t)
    F(k)=f(t(k));
end
plot(t,F)