Programacion


Listado 2.1
x= 0:0.05:10;
y= sin(x).*exp (-0.4*x);
plot (x,y)
xlabel (‘x’); ylabel (‘y’)

Listado 2.2
x= (0:0.05:10)´;
y= sin(x).*exp (-0.4*x);
plot (x,y)
xlabel (‘x’); ylabel (‘y’)Listado 2.3
p= 0:0.05:8*pi;
z= (cos(p) + i*sin(2*p)).*exp(-0.05*p) + 0.01*p;
plot (real(z), imag(z))
xlabel (´Re(z)´); ylabel (‘Im(z)’)

Listado 2.4
x= (0:0.4:10)’;
y= sin(x).*exp (-0.4*x);plot (x,y,’+’)
xlabel (‘x’); ylabel (‘y’)
Listado 2.5
x= (0:0.2:10)’;
y= sin(x).*exp (-0.4*x);
plot (x,y)
grid on
xlabel (‘x’); ylabel (‘y’)

Listado 2.6
t= 0:0.05: pi+ 0.1;
y= sin(3*t).*exp(-0.3*t);
polar (t,y)
title(‘Grafica polar´)
grid

Listado 2.7
t= .1:.1:3;
x= exp(t);
y= exp(t.*sinh(t));
loglog(x,y)
grid
xlabel (‘x’); ylabel (‘y’)



Listado 2.28
t= .1:.1:3semilogy(t, esp(t.*t))
grid
xlabel (‘t’); ylabel (‘exp(t.*t)’);

Listado 2.9
t= .1:.1:3;
semilogx(t,exp(t.*t))
grid
xlabel (‘t’); ylabel (‘exp(t.*t)’);

Listado 2.10
x= 0:0.05:5;
y= sin(x);
z=cos(x);
plot (x,y, x,z)

Listado 2.11
x= 0:0.05:5;
y(1,: )= sin(x);
y(2,: )= cos(x);
plot (x,y)

Listado 2.12
x= (0:0.05:5)’;
y( :, 1)= sin(x);
y( :, 2)= cos(x);
plot (x,y)

Listado2.13
x= 0:0.05:5;
y= sin(x);
plot (x,y)
hold on
z= cos(x);
plot(x,z,’--‘)
xlabel (‘x’); ylabel (‘y(-)’, z(--)’);

Listado 2.14
clear; clf; hold off
x= 0:0.05:5;
y= sin(x);
plot (x,y)
holdon
z= cos(x);
plot(x,z)
hold off
Listado 2.15A
M= [0:0.01:1]’; k=1.4;
P0_entre_p= (1+(k-1)/2+M.^20).^(k/k-1));
plot (M, P0_entre_p )
xlabel(‘M, numero de Mach’)
title(‘ relacion de presion,p(estancamiento)p(estatica)’)

Listado 2.15B
% Relación de presión vs. Numero de March
Clear; clf; hold off;
M=[ 0:0.01:1]’;
K= 1.4;
P0_entre_p= (1+(k-1)/2+M.^2).^(k/(k-1));
hold onaxis(‘square’); % hace que la grafica sea cuadrada
plot(M, p0_entre_p)
xlabel(‘M, numero de Mach’)
ylabel(‘p0/p’)
title(‘Relacion de presión, p(estancamiento)/p(estatica)’)
text(0.45,0.01, ‘comprensible’)...