MATLAB: How to calculate curl of Laguerre-Gaussian modes

curlgaussian beam

curl(LG mode)

Refer below link:

The answer defined in this way:

syms Fx(x,y,z) Fy(x,y,z) Fz(x,y,z)
curl([Fx Fy,Fz],[x, y, z])

Result:

ans(x, y, z) =
 diff(Fz(x, y, z), y) - diff(Fy(x, y, z), z)
 diff(Fx(x, y, z), z) - diff(Fz(x, y, z), x)
 diff(Fy(x, y, z), x) - diff(Fx(x, y, z), y)

Write the first line as:

function [y,z] = myfunc()

If you only request one output, only the first output is returned by default.