program sinetheta
real :: thd,thr,disc,tol
!tol
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
write(*,*)"discrepency>tolerence for angle=",thd,"degree"
end program sinetheta