! 求解二次方程 a*x**2+b*x+c=0
print *, '輸入3個系數'; read *, a,b,c

if(abs(a).le.tiny(a)) then
	print *, 'a的值小于允許的最小實數，認為是零'
	if(abs(b).le.tiny(b)) then
		print *, 'b的值小于允許的最小實數，認為是零'
		if(abs(c).le.tiny(c)) then
			print *, 'c的值小于允許的最小實數，認為是零'
			print *, '恒等式0=0，無需解'
		else
			print *, '無解'
		end if
	else
		print *, '一個解：x=', -c/b
	end if
else
	p1=-b/(2*a); s=b**2-4*a*c; p2=sqrt(abs(s))/(2*a)
	if(s<0.) then
		print *, '兩個復數解：x=', p1, ' +-', p2, 'i'
	else if (abs(s).le.tiny(s)) then
		print *, '重根      ：x=', p1
	else
		print *, '兩個實數解：x=', p1, ' +-', p2
	end if
end if
end