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

Go to the documentation of this file.

function y=myFp(x,n)

global unk


h=1/(n-1.0);


y=zeros(n*n*unk,1);

for j=2:n-1
    
    for i=2:n-1
        
        index=((j-1)*n+i-1)*unk+1;

        sc=x(index);
        sw=x(index-unk);
        se=x(index+unk);
        sn=x(index+n*unk);
        ss=x(index-n*unk);
        vc=x(index+1);
        vw=x(index+1-unk);
        ve=x(index+1+unk);
        vn=x(index+1+n*unk);
        vs=x(index+1-n*unk);
        
        dSTFdX=(se-sw)/(2.0*h);

                dSTFdY=(sn-ss)/(2.0*h);
        
                myU=dSTFdY;
                myV=-dSTFdX;
        
        dVORdX=(ve-vw)/(2.0*h);

                dVORdY=(vn-vs)/(2.0*h);
        
        
        UdVORdX=myU*dVORdX;
                
        VdVORdY=myV*dVORdY;
        
        
        y(index)=0;
        y(index+1)=UdVORdX+VdVORdY;
 
    end
   
end

---