Bài tập số 1 (Chương 3)

#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <process.h>

float a,b,c,x12, x1, x2, X1, X2, X3, X4, Y1, Y2, delta;
int chon;
int menu()
{
  int so;
  printf("\n---------Menu---------\n");
  printf("\n1. Giai phuong trinh ax^2 + bx + c = 0");
  printf("\n2. Giai phuong trinh aX^4 + bX^2 + c = 0");
  printf("\n0. Thoat");
  printf("\n----------------------");
  printf("\nChon cong viec: ");
  do
    {
      so = getch();

    }
  while (so>'2'||so<'0');
  putc(so, stdout);
  return so-'0';
}
void nhap()
{
  printf("\nNhap vao cac he so cua phuong trinh");
  do { printf("\n  a = "); scanf("%f",&a);
  } while (a==0)  ;
  printf("\n  b = ");     scanf("%f",&b);
  printf("\n  c = ");     scanf("%f",&c);
}
void hienthi()
{
  if (delta<0) { printf("\nPhuong trinh vo nghiem"); getch(); return;}
  if (chon==1)
    if (delta==0)
      if (c==0) printf("\nPhuong trinh co nghiem x = 0");
      else printf("\nPhuong trinh co nghiem kep x = %f",x12);
    else printf("\nPhuong trinh co 2 nghiem phan biet: \nx1 = %f\nx2 = %f",x1,x2);
  else
    if (delta==0)
      if (x12<0)  { printf("\nPhuong trinh vo nghiem"); getch(); return;}
   else if (x12==0)  printf("\nPhuong trinh co nghiem X = 0");
     else printf("\nPhuong trinh co 2 nghiem phan biet: \nX1 = %f\nX2 = %f",Y1,Y2);
    else
      if (x1*x2==0)
      {
   if ((x1<0)||(x2<0)) printf("\nPhuong trinh co nghiem X = 0\n");
   if ((x1>0)||(x2>0)) printf("\nPhuong trinh co 3 nghiem phan biet: \nX1 = 0\nX2 = %f\nX3 = %f",sqrt(x1+x2), -sqrt(x1+x2));
      }
      if (x1*x2>0)
   if (x1<0)  { printf("\nPhuong trinh vo nghiem"); getch(); return;}
   else printf("\nPhuong trinh co 4 nghiem phan biet:\nX1 = %f\nX2 = %f\nX3 = %f\nX4 = %f",X1,X2,X3,X4);
      if (x1*x2<0)
   if (x1<0)  printf("\nPhuong trinh co 2 nghiem phan biet:\nX1 = %f\nX2 = %f",X3,X4);
   else printf("\nPhuong trinh co 2 nghiem phan biet:\nX1 = %f\nX2 = %f",X1,X2);
  getch();
}
void ptbac2()
{
    delta=b*b-4*a*c;
    if (delta==0) x12=-b/2/a;
    if (delta>0) { x1=(-b+sqrt(delta))/2/a;  x2=(-b-sqrt(delta))/2/a; }
}
void ptbac4()
{
    ptbac2();
    if (x12>0) { Y1=sqrt(-b/a); Y2=-sqrt(-b/a); }
    if (x1>=0) {X1=sqrt(x1); X2=-sqrt(x1);}
    if (x2>=0) {X3=sqrt(x2); X4=-sqrt(x2);}
}
void main()
{
  do {
    clrscr();
    chon=menu();
    switch (chon)
    {
      case 0 : exit(0);
      case 1 : nhap();  ptbac2();    break;
      case 2 : nhap();  ptbac4();    break;
    }
    hienthi();
  }  while (1);

}