Hi Clayton Holmes,
Here are some assumptions i have made for writing the code :-
- All the variables and the elements in matices used in the summation have been initialised to 1.
- value of infinity is assigned with 10 ,which can be changed accordingly.
Below i have added the code for the given summation.
% initialising constant
inf =10;
x=1;
y=1;
z=1;
lambda=1;
t=1;
a=1;
B=ones([1,inf+1]);
N=ones([1,inf+1]);
U=ones([1,inf+1]);
C=ones([inf+1,inf+1,inf+1]);
sum=0;
for n=1:inf+1
for m=1:inf+1
for p=1:inf
sum=sum+cos(B(n)*z)*sin(N(m)*y)*sin(U(p)*y)*exp(-(lambda*lambda)*a*t)*C(n,m,p);
end
end
end
%value of sum
sum
