Квадратное уравнение

function [x1,x2]=korni(a,b,c)
 
if (a==0)&&(b==0)
    fprintf('eto ne yravnenie!')
    x1 = [];
    x2 = [];
else
if (a==0)&&(b~=0)
    fprintf('yravnenie ne kvadratnoe')
    x1 = -c/b;
    x2 = x1;
else
 
d = (b^2 - 4*a*c);
if (d<0)
    fprintf('discriminant otricatelnui, korni kompleksnue');
end
d=sqrt(d);
x1 = (-b + d)/2*a;
x2 = (-b - d)/2*a;
end
end
x=[(x1+x2)/2-10:0.1:(x1+x2)/2+10];
y=a*x.^2+b*x+c;
plot(x,y)
end

function[x1,x2] = root2 (a,b,c)
if (nargin<3)
error('задайте 3 коэффициента')
end
[Na, Ma]=size(a);
[Nb,Mb]=size(b);
[Nc,Mc]=size(c);
if (Na~=1)|(Ma~=1)|(Nb~=1)|(Mb~=1)|(Nc~=1)|(Mc~=1)
error(' аргументы должны быть числами')
end
D=b^2-4*a*c;
if D<0
error('комплексные корни')
end
x1=(-b+sqrt(D))/2;
x2(-b-sqrt(D))/2;
if(nargout ~=2
warning('это только один корень')
disp('для получения 2 корней используйте')
disp('[x1,x2]=root2(a,b,c)')
end

 
function [x1,x2]=korni(a,b,c)

[x1,x2]=korni(a,b,c)

function [x1,x2]=korni(a,b,c)
...
end

function [x1,x2]=korni(a,b,c)

function[x1,x2] = root2 (a,b,c)
 
if (nargin<3)
error('задайте 3 коэффициента')
end
[Na, Ma]=size(a);
[Nb,Mb]=size(b);
[Nc,Mc]=size(c);
if (Na~=1)|(Ma~=1)|(Nb~=1)|(Mb~=1)|(Nc~=1)|(Mc~=1)
error(' аргументы должны быть числами')
end
D=b^2-4*a*c;
if D<0
error('комплексные корни')
end
x1=(-b+sqrt(D))/2;
x2=(-b-sqrt(D))/2;
if  nargout ~= 2
warning('это только один корень')
disp('для получения 2 корней используйте')
disp('[x1,x2]=root2(a,b,c)')
end
end

>> [x1,x2] = root2(1,0,-4)
x1 =
     2
x2 =
    -2