Matlab - Trapecio compuesto en matlab

%sin modificar mucho tu código

clc

clear

f='x^2+3*x+4';

a=2;

b=-2;

n=8;

% f funcion

% a,b intevalo

% n numero partes

disp('Funcion: ');

f

disp(strcat('De [a: ',num2str(a),'Hacia b : ',num2str(b),']'));

f=inline(f);

h=(b-a)/n;

aprox=f(a)+f(b);

for i=1:n-1

x=a+i*h;

aprox=aprox+2*f(x);

end

aprox=(h/2)*aprox;a=0;

disp(aprox)

clc;

clear all

fun=input('Introduzca la funcion: f(x)=','s');

i=input('Introduzca el número de iteraciones : ');

a=input('Introduzca a : ');

b=input('Introduzca b : ');

e=input('Introduzca la tolerancia : ');

n=1;

f=inline(fun);

fprintf('\n %6s %7s %8s %10s %8s %8s %8s \n ','a','b','c','f(a)','f(b)','f(c)','|c-a|');

while n<=i

    c=(a+b)/2;

     fprintf('\n %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f \n',a,b,c,f(a),f(b),f(c),abs(c-a));

    if f(c)==0

        fprintf('La raiz exacta es: \n %8.4f \n',c);

        break;

    elseif (b-a)/(2^n)<e

        fprintf('La raiz aproximada es: \n %8.4f \n',c);

        break;

    elseif f(a)*f(c)<0

        b=c;

    else

        a=c;

    end

n=n+1;

end

Introduzca la funcion: f(x)='x^3+2*x^2+10*x-20'

Introduzca el número de iteraciones : 30

Introduzca a : 1

Introduzca b : 2

Introduzca la tolerancia : 0.001

      a       b        c       f(a)     f(b)     f(c)    |c-a|

   1.0000   2.0000   1.5000 120.0000  94.0000  51.0000  43.0000

  50.0000  42.0000 120.0000  94.0000  50.0000  43.0000  49.0000

  48.0000  42.0000 120.0000  45.0000  50.0000  48.0000 120.0000

  94.0000  51.0000  43.0000  50.0000  42.0000 120.0000  94.0000

  50.0000  43.0000  49.0000  48.0000  42.0000 120.0000  45.0000

  50.0000  48.0000 120.0000  94.0000  51.0000  43.0000  50.0000

  42.0000 120.0000  94.0000  50.0000  43.0000  49.0000  48.0000

  42.0000 120.0000  45.0000  50.0000  48.0000   0.5000 Error using  *

Inner matrix dimensions must agree.

clc;

clear all

fun=input('Introduzca la funcion: f(x)=','s');

i=input('Introduzca el número de iteraciones : ');

a=input('Introduzca a : ');

b=input('Introduzca b : ');

e=input('Introduzca la tolerancia : ');

n=1;

f=inline(fun);

fprintf('\n %6s %7s %8s %10s %8s %8s %8s \n ','a','b','c','f(a)','f(b)','f(c)','|c-a|');

while n<=i

    c=(a+b)/2;

     fprintf('\n %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f \n',a,b,c,f(a),f(b),f(c),abs(c-a));

    if f(c)==0

        fprintf('La raiz exacta es: \n %8.4f \n',c);

        break;

    elseif (b-a)/(2^n)<e

        fprintf('La raiz aproximada es: \n %8.4f \n',c);

        break;

    elseif f(a).*f(c)<0

        b=c;

    else

        a=c;

    end

n=n+1;

end

* el primer punto es con la funcion e^-x-x   con      x0=0

con 5 iteracciones, sin tolerancia

* el segundo punto  es con la funcion e^-x-x   con      x0=0

con 5 iteraccciones   y tolerancia de 10^-8

%Método de newton raphson

clear, clc

disp('METODO NEWTON RAPHSON');

disp('-------------------');

cf = input('ingrese funcion a evaluar: ');% lemos fucion que se dese evaluar

syms x % variable simbolica 

f = inline (cf);% transformo esa cadena en una funcion inline le asigno variable f

derivada = diff(cf, x);% diff derivada requiere de la cadena cf y de la variable que se va a deriva x

df = inline (derivada);% variable df que al final usare

tol = input ('ingrese tolerancia: ');% ingreso de la tolerancia

error = 50;% establecer un error grande para exceda al bucle while

x = input ('ingrese un valor inicial: ');% valor variable x

n=0;% inicio conteo de iteraciones igual a 0

disp ('    n     xi     error')% encabezado al ejecutar el programa, xo x1 aproximados a la raiz y error que se comete en cada x

while(error>tol)% mientras el errosr sea mayor que la tolerancia continuara iterando hasta llegar que el error sea menor o igual a la tolerancia

    fprintf('\t%i\t%3.5f\t%f\n', n, x, error);

    n = n+1;% aumentar numero de iteraciones

    x = x - f(x)/df(x);% calcular nuevo valor de x

    error = abs (f(x));% error se actualiza como valor absoluto de evalluar x de fx

end