METODOS NUMERICOS PRIMERA PRACTICA DE MATLAB
ALUMNO: OCADIO ACOSTA DENIS
i.a >> x=linspace(0,2*pi,200),plot(x,sin(x.^4))
ii.b >>x=linspace(0,2*pi,500),plot(x,sin(x.^4),'r'),grid on, axis equal
ii. >>x=0:pi/100:2*pi;plot(3*cos(x)-2,3*sin(x)+1,'r'), grid on, axis equal
iii. >>
x = -4:pi/100:4; y = (13-2*x)/3; y2 = 2*x-1; plot(x,y,x,y2 legend('(13-2*x)/3)','2*x-1')
iv.
>> hold on x=linspace(pi/2,5*pi/2,4);plot(cos(x),sin(x),'r'); x1=linspace(0,2*pi,100);plot(cos(x1),sin(x1),0.5*cos(x1),0.5*sin(x1),' b') hold off
v. hold
on x=linspace(0,2*pi,500);plot(2.5*cos(x),2.5*sin(x),'r*'); x1=linspace(0,2*pi,100); fill(0.3*cos(x1)+1,0.3*sin(x1)+1,'b',0.3*cos(x1)-1,0.3*sin(x1)+1,'r') x2=linspace(8.3*pi/6,9.7*pi/6,100);plot(4*cos(x2),4*sin(x2)+2.5,'g*') hold off
vi.
>> n=input('ingrese el valor de n: ') fprintf('la suma de los %g primeros números es: %8.0f\n',n,n*(n+1)/2) ingrese el valor de n: 5 n = 5 la suma de los 5 primeros números es: 15 >>
vii. >>
n=input('ingrese el valor de n: ') fprintf('la suma de cuadrados de los %g primeros números es: %8.0f\n',n*(n+1)*(2*n+1)/6) ingrese el valor de n: 10 n = 10 la suma de cuadrados de los 10 primeros números es: 385 >>
viii.
>> n=input('ingrese el valor de n: ') fprintf('la suma de los %g primeros números de esta forma es: %8.0f\n',n*(n+1)*(2*n+1)/6+n*(n+1)/2) ingrese el valor de n: 10 n = 10 la suma de los 10 primeros números de esta forma es: 440 >>
ix.
>> n=input('ingrese el valor de n: ') fprintf('la suma de los primeros %g números impares es: %8.0f\n',n*n) Ingrese el valor de n: 10 n = 10 la suma de los 10 primeros números impares es: 100 >>
x. >>
a=input('ingrese el valor del ancho: ') h=input('ingrese el valor de la altura: ') fprintf('el area del rectangulo es: %1.0f',a*h) x=linspace(0,a,200);y=linspace(0,h,200);y1=h;plot(x,0,'r*'); hold on;plot(x,y1,'r*');plot(0,y,'r*');plot(a,y,'r*'); grid on;title('AREA DE UN RECTANGULO'); xlabel('ANCHO');ylabel('ALTURA') axis square ingrese el valor del ancho: 8 a = 8 ingrese el valor de la altura: 2 h = 2 el área del rectángulo es: 16>>
xi.
a=input('ingrese el valor del primer lado: ') b=input('ingrese el valor del segundo lado: ') c=input('ingrese el angulo en sexagesimales: ') d=b*cos(c*pi/180) e=b*sin(c*pi/180) x=linspace(0,a,200);plot(x,0,'r*-'),axis equal; hold on; r=linspace(0,b,200);plot(r*cos(c*pi/180),r*sin(c*pi/180),'r*-'),axis equal; x1=d:pi/100:a;y=(e/(d-a))*(x1-a) ;plot(x1,y,'r*-'),axis equal; x2=a:pi/100:d;y=(e/(d-a))*(x2-a) ;plot(x2,y,'r*-'),axis equal grid on;title('AREA DE UN TRIANGULO'); axis square;
fprintf('el area del triangulo es = %8.0f\n',a*b*(sin(c*pi/180))/2) ingrese el valor del primer lado: 4 a = 4 ingrese el valor del segundo lado: 10 b = 10 ingrese el angulo en sexagesimales: 45 c = 45 d = 7.0711 e = 7.0711 el area del triangulo es = >>
14
xii. >> d1=input('ingrese
el valor del primer diagonal: ') d2=input('ingrese el valor del segundo diagonal: ') a=d1/2 b=d2/2 x=linspace(0,2*pi,5);plot(a*cos(x),b*sin(x)),axis equal; grid on;title('AREA DEl ROMBO'); xlabel('DIAGONAL 1');ylabel('DIAGONAL 2') axis square fprintf('el area del rombo es = %8.0f\n',d1*d2/2) ingrese el valor del primer diagonal: 8 d1 = 8 ingrese el valor del segundo diagonal: 10 d2 = 10 a = 4 b = 5 el area del rombo es =
>>
40
xiii.>> a=input('ingrese
el valor del # real: ') b=input('ingrese el valor del # imaginario: ') c=atan(b/a) r=sqrt(a.^2+b.^2) d=r*cos(c) e=r*sin(c) disp('el # en coord. pol. es'); z=[d,e]% ingrese el valor del # real: 3 a = 3 ingrese el valor del # imaginario: 4 b = 4 c = 0.9273 r = 5 d = 3.0000 e = 4.0000 el # en coord. pol. es z = 3.0000 >>
4.0000
xiv. >>
L=input('ingrese el valor del lado del cuadrado: ') fprintf('el area total es = %8.0f\n',2*L.^2) fprintf('el perimetro total es = %8.0f\n',L*4*(2+sqrt(2))) ingrese el valor del lado del cuadrado: 10 L = 10 el area total es = 200 el perimetro total es =
137
>>
xv.
>> a=input('ingrese dia y mes de nacimiento juntos: ') b=input('ingrese su edad: ') fprintf('el numero es = %8.0f\n',(2*a+5)*50+b+365-615) ingrese dia y mes de nacimiento juntos: 2704 a = 2704 ingrese su edad: 24 b = 24 el numero es = >>
270424