Matlab - PROBLEMA CON GUIDE EN UNA INTERFAZE GRAFICA!!! AYUDA!!

function varargout = heun(varargin)

% HEUN MATLAB code for heun.fig

%      HEUN, by itself, creates a new HEUN or raises the existing

%      singleton*.

%

%      H = HEUN returns the handle to a new HEUN or the handle to

%      the existing singleton*.

%

%      HEUN('CALLBACK',hObject,eventData,handles,...) calls the local

%      function named CALLBACK in HEUN.M with the given input arguments.

%

%      HEUN('Property','Value',...) creates a new HEUN or raises the

%      existing singleton*.  Starting from the left, property value pairs are

%      applied to the GUI before heun_OpeningFcn gets called.  An

%      unrecognized property name or invalid value makes property application

%      stop.  All inputs are passed to heun_OpeningFcn via varargin.

%

%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one

%      instance to run (singleton)".

%

% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help heun

% Last Modified by GUIDE v2.5 20-May-2014 23:57:03

% Begin initialization code - DO NOT EDIT

gui_Singleton = 1;

gui_State = struct('gui_Name',       mfilename, ...

                   'gui_Singleton',  gui_Singleton, ...

                   'gui_OpeningFcn', @heun_OpeningFcn, ...

                   'gui_OutputFcn',  @heun_OutputFcn, ...

                   'gui_LayoutFcn',  [] , ...

                   'gui_Callback',   []);

if nargin && ischar(varargin{1})

    gui_State.gui_Callback = str2func(varargin{1});

end

if nargout

    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});

else

    gui_mainfcn(gui_State, varargin{:});

end

% End initialization code - DO NOT EDIT

% --- Executes just before heun is made visible.

function heun_OpeningFcn(hObject, eventdata, handles, varargin)

% This function has no output args, see OutputFcn.

% hObject    handle to figure

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

% varargin   command line arguments to heun (see VARARGIN)

% Choose default command line output for heun

handles.output = hObject;

% Update handles structure

guidata(hObject, handles);

% UIWAIT makes heun wait for user response (see UIRESUME)

% uiwait(handles.figure1);

% --- Outputs from this function are returned to the command line.

function varargout = heun_OutputFcn(hObject, eventdata, handles) 

% varargout  cell array for returning output args (see VARARGOUT);

% hObject    handle to figure

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure

varargout{1} = handles.output;

function cajan_Callback(hObject, eventdata, handles)

% hObject    handle to cajan (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of cajan as text

%        str2double(get(hObject,'String')) returns contents of cajan as a double

% --- Executes during object creation, after setting all properties.

function cajan_CreateFcn(hObject, eventdata, handles)

% hObject    handle to cajan (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.

%       See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

    set(hObject,'BackgroundColor','white');

end

function cajato_Callback(hObject, eventdata, handles)

% hObject    handle to cajato (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of cajato as text

%        str2double(get(hObject,'String')) returns contents of cajato as a double

% --- Executes during object creation, after setting all properties.

function cajato_CreateFcn(hObject, eventdata, handles)

% hObject    handle to cajato (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.

%       See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

    set(hObject,'BackgroundColor','white');

end

function cajatf_Callback(hObject, eventdata, handles)

% hObject    handle to cajatf (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of cajatf as text

%        str2double(get(hObject,'String')) returns contents of cajatf as a double

% --- Executes during object creation, after setting all properties.

function cajatf_CreateFcn(hObject, eventdata, handles)

% hObject    handle to cajatf (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.

%       See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

    set(hObject,'BackgroundColor','white');

end

function condiciones_Callback(hObject, eventdata, handles)

% hObject    handle to condiciones (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of condiciones as text

%        str2double(get(hObject,'String')) returns contents of condiciones as a double

% --- Executes during object creation, after setting all properties.

function condiciones_CreateFcn(hObject, eventdata, handles)

% hObject    handle to condiciones (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.

%       See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

    set(hObject,'BackgroundColor','white');

end

function funcion_Callback(hObject, eventdata, handles)

% hObject    handle to funcion (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of funcion as text

%        str2double(get(hObject,'String')) returns contents of funcion as a double

% --- Executes during object creation, after setting all properties.

function funcion_CreateFcn(hObject, eventdata, handles)

% hObject    handle to funcion (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.

%       See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

    set(hObject,'BackgroundColor','white');

end

% --- Executes on button press in pushbutton1.

function pushbutton1_Callback(hObject, eventdata, handles)

% hObject    handle to pushbutton1 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

syms y

 n=str2double(get(handles.cajan,'string'));

 t0=str2double(get(handles.cajato,'string'));

 tf=str2double(get(handles.cajatf,'string'));

for i=1:n

    y0(i)=str2double(get(handles.condiciones,'string'));    %AQUI QUIZAS SEA EL PROBLEMA

end

 fii=get(handles.funcion,'string');

 fi=inline(fii,'t','y');

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 %sistemas de EDO! Metodo de Heun (RK 23)

[X, Y] = ode23(fi,[t0 tf], y0);

[X,Y]

figure(1)

plot (X,Y)

title('Solucion con Metodo de Heun')

xlabel('intervalos del tiempo final')

grid on

if n==2

    figure(2)

    plot3(X(:,1),Y(:,1),Y(:,2))

    title('solucion en 3D metodo de Heun')

    xlabel('Intervalos del Tiempo en x')

    grid on

    figure(3)

    plot(X(:,1),Y(:,1))

    title('solucion en metodo de Heun solucion de X vs Y1')

    xlabel('Intervalos del Tiempo en x')

    grid on

    figure(4)

    plot(X(:,1),Y(:,2))

    title('solucion en metodo de Heun solucion de X vs Y2')

    xlabel('Intervalos del Tiempo en x')

    grid on

    figure(5)

    plot(Y(:,1),Y(:,2))

    title('solucion en metodo de Heun solucion de Y1 vs Y2')

    xlabel('Intervalos del Tiempo en x')

    grid on

end

if n == 3

%con una dimension mayor no sale la grafica

figure(2)

plot3(Y(:,1),Y(:,2),Y(:,3))

title('solucion en 3D metodo de Heun')

xlabel('Intervalos del Tiempo en x')

grid on

%en esta seccion podemos ver la interaccion con cada variable

figure(3)

plot (Y(:,1),Y(:,2))

title('Solucion con Metodo de Heun solucion de Y1 vs Y2')

xlabel('intervalos del tiempo final')

grid on

figure(4)

plot (Y(:,1),Y(:,3))

title('Solucion con Metodo de Heun solucion de Y1 vs Y3')

xlabel('intervalos del tiempo final')

grid on

figure(5)

plot (Y(:,2),Y(:,3))

title('Solucion con Metodo de Heun solucion de Y2 vs Y3')

xlabel('intervalos del tiempo final')

grid on

end

if n>3

        nn=2;

    for i=1:length(X)

      for j=1:length(Y)

          for k=1:length(Y)

             figure(nn)

             p(nn,:)=plot([X(:,i),Y(:,j),Y(:,k)]');

             title('soluciones con Metodo de Heun entre las variables combinadas')

             xlabel('paso del tiempo x')

             grid on

             nn=nn+1;

          end

      end

    end

end

ans =

     0   NaN   NaN   NaN

    20   NaN   NaN   NaN

    40   NaN   NaN   NaN

    60   NaN   NaN   NaN

    80   NaN   NaN   NaN

   100   NaN   NaN   NaN

   120   NaN   NaN   NaN

   140   NaN   NaN   NaN

   160   NaN   NaN   NaN

   180   NaN   NaN   NaN

   200   NaN   NaN   NaN

clc

clear all

disp('----------------------------')

disp(' Metodo de Heun     ')

disp('----------------------------')

disp('Introdusca la Dimension del Sistema, n')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

n=input('Introdusca la dimension del sistema es')

x0=zeros(n,1) %genera el vector x0 con la dimension n

t0=input('Tiempo Inicial To')

for i=1:n

    fprintf('Valor o condicion inicial de la variable',i)

    x0(i)=input('xi0=')

end

N=input('tamaño de paso tao')

T=input('tiempo maximo T=')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

disp('Introdusca el Sistema de EDO en forma de x(i)')

disp(' ')

SIS_DIF=cell(n,1) %es la ()

disp(' ')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for i=1:n

    fprintf('Lado derecho de la Ecuacion t(1)=x,x(1)=y,x(2)=z',i)

    fii=input('fi=','s'); %el s es porque es una cadena de caracter

    fi=inline(fii,'t','x');

    SIS_DIF{i}=fi % es un {}

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

metodos=input('introduce el tipo de metodo que usted quiere,1),2) o 3)')

      %Metodo de Heun

        t=t0:N:T

        y=zeros(n,length(t));

        y(:,1)=x0;

        h=(T-t0)/N;

        for k=1:length(t)

            for i=1:n

                fi=SIS_DIF{i};

                yp(i,k+1)= y(i,k)+h*fi(t(k),y(i,k));

                y(i,k+1)=y(i,k+1)+(h/2)*(fi(t(k),y(i,k))+fi(t(k+1),yp(i,k+1)));

            end

        end