program sinetheta
real::thd,thr,disc,tol
! tolerance in percentage
tol=0.1
thd=0.1
do while(disc<tol)
thr=thd*3.1416/180
disc=((thr-sin(thr))/thr)*100
write(*,100)thd,sin(thr),disc
thd=thd+0.1
end do
100 format(f8.2,8x,f8.6,8x,f8.6)
end program sinetheta