1.
- HACER UN PROGRAMA QUE PERMITA DETERMINAR SI UN NÚMERO N ENTERO POSITIVO ES PERFECTO, SI NO LO ES, ENCONTRAR EL NÚMERO PERFECTO MÁS
CERCANO, INDICANDO CUANTO SE TIENE QUE RESTAR O SUMAR PARA HALLAR DICHO NÚMERO PERFECTO.
clc , clear all , format compact
disp('=================================')
n=input('INGRESE EL NUMERO = ');
disp('=================================')
while n<=1
n=input('INGRESE EL NUMERO = ');
end
D=n;
sum=0;
for d=1:n;
resto=rem(D,d);
if resto==0;
sum=sum+d;
disp([d])
end
end
if sum==2*d;
disp( 'PERFECTO')
else
disp('IMPERFECTO')
end
disp('=================================')
2. Hacer un programa para determinar la cantidad de números curiosos existentes en el rango 0 y
100000 ó 0 y 1000000
NOTA: Número curioso es: todo número natural n que cumple que n2 tiene al propio n como
última cifra. Por ejemplo, 25 y 36 son números curiosos.
clc, clear all, format compact
disp('==========================================================')
disp('¿CUALES SON LOS NUMEROS CURIOSOS QUE HAY ENTRE 0 A 1000000?')
disp('Son...')
disp('==========================================================')
for t=1:10;
r=t.^2;
i(r)=mod(r,10);
if i(r)==t;
s=(i(r)).^2;
disp(s)
end
end
for t=11:100;
r=t.^2;
i(r)=mod(r,100);
if i(r)==t;
s=(i(r)).^2;
disp(s)
end
end
for t=101:1000;
r=t.^2;
i(r)=mod(r,1000);
if i(r)==t;
s=(i(r)).^2;
disp(s)
end
end
disp('===========================================================')
1. Ejemplos: 22 octubre 2001 es correcta; 31 abril 2001 no es correcta porque abril no tiene 31 días;
29 de febrero de 2001 no es correcto porque el año 2001 no es bisiesto.
clc,clear all, format compact
d=input('Ingrese el número de día:')
m=input('Ingrese el mes:')
a=input('Ingrese el año:')
switch m
case 1
m='ENERO'
for d=1:31
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & (a/100)~=0
disp ('El año es bisiesto')
else if a/100==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 2
m='FEBRERO'
if (d>1 & d<29) | (d>1 & d<30)
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 3
m='MARZO'
for d=1:31
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 4
m='ABRIL'
for d=1:30
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 5
m='MAYO'
for d=1:31
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 6
m='JUNIO'
for d=1:30
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 7
m='JULIO'
for d=1:31
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 8
m='AGOSTO'
for d=1:31
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 9
m='SETIEMBRE'
for d=1:30
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 10
m='OCTUBRE'
for d=1:31
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 11
m='NOVIEMBRE'
for d=1:30
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 &a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
case 12
m='DICIEMBRE'
for d=1:31
end
if a>999 & a<10000
r=rem(a,4)
if r==0 & a/100~=0
disp ('El año es bisiesto')
else if a/100 ==0 & a/400~=0
disp('El año no es bisiesto')
else
disp('El año no es bisiesto')
end
end
else
disp('Introduce un número de 4 dígitos')
end
end
clc
clear all
disp('-----------------------------------')
disp('PROGRAMA PARA HALLAR LA ALTURA MAX.')
disp('-----------------------------------')
vi=input('ingrese la velocidad inicial:');
x=input('ingrese angulo:');
g=9.8;
hmax=(vi^2*(sin(x))^2)/2*g;
for h=(0:1:hmax);
v=((h*2*g)/(sin(x))^2)^0.5;
disp([h v])
end
clc, clear all;
disp('---------------------------------------------------------------------');
disp('-----PROGRAMA QUE IMNGRESE LOS COEFICIENTES DE LA ECUACION-----------');
disp('-------------------DE LA FORMA X=At^3+Bt^2+Ct+d----------------------');
disp('---------------------------------------------------------------------');
disp('despalsamiento:X=At^3+Bt^2+Ct+d ')
disp('velocida: (d(x)/d(t)) v=3At^2+2Bt+C')
disp('aceleracion: (d(v)/d(t)) a=6At+2B')
A=input('coef... A: ');
B=input('coef... B: ');
C=input('coef... C: ');
D=input('coef... D: ');
disp(' T(s) D(m) V(m/s) A(m/s2)')
t=0;
format short
for t=1:10;
x=A*t^3+B*t^2+C*t+D;
v=3*A*t^2+2*B*t+C;
a=6*A*t+2*B;
disp ([ t , x, v, a ])
subplot(2,2,1),plot(t,x,'r*'),title('grafica D vs T'),xlabel('T(s)'),ylabel('D(m)'),hold on;
subplot(2,2,2),plot(t,v,'b*'),title('grafica V vs T'), xlabel('T(s)'),ylabel('V(m/s)'),hold on;
subplot(2,2,3),plot(t,a,'g*'),title('grafica A vs T'), xlabel('T(s)'),ylabel('A(m/s2)'),hold on;
end
clc, clear all;
disp('--------------------==============-')
disp('PROGRAMA PARA CONVERTIR TEMPERATURA')
disp('--------------------==============-')
temp=input('ingrese la temperatura:');
disp 'Centigrado - Kelvin - farenheit'
for x=temp-5:1:temp+5;
kelvin=x+273;
Farenheit=(x-32)/9;
disp([x kelvin Farenheit])
end
clc
clear all
disp('-------------------------')
disp('cuadro de multiplicacion')
disp('-------------------------')
for n=1:1:12
n;
m=n*1;
p=n*2;
q=n*3;
s=n*4;
t=n*5;
disp([n,m,p,q,s,t])
end
clc
clear all
disp('----------------------------------------')
disp('programa que produce graficos de petalos')
disp('----------------------------------------')
theta=-pi:0.01:pi;
rho(1,:) = 2*sin(5*theta).^2;
rho(2,:) = cos(10*theta).^3;
rho(3,:) = sin(theta).^2;
rho(4,:) = 5*cos(3.5*theta).^3;
for i = 1:4
polar(theta,rho(i,:))
pause
end
clc
clear all
disp('------------------------------------------------------------------------------------------------------------')
disp('PROGRAMA PARA CONVERTIR GRADOS CELCIUS A KELVIN Y FARENHEIT')
disp('------------------------------------------------------------------------------------------------------------')
disp(' celcius kelvin farenheit')
for C=1:10;
K=C+273;
F=9/5*C+32;
disp([C K F])
end
clc
clear all
disp('------------------------------------------------------')
disp('programa para hallar factoriales hasta n')
disp('------------------------------------------------------')
e=input('ingrese el ultimo factorial:');
N=0;
fact=1;
disp(' n fact')
disp([N fact])
for N=1:e
fact=fact*N;
disp([N fact])
end
clc
clear all
disp('----------------------------------------------------------------------')
disp('^^^^converción de grados centrigrados a grados kelvin y Farenheit^^^^')
disp('----------------------------------------------------------------------')
disp('celsius kelvin farenheit')
for n=1:10;
C=n^2; %muestra los cuadrados de los numeros
K=n+273; % para convertir en grados Kelvin
F=(9*n+150)/5; %para convertir en grados Farenheit
disp([n,C,K,F])
end
clc
clear all
disp('----------------------------------------')
disp(' SUMA DE LOS NUMEROS NATURALES HASTA M')
disp('----------------------------------------')
M=input('ingrese numero:');
s=0;
disp(' M SUAMATORIA ')
for n=1:M
s=n+s;
disp([n,s])
end
clc
clear all
disp('-----------------------------')
disp(' TABLA DE SUMA ')
disp('-----------------------------')
n=input('ingrese n:');
disp(' n numeros resultado ')
for i=0:10;
s=n+i;
disp([n,i,s])
end
clc
clear all
disp('------------ -----')
disp(' TABLA DE RESTA ')
disp('------------------')
n=input('ingrese el numero:');
disp('numeros n resultado ')
for m=n:n+10;
d=m-n;
disp([m n d])
end
clc
clear all
disp('-------------------')
disp('metodo del trapecio')
disp('--------------------')
x=input('INGRESE VALOR X')
y=input('INGRESE VALOR F(X)')
a=input('ingrese limite inferior')
b=input('ingrese limite superior')
n=length(x)
N=n-1
h=(b-a)/N
suma=0
g=zeros(n,1);
for i=1:n
g(i)=polyval(y,x(i));
suma=suma + g(i)
end
area=(h/2)*(2*suma-g(1)-g(n))