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(*,*)thd,sin(thr),disc
thd=thd+0.1
end do
end program sinetheta