dw.m

用无网格法计算VAN DE POL 振子

clear all,clc
syms r b
gau = (exp(-r^2*b^2)-exp(-b^2))/(1-exp(-b^2));
gau1 = simplify(diff(gau,r,3))
gau2 = simplify(diff(gau,r,4))


cub = 1-6*r^2+8*r^3-3*r^4;
cub1 = simplify(diff(cub,r,3))
cub2 = simplify(diff(cub,r,4))


spla =  2/3 - 4*r^2 + 4*r^3;
spla1 = simplify(diff(spla,r,3))
spla2 = simplify(diff(spla,r,4))

splb = 4/3 - 4*r + 4*r^2 - 4*r^3/3;
splb1 = simplify(diff(splb,r,3))
splb2 = simplify(diff(splb,r,4))


CSR = (1-r)^6*(6+36*r+82*r^2+72*r^3+30*r^4+5*r^5);
CSR1 = simplify(diff(CSR,r,3))
CSR2 = simplify(diff(CSR,r,4))