amoniaco

Esta es una vista previa del archivo. Inicie sesión para ver el archivo original

clear
clc
% DATOS
v = [0.5 0.6 0.72 0.84 0.96 1.08 1.25]; % volumen
p = [1400 1248 1100 945 802 653 500]; % presión 
% Pregunta 1
W1 = (v(end) - v(1))*(p(end) + p(1))/2
% Pregunta 2
W2 = 0;
for i = 1:6
 W2 = W2 + (v(i+1) - v(i))*(p(i+1) + p(i))/2;
endfor
W2
% Pregunta 3
D = abs(W2 - W1)
% Pregunta 4
s = spline(v,p);
f = @(v) ppval(s,v);
n = (v(end) - v(1))/0.05;
W4 = trapecios(f,v(1),v(end),n)
% Pregunta 5
W5_0 = W4;
while true
 n = 2*n;
 W5_1 = trapecios(f,v(1),v(end),n);
 if abs(W5_0 - W5_1) <= 1e-2
 break;
 endif
 W5_0 = W5_1;
endwhile
W5 = W5_1