请教求和及画图

charlyluo

我的目的是画一个与g有关的求和表达式的图
程序如下：
clear
syms g
a=7;
n=100;
I1=0;
for i=1:n
f1=(cos(sqrt(2*i+1)*g))^2;
f1=a^(4*i+2)/gamma(2*i+2)*f1;
I1=I1+f1;
end
g=0:20;
plot(g,I1);

提示的错误：
??? Error using ==> plot
Conversion to double from sym is not possible.
希望高手解答

charlyluo

在画图之前，我可以求出I1的表达式，可以一旦我要去画图的时候，就有那个错误提醒

friendchj

参考（注意 n过大会出现NAN）：
clc
clear
syms g
a=7;
n=10;
I1=0;
for i=1:n
f1=(cos(sqrt(2*i+1)*g))^2;
f1=a^(4*i+2)/gamma(2*i+2)*f1;
I1=I1+f1;
end
x=0:20;
plot(x,subs(I1,g,x))