[問題] 請問function迴圈的問題已回收
副程式
function[w,dwdx,dwdy,dwdxx,dwdxdy,dwdyy]=plot1(type,para,fx,fy,dm)
%weight('SPLIN',2,2,2,1)
%type-有二種指數(Gauss)和多項式(Spline)
%para-常數如c
%w,dwdx,dwdy...表示微分
%--------fx=(x-Nx)
%--------fy=(y-Ny)
f=sqrt(fx^2+fy^2);
if (fx>=0)
dfdx=(2/f)*fx;
else
dfdx=(-2/f)*fx;
end
if (fy>=0)
dfdy=(2/f)*fy;
else
dfdy=(-2/f)*fy;
end
-------------------------------------------主程式的一部分
for a = 1 : npoints %高斯節點
for b = 1 : npoints
for i = 1 : nnodes %自已畫的節點
for j = 1 : nnodes
fx=xx(a)-xxi(i);
fy=yy(b)-yyi(j);
f=sqrt(fx^2+fy^2);
[wn(i),dwndx(i),dwndy(i),dwndxx(i),dwndxdy(i),dwndyy(i)]=plot1('GAUSS',4.2,fx,fy,dm);
end
end
但是他不會帶主程式的每個fx和fy進去耶
請問要怎麼修改
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.134.88