Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Community
Ask the community for help and clear up your study doubts
Discover the best universities in your country according to Docsity users
Free resources
Download our free guides on studying techniques, anxiety management strategies, and thesis advice from Docsity tutors
Mathemtathic Methods and Paleontologhy of imperal use in the divine comedy
Typology: Summaries
1 / 2
This page cannot be seen from the preview
Don't miss anything!
Xr=(Xa+Xb)/2;
% Programa para cálculo de Raíces de Ecuaciones no Lineales % Método de Bisección clear all; % Limpia la memoria clc; % Limpia la ventana de comandos % Ingresar la función, intervalo de evaluación y porcentaje de error. fprintf('\n Cálculo de la raíz de una ecuación por el método de Bisección \n\n'); Y=input('Dame la función : ','s'); Xa=input('Dame el intervalo inferior : '); Xb=input('Dame el intervalo superior : '); error=input('Dame el porciento del error : '); % Determinar si en el intervalo entre [Xa,Xb] está comprendida la raíz de la ecuación x=Xa; Ya=eval(Y); % f(Xa) x=Xb; Yb=eval(Y); % f(Xb) if (YaYb)> fprintf('\n\n No existe una raíz en el intervalo [Xa,Xb] \n\n'); fprintf('\Execute de nuevo el programa, por favor.\n\n'); % break end % Se realizan los cálculos para determinar la raíz en la siguiente sección. fprintf('\n\n') disp(' N Xa Xb Xr F(Xa) F(Xb) F(Xr) Error '); disp('|---|-----------------|------------------|-------------------|-------------------|----------- -------|--------------------|-----------'); Xant=0; % X anterior N=0; % Contador de Iteraciones while N< Xr=(Xa+Xb)/2; Xact=Xr; % X actual x=Xr; Yr=eval(Y); Ea=abs((Xact-Xant)/Xact)100; % Cálculo del Error ds_i=sprintf('%4d %17.8g %17.8g %17.8g %17.8g %17.8g %17.8g %10.5g', N, Xa, Xb, Xr, Ya, Yb, Yr, Ea); disp(ds_i); % Impresión de los Resultados if Ea<error fprintf('\n\n La Raíz Exacta es: %d',Xr) fprintf('\n\n Número de iteraciones: %d \n\n',N); break end % Se determina el nuevo intervalo de evaluación if (YaYr)< Xb=Xr; elseif (YaYr)== fprintf('\n\n La raíz exacta es: %17.8g',Xr) fprintf('\n\n Número de iteraciones: %d',N); break else Xa=Xr; end Xant=Xr; N=N+1; end
