clc;
clear;
display('                   Programa para calcular diagrama de Corte y Momento en una viga        ');
display('Las condiciones de apoyo son las siguientes: Empotrado: 1 , Rotulado Deslizante: 2 ,  Rotulado Fijo: 3   ');
L=input('Ingrese largo de la viga [m]     ');
q=input('Ingrese el valor de la carga distribuida rectangular [ton]     ');
b=input('Ingrese el valor de la base en [m]    ');
h=input('Ingese el valor de la altura en [m]    ');
E=input('Ingrese el valor del módulo de elasticidad E en Ton/m2    ');
I=b*h^3/12;
comb=input('Ingrese la combinación de apoyos, estas son: 11,12,21,23,32     ');
display('Figura 1 muestra diagrama de Momento');
display('Figura 2 muestra diagrama de Corte');

if comb==11
    syms x
    M=-q*L/2*(x^2/L+L/6-x);     %Ecuación de momento para viga con doble empotrado
    V=q*L/2*(2*x/L-1);          %Ecuación de corte para viga con doble empotrado
    figure(1);
    ezplot(M,[0 L]);
    figure(2);
    ezplot(V,[0 L]);
    delta=q*L^4/(384*E*I);      %Ecuación de la deformada
    display(['La deformada máxima es de ' num2str(delta) '[m]=']);%Ecuación de la deformada
    
    
elseif comb==12
    syms x
    M=-1/8*q*L^2-q*x^2/2+5/8*q*L*x;      %Ecuación de momento para viga con empotrado - deslizante
    V=-q*x+5*q*L/8;                      %Ecuación de corte para viga con empotrado - deslizante
    figure(1);
    ezplot(M,[0 L]);
    figure(2);
    ezplot(V,[0 L]);
    delta=q*L^4/(185*E*I);               %Ecuación de la deformada
    display(['La deformada máxima es de ' num2str(delta) '[m]=']);

elseif comb==21
    
    syms x
    M=-q*L*(x^2/(2*L)-3*x/8);             %Ecuación de momento para viga con deslizante - empotrado
    V=q*L*(x/L-3/8);                      %Ecuación de corte para viga con empotrado - deslizante
    figure(1);
    ezplot(M,[0 L]);
    figure(2);
    ezplot(V,[0 L]);
    delta=q*L^4/(185*E*I);                %Ecuación de la deformada
    display(['La deformada máxima es de ' num2str(delta) '[m]=']);
    
elseif comb==23
    
    syms x
    M=-q*x^2/2+q*L*x/2;              %Ecuación de momento para viga con rotulado deslizante - rotulado fijo
    V=-q*x+q*L/2;                    %Ecuación de corte para viga con rotulado deslizante - rotulado fijo
    figure(1);
    ezplot(M,[0 L]);
    figure(2);
    ezplot(V,[0 L]);
    delta=5*q*L^4/(384*E*I);         %Ecuación de la deformada
    display(['La deformada máxima es de ' num2str(delta) '[m]=']);
        
elseif comb==32
    
    syms x
    M=q*L/2*(x-x^2/L);               %Ecuación de momento para viga con rotulado fijo - rotulado deslizante
    V=-q*L/2*(1-2*x/L);              %Ecuación de corte para viga con rotulado fijo - rotulado deslizante
    figure(1);
    ezplot(M,[0 L]);
    figure(2);
    ezplot(V,[0 L]);
    delta=5*q*L^4/(384*E*I);        %Ecuación de la deformada
    display(['La deformada máxima es de ' num2str(delta) '[m]=']);
   
end