Biot-Savart direct integration on a generic curve

Biot-Savart integration on a generic curve, by numerical integration.
You can insert the points of a generic curve and the script will calculate the B induction field everywhere.

You can choose between three test cases: simple straight wire, bent wire or solenoid, or you can insert the coordinates of your generic curve.
The script will also plot the quiver3 of the field.

Image in the screenshot has been simply taken using streamtube() function:

[xm, ym] = meshgrid(linspace(Dom(1,1), Dom(1,2), ND), linspace(Dom(2,1), Dom(2,2), ND));
zm = xm * 0;

streamtube(stream3(X,Y,Z,U,V,W,xm,ym,zm+6), 0.1);

You can find more info (in Italian) here:
http://www.wirgilio.it/blog/2013/06/19/campo-magnetico-e-velocita-indotta-da-un-vortice-lungo-una-curva-generica-con-matlab/