prg calculo

Click here to load reader

Post on 26-Oct-2014

135 views

Category:

Documents

0 download

Embed Size (px)

TRANSCRIPT

UNIVERSIDAD NACIONAL DEL CENTRO DEL PER

FACULTAD DE INGENIERA QUMICA

LENGUAJES DE PROGRAMACION

CATEDRTICO: PRESENTADO POR:

Ing. LPEZ GUTIERREZ, Helmer MARAVI RIOS, Ramiro Alex HUANCAYO PERU MARZO-2009

DEDICATORIAEl presente trabajo lo dedico a todos los ingenieros de la Facultad de Ingeniera Qumica por la enseanza que nos brindan.

Ejercicio 1%PROGRAMA PARA OBTENER LA DERIVADA clc; disp('==================================================') disp('-----------------***OBTENCION DE LA DERIVADA***---------------') disp('---Y GRAFICAS DE LAS FUNCIONES: INICIAL Y DERIVADA---) disp('==================================================') syms x y %syms define a X e Y como variables

A=input('ingrese funcion:') %pide la funcion a ser evaluada B=diff(A); %determina la derivada de la funcion inicial disp('la derivada de la funcion ser:') disp(B) %muestra en la pantalla el resultado de la derivada de la funcion inicial disp('A continuacion se mostrarn las graficas: ') disp('Primero de la Funcion inicial') disp('Luego de la funcion derivada') subplot(1,2,1), ezplot(y-A) %la funcion ezplot muestra el grafico panoramico de la funcion subplot(1,2,2), ezplot(y-B) %El ezplot muestra el titulo por defecto

SALIDA:=================================================== ----------------***OBTENCION DE LA DERIVADA***---------------------Y GRAFICAS DE LAS FUNCIONES: INICIAL Y DERIVADA---=================================================== ingrese funcion:3*x+x^2 A= 3*x+x^2 la derivada de la funcion ser: 3+2*x A continuacion se mostrarn las graficas: Primero de la Funcion inicial Luego de la funcion derivada

Ejercicio 2%PROGRAMA PARA HALLAR DERIVADAS PARCIALES y SU GRAFICA clc; disp('=================================================') disp('---------------------OBTENCION DE LA DERIVADA-------------------') disp('----------------------PARCIAL RESPECTO A X e Y----------------------') disp('=================================================') syms x y z %la funcion syms define como variables a X e Y

A=input('ingrese funcion:');

B=diff(A,x) C=diff(A,y) disp('la derivada respecto a x es:'),; disp(diff(A,x)) %calcula la derivada parcial respecto a X disp('la derivada respecto a y es:'); disp(diff(A,y)) %calcula la derivada parcial respecto a y subplot(1,2,1), ezplot(z-B) %la funcion ezplot muestra el grafico panoramico de la funcion subplot(1,2,2), ezplot(z-C) %El ezplot muestra el titulo por defecto

SALIDA:=================================================== ---------------------OBTENCION DE LA DERIVADA--------------------------------------------PARCIAL RESPECTO A X e Y-----------------------=================================================== ingrese funcion:x^2+y^3 la derivada respecto a x es: 2*x la derivada respecto a y es: 3*y^2

Ejercicio 3 %METODO de NEWTON RAPHSON mas rpido converge la funcin(llega a la solucion) clc; disp('===========================================') disp('-------------* METODO DE NEWTON RAPHSON *---------') disp('===========================================') nombre_f=input(' ingrese funcin asociada f(x)=','s'); x0=input(' ingrese valor inicial:'); fprintf('\n'); fprintf('iter aproximacin g(x) error\n');

i=1;e=1;delta=0.001; while e>=3E-12 & i>

Ejercicio 4 %METODO de VONN MISSES clc; disp('==========================================') disp('---------------* METODO DE VONN MISSES *-------------') disp('==========================================') nombre_f=input(' ingrese funcion asociada f(x)=','s'); x0=input(' ingrese valor inicial:'); fprintf('\n'); fprintf('iter aproximacion g(x) error\n'); i=1;e=1;delta=0.001;

x=x0; fx0=eval(nombre_f); x=x0-delta; df1=eval(nombre_f); x=x0+delta; df2=eval(nombre_f); dfx0=(df2-df1)/(2*delta); while e>=3E-12 & i