|
发表于 2023-5-25 08:13:10
|
显示全部楼层
for i=-0.4:0.1:0.4
fimplicit(@(x,y)exp(a*a*t)*y+A00*y./(x.^2+y.^2)+B00*y.*besselk(1,a*sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2)+Re*c*x.*y./(x.^2+y.^2).*A11./(x.^2+y.^2)+B11.*besselk(2,sqrt(2)*a*sqrt(x.^2+y.^2))+B00*besselk(2,a*sqrt(x.^2+y.^2))/a/2))-i,[-15,15,-15,15]);
end |
|