EJEMPLO A CON BISECCION

clc
clear
x=[-1:0.01:2];
y1=1./(x+1);
y2=exp(x-1);
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off

function y=ejerAbis(x)
y=exp(x-1)-1/(x+1);
end
>> a=0, b=1, tol=5*10^(-6)
a=0
b=1
tol = 5.0000e-06
>> [c] = bisec('ejerAbis',a,b,tol)
1 0.000000000 1.000000000 0.500000000 -0.060136007
2 0.500000000 1.000000000 0.750000000 0.207372212
3 0.500000000 0.750000000 0.625000000 0.071904663
4 0.500000000 0.625000000 0.562500000 0.005648526
5 0.500000000 0.562500000 0.531250000 -0.027277215
6 0.531250000 0.562500000 0.546875000 -0.010825973
 7 0.546875000 0.562500000 0.554687500 -0.002592031
8 0.554687500 0.562500000 0.558593750 0.001527371
9 0.554687500 0.558593750 0.556640625 -0.000532543
10 0.556640625 0.558593750 0.557617188 0.000497360
11 0.556640625 0.557617188 0.557128906 -0.000017605
12 0.557128906 0.557617188 0.557373047 0.000239875
13 0.557128906 0.557373047 0.557250977 0.000111134
14 0.557128906 0.557250977 0.557189941 0.000046765
15 0.557128906 0.557189941 0.557159424 0.000014580
16 0.557128906 0.557159424 0.557144165 -0.000001512
17 0.557144165 0.557159424 0.557151794 0.000006534
c = 0.5572
EJEMPLO B CON BISECCION
clc
clear
x=[0:0.01:2];
y1=cos(x);
y2=x.^(1/2);
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off

function y=ejerBbis(x)
y=x^(1/2)-cos(x);
end
a=0.6, b=0.8, tol=5*10^(-6)

a = 0.6000
b = 0.8000
tol = 5.0000e-06
>> [c] = bisec('ejerBbis',a,b,tol)

1 0.600000000 0.800000000 0.700000000 0.071817839

2 0.600000000 0.700000000 0.650000000 0.010141976
3 0.600000000 0.650000000 0.625000000 -0.020393704
4 0.625000000 0.650000000 0.637500000 -0.005150267
5 0.637500000 0.650000000 0.643750000 0.002489686
6 0.637500000 0.643750000 0.640625000 -0.001331824
7 0.640625000 0.643750000 0.642187500 0.000578546
8 0.640625000 0.642187500 0.641406250 -0.000376735
9 0.641406250 0.642187500 0.641796875 0.000100882
10 0.641406250 0.641796875 0.641601562 -0.000137933
11 0.641601562 0.641796875 0.641699219 -0.000018527
12 0.641699219 0.641796875 0.641748047 0.000041177
13 0.641699219 0.641748047 0.641723633 0.000011325
14 0.641699219 0.641723633 0.641711426 -0.000003601
15 0.641711426 0.641723633 0.641717529 0.000003862
c = 0.6417
EJEMPLO C CON BISECCION
clc
clear
x=[0:0.01:3];
y1=(x.^3);
y2=x+2;
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off
function y=ejerCbis(x)
y=x^3-x-2
end
>> a=1.3, b=1.6, tol=5*10^(-6)
a = 1.3000
b = 1.6000
tol = 5.0000e-06
>> [c] = bisec('ejerCbis',a,b,tol)
y = -1.1030
y = 0.4960
y = -0.4014
1 1.300000000 1.600000000 1.450000000 -0.401375000
y = 0.021578
2 1.450000000 1.600000000 1.525000000 0.021578125
y = -0.1962
3 1.450000000 1.525000000 1.487500000 -0.196173828
y = -0.088886
4 1.487500000 1.525000000 1.506250000 -0.088886475
y = -0.034054
5 1.506250000 1.525000000 1.515625000 -0.034053802
y = -6.3381e-03
6 1.515625000 1.525000000 1.520312500 -0.006338055
y = 7.5949e-03
7 1.520312500 1.525000000 1.522656250 0.007594943
y = 6.2218e-04
8 1.520312500 1.522656250 1.521484375 0.000622176
y = -2.8595e-03
9 1.520312500 1.521484375 1.520898438 -0.002859506
y = -1.1191e-03
10 1.520898438 1.521484375 1.521191406 -0.001119057
y = -2.4854e-04
11 1.521191406 1.521484375 1.521337891 -0.000248539
y = 1.8679e-04
12 1.521337891 1.521484375 1.521411133 0.000186794
y = -3.0878e-05
13 1.521337891 1.521411133 1.521374512 -0.000030878
y = 7.7956e-05
14 1.521374512 1.521411133 1.521392822 0.000077956
y = 2.3539e-05
15 1.521374512 1.521392822 1.521383667 0.000023539
c = 1.5214
EJEMPLO D CON BISECCION
clc
clear
x=[2:0.01:3];
y1=(x.^4);
y2=(12*x)+8;
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off

function y=ejerDbis(x)
y=x^4-12*x-8
end
>> a=2.4, b=2.6, tol=5*10^(-6)

a = 2.4000
b = 2.6000
tol = 5.0000e-06
>> [c] = bisec('ejerDbis',a,b,tol)
y = -3.6224
y = 6.4976
y = 1.0625
1 2.400000000 2.600000000 2.500000000 1.062500000
y = -1.3700
2 2.400000000 2.500000000 2.450000000 -1.369993750
y = -0.1767
3 2.450000000 2.500000000 2.475000000 -0.176718359
y = 0.4371
4 2.475000000 2.500000000 2.487500000 0.437089868
y = 0.1287
5 2.475000000 2.487500000 2.481250000 0.128742799
y = -0.024348
6 2.475000000 2.481250000 2.478125000 -0.024347610
y = 0.052108
7 2.478125000 2.481250000 2.479687500 0.052107523
y = 0.013857
8 2.478125000 2.479687500 2.478906250 0.013857453
y = -5.2507e-03
9 2.478125000 2.478906250 2.478515625 -0.005250703
y = 4.3020e-03
10 2.478515625 2.478906250 2.478710938 0.004301969
y = -4.7472e-04
11 2.478515625 2.478710938 2.478613281 -0.000474719
y = 1.9135e-03
12 2.478613281 2.478710938 2.478662109 0.001913537
y = 7.1939e-04
13 2.478613281 2.478662109 2.478637695 0.000719387
y = 1.2233e-04
14 2.478613281 2.478637695 2.478625488 0.000122329
y = -1.7620e-04
15 2.478613281 2.478625488 2.478619385 -0.000176196
c = 2.4786
EJEMPLO E CON BISECCION
clc
clear
x=[2:0.01:3];
y1=4*((x-2).^2)-1;
y2=log(x);
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off
function y=ejerEbis(x)
y=log(x)+1-4*(x-2)^2;
end
>> a=2.5, b=3, tol=5*10^(-6)
a = 2.5000
b=3
tol = 5.0000e-06
>> [c] = bisec('ejerEbis',a,b,tol)
1 2.500000000 3.000000000 2.750000000 -0.238399088
2 2.500000000 2.750000000 2.625000000 0.402580896
3 2.625000000 2.750000000 2.687500000 0.097986393
4 2.687500000 2.750000000 2.718750000 -0.066234034
5 2.687500000 2.718750000 2.703125000 0.016869449
6 2.703125000 2.718750000 2.710937500 -0.024434000
7 2.703125000 2.710937500 2.707031250 -0.003720199
8 2.703125000 2.707031250 2.705078125 0.006590144
9 2.705078125 2.707031250 2.706054688 0.001438852
10 2.706054688 2.707031250 2.706542969 -0.001139703
11 2.706054688 2.706542969 2.706298828 0.000149817
12 2.706298828 2.706542969 2.706420898 -0.000494883
13 2.706298828 2.706420898 2.706359863 -0.000172518
14 2.706298828 2.706359863 2.706329346 -0.000011347
15 2.706298828 2.706329346 2.706314087 0.000069236
16 2.706314087 2.706329346 2.706321716 0.000028945
c = 2.7063
EJEMPLO F CON BISECCION
clc
clear
x=[1:0.01:3];
y1=(x.^2)
y2=(x+5)./2
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off

function y=ejerFbis(x)
y=2*x^2-x-5
end
>> a=1.5, b=2.5, tol=5*10^(-6)

a = 1.5000
b = 2.5000
tol = 5.0000e-06
>> [c] = bisec('ejerFbis',a,b,tol)

y = -2
y=5
y=1
1 1.500000000 2.500000000 2.000000000 1.000000000
y = -0.6250
2 1.500000000 2.000000000 1.750000000 -0.625000000
y = 0.1562
3 1.750000000 2.000000000 1.875000000 0.156250000
y = -0.2422
4 1.750000000 1.875000000 1.812500000 -0.242187500
y = -0.044922
5 1.812500000 1.875000000 1.843750000 -0.044921875
y = 0.055176
6 1.843750000 1.875000000 1.859375000 0.055175781
y = 5.0049e-03
7 1.843750000 1.859375000 1.851562500 0.005004883
y = -0.019989
8 1.843750000 1.851562500 1.847656250 -0.019989014
y = -7.4997e-03
9 1.847656250 1.851562500 1.849609375 -0.007499695
y = -1.2493e-03
 10 1.849609375 1.851562500 1.850585938 -0.001249313
y = 1.8773e-03
11 1.850585938 1.851562500 1.851074219 0.001877308
y = 3.1388e-04
12 1.850585938 1.851074219 1.850830078 0.000313878
y = -4.6775e-04
13 1.850585938 1.850830078 1.850708008 -0.000467747
y = -7.6942e-05
14 1.850708008 1.850830078 1.850769043 -0.000076942
y = 1.1847e-04
15 1.850769043 1.850830078 1.850799561 0.000118466
y = 2.0762e-05
16 1.850769043 1.850799561 1.850784302 0.000020762
y = -2.8090e-05
17 1.850769043 1.850784302 1.850776672 -0.000028090
c = 1.8508
EJEMPLO G CON BISECCION
clc
clear
x=[-2:0.01:1];
y1=(x.^3);
y2=exp(x)-3;
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off

function y=ejerGbis(x)
y=x^3+3-exp(x)
end
>> a=-1.5, b=-1, tol=5*10^(-6)
a = -1.5000
b = -1
tol = 5.0000e-06
>> [c] = bisec('ejerGbis',a,b,tol)
y = -0.5981
y = 1.6321
y = 0.7604
1 -1.500000000 -1.000000000 -1.250000000 0.760370203
y = 0.1476
2 -1.500000000 -1.250000000 -1.375000000 0.147551029
y = -0.2080
3 -1.500000000 -1.375000000 -1.437500000 -0.207979803
y = -0.025975
4 -1.437500000 -1.375000000 -1.406250000 -0.025974846
y = 0.061837
5 -1.406250000 -1.375000000 -1.390625000 0.061837002
y = 0.018195
6 -1.406250000 -1.390625000 -1.398437500 0.018194677
y = -3.8240e-03
7 -1.406250000 -1.398437500 -1.402343750 -0.003824013
y = 7.2018e-03
8 -1.402343750 -1.398437500 -1.400390625 0.007201828
y = 1.6930e-03
9 -1.402343750 -1.400390625 -1.401367188 0.001693034
y = -1.0645e-03
10 -1.402343750 -1.401367188 -1.401855469 -0.001064457
y = 3.1455e-04
11 -1.401855469 -1.401367188 -1.401611328 0.000314546
y = -3.7489e-04
12 -1.401855469 -1.401611328 -1.401733398 -0.000374891
y = -3.0156e-05
13 -1.401733398 -1.401611328 -1.401672363 -0.000030156
y = 1.4220e-04
14 -1.401672363 -1.401611328 -1.401641846 0.000142199
y = 5.6022e-05
15 -1.401672363 -1.401641846 -1.401657104 0.000056022
y = 1.2933e-05
16 -1.401672363 -1.401657104 -1.401664734 0.000012933
c = -1.4017
EJEMPLO H CON BISECCION
clc
clear
x=[0:0.01:1];
y1=(x.^3);
y2=(-3*x.^2)+1;
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off
function y=ejerHbis(x)
y=x^3+3*x^2-1
end
>> a=0.4, b=0.6, tol=5*10^(-6)

a = 0.4000
b = 0.6000
tol = 5.0000e-06
>> [c] = bisec('ejerHbis',a,b,tol)

y = -0.4560
y = 0.2960
y = -0.1250
1 0.400000000 0.600000000 0.500000000 -0.125000000
y = 0.073875
2 0.500000000 0.600000000 0.550000000 0.073875000
y = -0.028422
3 0.500000000 0.550000000 0.525000000 -0.028421875
y = 0.022006
4 0.525000000 0.550000000 0.537500000 0.022005859
y = -3.3875e-03
5 0.525000000 0.537500000 0.531250000 -0.003387451
y = 9.2643e-03
6 0.531250000 0.537500000 0.534375000 0.009264252
y = 2.9272e-03
7 0.531250000 0.534375000 0.532812500 0.002927174
y = -2.3294e-04
8 0.531250000 0.532812500 0.532031250 -0.000232944
y = 1.3464e-03
9 0.532031250 0.532812500 0.532421875 0.001346413
y = 5.5656e-04
10 0.532031250 0.532421875 0.532226562 0.000556559
y = 1.6176e-04
11 0.532031250 0.532226562 0.532128906 0.000161764
y = -3.5601e-05
12 0.532031250 0.532128906 0.532080078 -0.000035601
y = 6.3079e-05
13 0.532080078 0.532128906 0.532104492 0.000063079
y = 1.3738e-05
14 0.532080078 0.532104492 0.532092285 0.000013738
y = -1.0932e-05
15 0.532080078 0.532092285 0.532086182 -0.000010932
c = 0.5321
EJEMPLO I CON BISECCION
clc
clear
x=[0:0.01:1];
y1=4./(x+1);
y2=exp(x);
plot(x,y1,'r','linewidth',2)
hold on
plot(x,y2,'b','linewidth',2)
grid on
hold off

function y=ejerIbis(x)
y=(x+1)*exp(x)-4
end
>> a=0.6, b=1, tol=5*10^(-6)
a = 0.6000
b=1
tol = 5.0000e-06
>> [c] = bisec('ejerIbis',a,b,tol)
y = -1.0846
y = 1.4366
y = 5.9737e-03
1 0.600000000 1.000000000 0.800000000 0.005973671
y = -0.5766
 2 0.600000000 0.800000000 0.700000000 -0.576620397
y = -0.2952
3 0.700000000 0.800000000 0.750000000 -0.295249971
y = -0.1472
4 0.750000000 0.800000000 0.775000000 -0.147198974
y = -0.071263
5 0.775000000 0.800000000 0.787500000 -0.071263018
y = -0.032809
6 0.787500000 0.800000000 0.793750000 -0.032808552
y = -0.013459
7 0.793750000 0.800000000 0.796875000 -0.013458572
y = -3.7528e-03
8 0.796875000 0.800000000 0.798437500 -0.003752754
y = 1.1079e-03
9 0.798437500 0.800000000 0.799218750 0.001107881
y = -1.3231e-03
10 0.798437500 0.799218750 0.798828125 -0.001323081
y = -1.0776e-04
11 0.798828125 0.799218750 0.799023438 -0.000107761
y = 5.0002e-04
12 0.799023438 0.799218750 0.799121094 0.000500019
y = 1.9612e-04
13 0.799023438 0.799121094 0.799072266 0.000196119
y = 4.4176e-05
14 0.799023438 0.799072266 0.799047852 0.000044176
y = -3.1793e-05
15 0.799023438 0.799047852 0.799035645 -0.000031793
y = 6.1915e-06
16 0.799035645 0.799047852 0.799041748 0.000006192
c = 0.7990