D:/demona/2.0/manual/doc/2-31set/matlab/fs.m

Go to the documentation of this file.

function y=fs(s,v,n)

y=zeros(n*n,1);

y=s;

h=1/(n-1.0);

for j=2:n-1
    
    for i=2:n-1
        
        index=(j-1)*n+i;
        
        sc=s(index);
        sw=s(index-1);
        se=s(index+1);
        sn=s(index+n);
        ss=s(index-n);
        vc=v(index);
                       
        y(index)=(ss+sn+sw+se-4*sc)/h^2+vc;
          
                
    end
    
    
end

---