I have a solution obtained by variational form of RT0 fespace

```
fespace Eh(Th,RT0Ortho);
Eh<complex> [Ax,Ay];
varf vK11([ux,uy],[vx,vy])
=int2d(Th)( nu * (dx(vy)-dy(vx)) * (dx(uy)-dy(ux))
+ 1i*2*pi*freqf*sig* [vx,vy]' * [ux,uy] )
-int2d(Th)( sig * [vx,vy]' * [E0x,E0y] )
+on(CoilTerminalIn,CoilTerminalOut, ux=0,uy=0)
+on(AirInBot1,AirOutBot,AirOutLeft, ux=0,uy=0);
matrix<complex> Kside = vK11(Eh,Eh);
complex[int] bside = vK11(0,Eh);
set(Kside, solver=UMFPACK);
complex[int] sol;
sol = Kside^-1 * bside;
```

how can I dispatch `Ax`

and `Ay`

from `sol`

it seems `[Ax[],Ay[]]=sol;`

is the wrong approach.