无网格方法解悬臂梁问题

无网格方法解悬臂梁问题，自己写的，收个辛苦分

function [wdwdxdwdy]=cubwgt（diftvdmaxdm）
l=length（v）;
for i=1:l
    drdx=sign（dif（1i））/dm（1v（i））;
    drdy=sign（dif（2i））/dm（2v（i））;
    rx=abs（dif（1i））/dm（1v（i））;
    ry=abs（dif（2i））/dm（2v（i））;
    if rx>0.5
        wx=（4/3）-4*rx+4*rx*rx-（4/3）*rx^3;
        dwx=（-4+8*rx-4*rx^2）*drdx;
    elseif rx<=0.5
        wx=（2/3）-4*rx*rx+4*rx^3;
        dwx=（-8*rx+12*rx^2）*drdx;
    end
    if ry>0.5
        wy=（4/3）-4*ry+4*ry*ry-（4/3）*ry^3;
        dwy=（-4+8*ry-4*ry^2）*drdy;
    elseif ry<=0.5
        wy=（2/3）-4*ry*ry+4*ry^3;
        dwy=（-8*ry+12*ry^2）*drdy;
    end
    w（i）=wx*wy;
    dwdx（i）=wy*dwx;
    dwdy（i）=wx*dwy;
end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-04-28 19:35  无网格方法解悬臂梁问题\
     文件         666  2014-04-28 19:35  无网格方法解悬臂梁问题\cubwgt.m
     文件         141  2014-04-28 19:30  无网格方法解悬臂梁问题\domain.m
     文件         695  2014-04-28 19:34  无网格方法解悬臂梁问题\egauss.m
     文件        4590  2014-04-28 19:29  无网格方法解悬臂梁问题\ff.m
     文件         221  2014-04-28 19:30  无网格方法解悬臂梁问题\gauss2.m
     文件        1252  2014-04-28 19:32  无网格方法解悬臂梁问题\mesh2.m
     文件         957  2014-04-28 19:31  无网格方法解悬臂梁问题\shape.m