B=.06//flux density
I=40D-3//current in coil
l=4D-2//length of coil side
F=B*I*l
N=50//no. of turns
mprintf("Force acting on each coil side=%f*10^-3 N",F*N*10^3)