%Operaciojnes Matemáticas
5*sin(30*pi/180)*sqrt(-4)

%Ayuda
help sqrt

%Asignación de Variables
A=3
B=2
C=A^B

%Matrices
A=[1 2 3;4 5 6;7 8 9]
B=[2 4 6;8 10 12;14 16 18]
C=A+B
D=A*B

%Variables Simbólicas
s=sym('s')
res=(s+2)*(s+3)
res=expand((s+2)*(s+3))
pretty(res)
conv([1 2],[1 3])%Otra forma de multiplicar polinomios

%break
%Graficas
t=0:0.0001:5e-2
y=sin(2*pi*60*t)
plot(t,y)
plot(t,y,'.r')
title('Grafica Sinoidal')
xlabel('tiempo')
ylabel('Valores')

%Sistemas de Ecuaciones Lineales
%Ecuaciones de Malla
% 3 I1 +    8 I2    -   3 I3=0
% - I1  +    6 I2   +  5 I3= 20
%  I1   +     I2     -  10 I3=-8
%A*I=B     --->  I = linsolve(A,B)
[I] = linsolve([3 8 -3;-1 6 5;1 1 -10],[0 20 -8]')% traspuesta es: '
[x,y]=solve('5*x+2*s*y=3','2*x+s*y=0')


%Raices de un  polinomio
%8*s^2+3*s-8=0
raices=roots([8 3 -8])


%Ecuación Diferencial Lineal
%d2V/dt2+200dV/dt+250KV=1e7
V=sym('V')
t=sym('t')
solucion=dsolve('D2V+200*DV+250e3*V=1e7','V(0)=40','DV(0)=5000','t')
solucion=simplify(solucion)
tiempo=0:10e-6:100e-3;
senal=double(subs(solucion,t,tiempo));
figure
plot(tiempo,senal)
axis([0 80e-3 34 50])

whos
dir

%Simulación desde Simulink
sim('circuito1_curso')%stop time=0.5, Relative_Tolerance 1e-5,  type:variable_step,  ODE45
figure
plot(Vout.time,Vout.signals.values)

%Programción for, if, while
sum=0;
for i=2:2:20
    sum=sum+i;
end

%Mostrar avisos de Pantalla
disp(strcat('La suma es:',num2str(sum)))


%Cerrar Graficas
close all

%Borrar Variables
clear all

