fun3 = fun1 – fun 2 , where;
fun1= @(x1) (i*w*x1./U)*(-1./r.^3)+(3*(x1*x1)./(r.^5));
fun2=fun2= @(x1) exp(i*w*x1./U).*(((k2.^2)*x1*x1./r.^3)*exp(-k2*r+k*x1)-(k2*k*x1*1./(r.^2))*exp(-k2*r+k*x1)-(k2*1./(r.^2))*exp(-k2*r+k*x1)…
+(2*k2*x1*x1./(r.^4))*exp(-k2*r+k*x1)-(k2*k*x1*1./(r.^2))*exp(-k2*r+k*x1)+((k.^2)*1./r)*exp(-k2*r+k*x1)…
-(k*1*x1./(r.^3))*exp(-k2*r+k*x1)-(k2*x1*x1./(r.^4))*exp(-k2*r+k*x1)+(k*x1*1./(r.^3))*exp(-k2*r+k*x1)…
+(1./(r.^3))*exp(-k2*r+k*x1)-(3*x1*x1./(r.^5))*exp(-k2*r+k*x1));
where: x2=0; ru=1; mu=1; w=0.01; U=1;
r = sqrt((x1).^2 + (x2).^2 + (x3).^2); k2 =sqrt((k.^2)+(h.^2)); h =sqrt(-i*ru*w./mu);
Best Answer