parameter[e_tot] = 15e6
parameter[geometry] = closed
parameter[particle] = electron

ldl = 0.20
lds = 0.05
lqf = 0.10
lqd = 0.10
ang = pi/21.0
g = 3.5

tdl = 2.0*atan(ldl*sin(0.5*ang)/(lds+ldl*cos(0.5*ang)))
tds = 2.0*atan(lds*sin(0.5*ang)/(ldl+lds*cos(0.5*ang)))
r = 0.5*ldl/sin(0.5*tdl)

dlh: sbend, l=0.5*r*tdl, g=1.0/r, g_err=-1.0/r
ds: sbend, l=r*tds, g=1.0/r, g_err=-1.0/r
qf: sbend, field_master=t, l=lqf, 
	b_field=0.0, 
	b1_gradient=-g
qd: sbend, field_master=t, l=lqd, b_field=0.0, b1_gradient=g,
	b_field_err=-2.0*sin(0.5*ang)*beginning[p0c]/(c_light*lqd) - g/cos(0.5*ang)*(0.5*lqf*sin(0.5*(tdl-tds))+0.5*lqd*sin(0.5*ang)+r*sin(0.5*tdl)*sin(0.5*tds))

cell: line=(dlh,qf,ds,qd,dlh)