Convert a for loop from MathCAD to Matlab

I have to solve for sigma c in the attached pdf and need some help converting it from Mathcad to Matlab. Can someone help me translate the code to Matlab? Here's the sizes of m,n,j:

%m = 1:71; 
%n = 1:41;
%j = 1:512;

Thanks!

[dpb]

And your attempt is where and you got stuck on what part, precisely???

[dpb]

What is the pseudo-code for the a conditional term?; I can't decipher the Matchcad-ese, sorry. Specifically, if I presume it's something like

a=dp/6;
a(0)=2*p(0)+p(1);
a(N-1)=2*p*(N-1)+p(N-2)*p(n-1)+4*p(n)+p(n+1);

then since the if implies n==N, then p(n+1) would seem to have a bounds problem???

And, of course, what are p, delta-p, k, etc., etc., ...

Answer 1

[Walter Roberson]

for m ε 0..M - 1

translates to the MATLAB code

for m = 0 : M - 1

However, keep in mind that in MATLAB, subscripts start at 1 instead of 0, so anywhere you wanted to use m as a subscript you would need to add 1 to it. Mostly you would end up writing as

for m = 1 : M

[Sean Raffetto]

I think the problem is not the code translation itself. There's more work ahead of this snippet that's not functioning correctly.

[Walter Roberson]

What is the question then?

[dpb]

"... work ahead of this snippet that's not functioning correctly."

And we were to be able to divine that how, precisely???