deg2rad=pi/180;
x0=[180 0]*deg2rad;
global um r ym s net en1 fn1
en1=0;
fn1=0;
