Jaja

Tarea de Práctica №2
Ejercicio №1
>> I=dsolve('DI=-0.19*(I-70)','I(0)=300')

I =

70+230*exp(-19/100*t)


>> ezplot('70+230*exp(-19/100*t)')>>

Ejercicio №2
>> dsolve('D2z=-9*Dz-25*z','z(0)=5, Dz(0)=0')

ans =

45/19*19^(1/2)*exp(-9/2*t)*sin(1/2*19^(1/2)*t)+5*exp(-9/2*t)*cos(1/2*19^(1/2)*t)Ejercicio №3
>> dsolve('D2z+Dz+12*z=20*cos(3*t)','z(0)=2, Dz(0)=0')

ans =-64/141*exp(-1/2*t)*sin(1/2*47^(1/2)*t)*47^(1/2)-4/3*exp(-1/2*t)*cos(1/2*47^(1/2)*t)+10/3*sin(3*t)+10/3*cos(3*t)


>> ezplot('-64/141*exp(-1/2*t)*sin(1/2*47^(1/2)*t)*47^(1/2)-4/3*exp(-1/2*t)*cos(1/2*47^(1/2)*t)+10/3*sin(3*t)+10/3*cos(3*t)')

Tareade Práctica №3
Ejercicio №2
function y=serie1(x)
global n;
y=0;
for i=1:n
y=y+pi/4+((1-(-1)^i)/(pi*i^2))*cos(n*x)+(1/n)*sin(n*x);
end

function y=funcion1(x)y=(pi-x)*heaviside(x);

global n;
ezplot('funcion1',[-pi/2,pi/2]);
hold on
m=15;
for n=1:m
ezplot('serie1',[-pi/2,pi/2]);
end

Ejercicio №3
N=20
S=0.02;X=rand(n,1)-0.5;
Y= rand(n,1)-0.5;
H=plot(x,y,’.’);
axis([-1 1 -1 1])
axis square
grid off
set (h,’erasemode’,’xor’,’’markersize’,18)
while 1
drawnow
x=x+s*randn(n,1);y=y+s*randn(n,1);
set(h,’xdata’,x,’ydata’,y)
end

Tarea de la práctica №4
function sised=tanques(t,x);sised=[-0.08*x(1)+0.1333*x(2)+20;0.08*x(1)-0.0666*x(2)+0.02*x(3);0.05333*x(2)-0.08*x(3)];

t0=0; tfinal=100;
x0=[5 7 6];
[t,x]=ode23('tanques',t0,tfinal,x0);
plot(t,x)

Universidad Nacional Autónoma de Honduras
(UNAH)

Tareade Ecuaciones Diferenciales

Catedrático: Lic. Joel Escoto

Alumno: Marcos Alfonso Hernández

№ de Cuenta: 20081902105

Sección: 0801

Fecha de Entrega: 04-12-2012