dta.m

来自「C++编译指导」· M 代码 · 共 4 行

function u=dta(a,b,c,eps)
   N=size(a,1);
   u=(a.^2+b.^2).*c+(36*(eps^2)).*ones(N,1);
   u=(c./(sqrt(u)))/2;