ООП: функция, которая решает квадратные уравнения.

classdef 1st
    properties
        D;
        X1;
        X;
    end
    methods
        function value = get.X(this)
            value = this.X; 
         end    
       
          function this = set.X(this,value)
             this.X = value; 
          end    
         
          function value = get.X1(this)
            value = this.X1; 
         end    
          function this = set.X1(this,value)
             this.X1 = value; 
          end    
 
    end
    methods
        function object=Untitled(varargin)
        A = input('Введите  А:');
        B = input('Введите  В:');
        C = input('Введите  С:');
        object.X1=(B+(B^2+4*A*C)^(0.5))/(2*A);
        object.X=(B-(B^2+4*A*C)^(0.5))/(2*A);
        fprintf('%s Корни уравнения: %s\n',object.X, object.X1);
        end
    end
       
end

uran
 
ans = 
 
  uran
 
  Properties:
     D: []
    X1: []
     X: []
 
  Methods

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

this.a=a;this.b=b;this.c=c;

classdef Kalk_1 < hgsetget
    properties
        A;
        B;
        C;
        D;
    end
    methods
         function value = get.A(this)
            value = this.A; 
         end    
          function set.A(this,value)
             this.A = value; 
          end    
         
          function value = get.B(this)
            value = this.B; 
         end    
          function  set.B(this,value)
             this.B = value; 
          end    
 
         function value = get.C(this)
            value = this.C; 
         end    
          function  set.C(this,value)
             this.C = value; 
          end    
         
          % Вычислимое свойство.Дискриминант
          function value = get.D(this)
             value = sqrt(this.B*this.B*this-4*this.A*this.C); 
          end 
    end
    methods
       function Vvod(this)
           this.A=input('Введите значение параметра А:');
       end
    end
 methods
       function Vvod2(this)
           this.B=input('Введите значение параметра B:');
       end
    end methods
       function Vvod1(this)
           this.C=input('Введите значение параметра C:');
       end
    end
    methods(Static)
        function [x1,x2]=Raschet(this)
 
if (this.A==0)&&(this.B==0)
    fprintf('Не квадратное уравнение')
    x1 = [];
    x2 = [];
else
if (this.A==0)&&(this.B~=0)
    fprintf('Не квадратное уравнение')
    x1 = -this.C/this.B;
    x2 = x1;
    formatSpec = 'Корень %4.2f';
    fprintf(formatSpec,x1)
else
if (this.D<0)
    fprintf('Дискриминант отрицательный');
end
this.D=sqrt(this.D);
x1 = (-this.B + this.D)/2*this.A;
x2 = (-this.B - this.D)/2*this.A;
formatSpec = 'Первый корень %4.2f второй %8.3f ';
fprintf(formatSpec,x1,x2)
end
end
x=((x1+x2)/2-10:0.1:(x1+x2)/2+10);
y=this.a*x.^2+this.b*x+this.c;
plot(x,y)
end
    end
    
end

roots([1 2 -3]) %где 1 2 -3 - коэффициенты, x^2+2x-3

global A B C D

roots([1 5 -6]) %где 1 5 -6 - коэффициенты, x^2+5x-6