Исходник программы - определение точек пересечения прямой с коружностью на Си(C++)

Код на C++
#include<conio.h>
#include<stdio.h>
#include<iostream.h>
#include<math.h>
 
void prover(float x,float y,float k,float b,float r)
{ //находим дискременант квадратного уравнения
      float d=(pow((2*k*b-2*x-2*y*k),2)-(4+4*k*k)*(b*b-r*r+x*x+y*y-2*y*b));
 
  //если он равен 0, уравнение не имеет решения
     if(d<0) {cout<<"\n\nПрямая и окружность не пересекаются";return;}
  
  //иначе находим корни квадратного уравнения
     float x1=((-(2*k*b-2*x-2*y*k)-sqrt(d))/(2+2*k*k));
     float x2=((-(2*k*b-2*x-2*y*k)+sqrt(d))/(2+2*k*k));
 //если абсциссы точек совпадают, то пересечение только в одной точке    
     if (x1==x2)  {printf("\n\nПрямая и окружность имеют точку касания:(%.2f,%.2f)",x1,y1); 
     return;}
 //иначе находим ординаты точек пересечения 
     float y1=k*x1+b;
     float y2=k*x2+b;
     printf("\n\nПрямая и окружность пересекаются в точках: (%.2f,%.2f), 
     (%.2f,%.2f)",x1,y1,x2,y2); 
     return;
}
 
void main()
{ float x,y,x1,x2,y1,y2,k,b,r;
  clrscr();
  cout<<"\nВведите координаты центра окружности:\n"<<"x=";
  cin>>x;
  cout<<"y=";
  cin>>y;
  cout<<"Введите радиус окружности:\n"<<"R=";
  cin>>r;
 
  cout<<"\nВведите координаты точек прямой:\n"<<"x1=";
  cin>>x1;
  cout<<"y1=";
  cin>>y1;
  cout<<"x2=";
  cin>>x2;
  cout<<"y2=";
  cin>>y2;
  k=(y1-y2)/(x1-x2);
  b=y1-k*x1;
  prover(x,y,k,b,r);
  getch();
}
При использовании обязательна ссылка на http://DMTSoft.ru