%  脚本文件：calc_roots.m
% 
% 目标：
% 该程序解决了二元一次方程组：a*x^2 + b*x + c =0的根的问题
% 
% 版本记录
%     日期           编者          描述
%    =====         =========     ================
% 2015-10-2 20:10     泡泡         源码
% 
% 定义变量：
% 
% a            --x^2的参数
% b            --x的参数
% c            --方程常量参数
% discriminant --方程判别式变量
% imag_part    --方程虚根部分
% real_part    --方程实根部分
% x1           --方程实根1
% x2           --方程实根2
% 
%清除变量或指令
clc;
% 允许用户输入三个参数
disp ('该功能是求解方程的根 ');
disp ('方程：A*X^2 + B*X + C = 0. ');
a = input ('输入参数 A: ');
b = input ('输入参数 B: ');
c = input ('输入参数 C: ');

%计算判别式
discriminant = b^2-4*a*c;

%根据判别式分情况讨论
if discriminant >0  %两个实根
     x1 = (-b+sqrt(discriminant ))/(2*a)
     x2 = (-b-sqrt(discriminant ))/(2*a)
     disp('二元一次方程有两个不相等实根！ ');
     fprintf (' x1=  %f\n' ,x1);
     fprintf (' x2=  %f\n' ,x2);
elseif discriminant == 0
     x1 =  -b/(2*a)
     disp('二元一次方程有两个相等实根！ ');
     fprintf (' x1=x2=  %f\n' ,x1);
else 
     real_part = -b/(2*a)
     image_part = sqrt( abs(discriminant) )/(2*a)
     disp('二元一次方程有两个复数根 ');
     fprintf (' x1=  %f +i %f\n' ,real_part ,image_part );
     fprintf (' x2=  %f -i %f\n' ,real_part ,image_part );
end