regbndpointsro.m

本代码为基于matlab的求解偏微分方程的利用径向基函数的无网格方法的源码

function bndpoints=regbndpointsRO(h) % defines "regular" boundary values 
                           % with approx. stepsize h on [0,2*pi]
global RBFdomain
 n=ceil(pi/h);  % the stepsize must be readjusted 
 h=2*pi/n;        % to make the angles lie at equidistant values
 phi=0:h:(2*pi) ;
 phi(length(phi))=[];% careful: no duplication of last=first point
 [nx,nphi]=size(phi);
 z=zeros(nphi,2);
 r=bndfct(phi);
 [x,y]=pol2cart(phi,r);
 bndpoints=[x' y'];