function circuit1(in)
    clc
  %  close all;
    global data;
    
    if nargin==0
        data=[]; 
        init;
    elseif nargin==1
        switch in
           %перечисление кнопок окна
            case 'Start', Start
        end
    end
    
    function init
       %название самого окна
        h.fig=figure('ToolBar','figure','NumberTitle','off','Name','Электрическая цепь второго порядка');
        %положение кнопки 
        uicontrol(h.fig,'style','pushbutton','String','Построить графики',...
            'Position', [70,10,150,50],'Callback','circuit1(''Start'')');
        
       %окна графиков
        h.previewI=axes('Parent',h.fig,'Position', [0.5,0.55,0.4,0.4]);
        title(h.previewI,'График тока на элементе');
        h.previewU=axes('Parent',h.fig,'Position', [0.5,0.05,0.4,0.4]);
        title(h.previewU,'График напряжения на элементе');
        h.graph=axes('Parent',h.fig,'Position', [0.05,0.65,0.4,0.4]);
         data.RGB=imread('схема2.jpg');%загрузка изображения схемы(изображение повернуто для корректной загрузки)
 
         
       image('Parent',h.graph,'CData',data.RGB)
        axis(h.graph,'off');
        
       
        %Поля ввода и надписи
         uicontrol(h.fig,'style','text','String','Параметры цепи:','Position', [10,250,100,20]);
         
        uicontrol(h.fig,'style','text','String','С,Ф:','Position', [10,190,100,20]);
        h.C=uicontrol(h.fig,'style','edit',...
            'String','0.003','Position', [110,190,100,20]);
        
        uicontrol(h.fig,'style','text','String','L,Гн:','Position', [10,170,100,20]);
        h.L=uicontrol(h.fig,'style','edit',...
            'String','1','Position', [110,170,100,20]);
        
        uicontrol(h.fig,'style','text','String','r1,Ом:','Position', [10,150,100,20]);
        h.r1=uicontrol(h.fig,'style','edit',...
            'String','40','Position', [110,150,100,20]);
        
        uicontrol(h.fig,'style','text','String','r2,Ом:','Position', [10,130,100,20]);
        h.r2=uicontrol(h.fig,'style','edit',...
            'String','100','Position', [110,130,100,20]);
        
         uicontrol(h.fig,'style','text','String','Ток I,А:','Position', [10,110,100,20]);
        h.I=uicontrol(h.fig,'style','edit',...
            'String','1','Position', [110,110,100,20]);
        
         uicontrol(h.fig,'style','text','String','элемент цепи:','Position', [10,90,100,20]);
         %переключатель 
        h.element=uicontrol(h.fig,'style','popup',...
            'String',{'r1' 'r2' 'L' 'C'},'Position', [110,90,100,20]);
        
        data.h=h;
        
%imagesc(X);
    end
    
    function Start
        h=data.h;
        %заргузка параметров цепи из окна
        r1=str2double(get(h.r1,'String'));
        C=str2double(get(h.C,'String'));
        L=str2double(get(h.L,'String'));
        r2=str2double(get(h.r2,'String'));
        I=str2double(get(h.I,'String'));
        %находим корни характеристического уравнения
        p1=-((L+r1*r2*C)/(2*r1*L*C))+sqrt(((L+r1*r2*C)/(2*r1*L*C))^2-(r1+r2)/(L*C*r1));
        p2=-((L+r1*r2*C)/(2*r1*L*C))-sqrt(((L+r1*r2*C)/(2*r1*L*C))^2-(r1+r2)/(L*C*r1));
        %находим постоянные
        A1=r2*p2*I/(p2-p1)/(r1+r2);
        A2=-r2*p1*I/(p2-p1)/(r1+r2);
        Iyst=I*r1/(r1+r2);%установившееся значение тока на катушке
        Tmax=max(abs(1/p1),abs(1/p2));
        t=linspace(0,Tmax,100);%вектор времени
        dT=Tmax/100;%шаг по времени
        %Ток и напряжение на элементах
        IL=Iyst+A1*exp(t.*p1)+A2*exp(t.*p2);
        Imax=Iyst+A1*exp(p1*(Tmax+dT))+A2*exp(p2*(Tmax+dT));
        IL1=[0 IL];
        IL2=[IL Imax];
        UL=L*(IL2-IL1)/dT;%Метод нахождения производной правой конечной разностью
        ULend=UL(end);
        UL=UL(2:end);
        
        U2=r2*IL;
        U2end=r2*Imax;
        I2=IL;
        
        UC=UL+U2;
        UCend=ULend+U2end;
        for i=1:(100-1)
            IC(i)=C*(UC(i+1)-UC(i))/dT;
        end;
        ICend=C*(UCend-UC(end))/dT;
        IC=[IC ICend];
        
        I1=UC/r1;
        U1=UC;
        
        %вывод графиков
        element=get(h.element,'Value');
        if element==1 
         plot(h.previewI,t,I1)
         title(h.previewI,'График тока на элементе');
         plot(h.previewU,t,U1)
         title(h.previewU,'График напряжения на элементе');
         
         elseif element==2
         plot(h.previewI,t,I2)
         title(h.previewI,'График тока на элементе');
         plot(h.previewU,t,U2)
         title(h.previewU,'График напряжения на элементе');
             
          elseif element==3
         plot(h.previewI,t,IL)
         title(h.previewI,'График тока на элементе');
         plot(h.previewU,t,UL)
         title(h.previewU,'График напряжения на элементе');     
              
                 elseif element==4
         plot(h.previewI,t,IC)
         title(h.previewI,'График тока на элементе');
         plot(h.previewU,t,UC)
         title(h.previewU,'График напряжения на элементе');
        end
        %вывод графиков в отдельное окно
        figure
        plot(t,IL)
        
        
     
%         cla(data.h.preview); image('Parent',data.h.preview,'CData',data.Data.RGB_)
%         cla(data.h.graph3D); temp=surf(data.h.graph3D,X,Y,Z);
%         set(temp,'LineStyle','none');
    end
 
end    
 
            
function out=const
    out.m=9.109383561111111111*10^-31;
    out.q=1.66*10^-19;
end