I understand that you want to write a script that automate the formula given below:
Psi_cx(m) = psi_0 + M_cvec(n,m)*cos(n*(2*pi*x/a))*pi^(-1/2)
You can do it using following script:
M_cvec = magic(9);
psi_0 = 2*pi^(-1/2);
%Create symbolic scalar variables
syms x;
syms a;
%lets m=1 and n=9
m=1;
n=9;
%storing result of each 'm' values.
Psi_cx=[];
for i=1:m
Psi_cxTemp=psi_0;
for j=1:n
Psi_cxTemp=Psi_cxTemp+M_cvec(j,i)*cos(j*(2*pi*x/a ))*pi^(-1/2);
end
Psi_cx=[Psi_cx Psi_cxTemp];
end
%Now use "subs" for substituting symbolic scalar variables with their respective values.
% Example : for x=1,m=1 and a=1
Psicx1=subs(Psi_cx(1),x,1);
Psicx1=subs(Psicx1,a,1);
Here, we are using "x" and "a" as a symbolic scalar variable so that we can substitute it with their respective values later using "subs" function.
Please refer to the following documentation for more information on "Symbolic scalar variables" and "Substituting symbolic scalar variables":
- For "Symbolic scalar variables": Create symbolic scalar variables and functions, and matrix variables and functions - MATLAB syms - MathWorks India
- For "Substituting symbolic scalar variables": Symbolic substitution - MATLAB subs - MathWorks India
