1.
HALLAR:
1−42 1 −4 2
1 1 1 1 1 1 1 1
I = ∭ ln 2+
−1−21
( )(z x− y )
dx dy dz=∫ ln 2+ (
−1 z x( )
− y
)dx ∫ ln 2+ (
−2 z x− y ( )
) dy ∫ ln 2+ (
1 z x− y
) dz ( )
SOLUCIÓN:
clc,clear
syms x y z
a=1;
b=2;
n=2;
h=(b-a)/n;
x1=a;
x2=a+h;
x3=x2+h;
fx1=log(2+1/z)*(1/(x1-y));
fx2=log(2+1/z)*(1/(x2-y));
fx3=log(2+1/z)*(1/(x3-y));
I1=h/3*(fx1+4*fx2+fx3);
a=-2;
b=-4;
n=2;
h=(b-a)/n;
y1=a;
y2=a+h;
y3=y2+h;
fy1=eval(fx1);
fy2=eval(fx2);
fy3=eval(fx3);
I2=h/3*(fy1+4*fy2+fy3);
a=-1;
b=1;
n=2;
h=(b-a)/n;
z1=a;
z2=a+h;
z3=y2+h;
fz1=eval(fy1);
fz2=eval(fy2);
fz3=eval(fy3);
I3=h/3*(fz1+4*fz2+fz3);
I=I1*I2*I3;
disp('INTEGRAL:')
disp(I)
INTEGRAL:
(Log (1/z + 2)/ (3*(y - 1)) + log (1/z + 2)/ (3*(y - 2)) + (4*log (1/z + 2))/ (3*(y - 3/2)))^2*(log (1/z
+ 2)/ (6*(y - 1)) + log (1/z + 2)/ (6*(y - 2)) + (2*log (1/z + 2))/ (3*(y - 3/2)))
2. Sea el sistema mecánico:
Las ecuaciones del movimiento son:
J 1 θ1 '' + ( k 1 +k 2 ) θ1 −k 2 θ 2= 0
J 2 θ2 ''− k 2 θ1 + k 2 θ2 =0
k 2 =2 k 1 =1
J 2 =2 J 1 =1 Se sabe que en el instante t=0 s, los discos
están en la posición de θ1 =0,1 rad y θ2 =0,2 rad , la velocidad angular inicial
de ambos discos es sin movimiento =0 rad/s. Aproximar las posiciones y
velocidades angulares de los discos en t=0,45 seg. Utilizar el método de Euler
modificado.
SOLUCION:
clc;
clear all;
syms t x y w z;
disp(' MÉTODO EULER MODIFICADO');
disp(' *************************');
disp('*De las siguientes EDOs:');
disp(' dw=(y-(1.5*x))/0.5');
disp(' dx=w');
disp(' dz=x-y');
disp(' dy=z;');
disp('*Realizando la iteración:');
dw=(y-(1.5*x))/0.5;
dx=w;
dz=x-y;
dy=z;
t0=0;
tf=0.45;
wi=0;
zi=0;
xi=0.1;
yi=0.2;
n=10;
h=(tf-t0)/n;
disp('--------------------------------------------------------------')
;
fprintf('i t(i) w z Teta1 Teta2\n');
disp('--------------------------------------------------------------')
;
fprintf('%2d %10.6f %10.6f %10.6f %10.6f %10.6f\n',i,t0,wi,zi,xi,yi);
while i<n;
t=t0; w=wi; z=zi; x=i; y=yi;
wwi=wi+h*eval(dw);
wi=wi+(h/2)*(eval(dw)+double(wwi));
xxi=xi+h*eval(dx);
xi=xi+(h/2)*(eval(dx)+double(xxi));
zzi=zi+h*eval(dz);
zi=zi+(h/2)*(eval(dz)+double(zzi));
yyi=yi+h*eval(dy);
yi=yi+(h/2)*(eval(dy)+double(yyi));
t0=t0+h;
i=i+1;
fprintf('%2d %10.6f %10.6f %10.6f %10.6f
%10.6f\n',i,t0,wi,zi,xi,yi);
end
disp('--------------------------------------------------------------')
3. Resolver el siguiente modelo matemático con el método implícito
∂C ∂2 C ∂C
=D( 2 + )
∂x ∂Z ∂Z
Establecer las condiciones necesarias para resolver la EDP considerando que
X representa el tiempo y z la dirección de la variación de la concentración C.
D=1.7*10^-2 cm^2/h
SOLUCION:
function inplicito_sinderivada
clear
D=1.7*10^-2; %cm^2/h
m=5;
n=10;%
Ta=10;%u(0,t)=ta
Tb=5;%u(L,t)=tb
To=0;%u(x,0)=t0
dx=0.25;%
dt=0.01;%
L=1;%
k=1;%
clear x u t
t(1)=Ta;
t(m)=Tb;
for i=2:D-1
t(i)=To;
end
landa=k*dt/(dx^2);%
r=landa;%
q=-1-2*landa;%
p=landa;%
hold on
title('curva de de concentracion');
xlabel('x distancia');
ylabel('t tiempo');
x=0:dx:L;
plot(x,t,'r')
for j=1:n
clear a b c d
for i=1:m-2
a(i)=p;
b(i)=q;
c(i)=r;
d(i)=-1*t(i+1);
end
d(1)=d(1)-a(1)*t(1);
d(m-2)=d(m-2)-c(m-2)*t(m);
n=length(d);
w(1)=b(1);
g(1)=d(1)/w(1);
for i=2:n
w(i)=b(i)-a(i)*c(i-1)/w(i-1);
g(i)=(d(i)-a(i)*g(i-1))/w(i);
end
u(n)=g(n);
for i=n-1:-1:1
u(i)=g(i)-c(i)*u(i+1)/w(i);
end
u;
t=[t(1) u t(m)]
if mod(j,5)==0
plot(x,t,'r');
end
end
inplicito_sinderivada
t = 10.0000 1.2396 0.2271 0.6336 5.0000
t = 10.0000 2.2217 0.5815 1.1565 5.0000
t = 10.0000 3.0165 1.0005 1.6035 5.0000
t = 10.0000 3.6725 1.4451 1.9960 5.0000
t = 10.0000 4.2236 1.8912 2.3474 5.0000
t = 10.0000 4.6936 2.3249 2.6662 5.0000
t = 10.0000 5.0998 2.7379 2.9578 5.0000
t = 10.0000 5.4545 3.1263 3.2258 5.0000
t = 10.0000 5.7672 3.4884 3.4727 5.0000
t = 10.0000 6.0447 3.8240 3.7004 5.0000
