s1=newton(x^2-10,2*x,3,1e-14) {??? Undefined function or variable 'x'. } s1=newton('x^2-10','2*x',3,1e-14) s1 = 3.0000 3.1667 3.1623 3.1623 3.1623 format long g s1 s1 = 3 3.16666666666667 3.16228070175439 3.16227766016984 3.16227766016838 s2=newton('x^3+4*x^2-10','3*x^2+8*x',1,1e-14) s2 = 1 1.45454545454545 1.36890040106952 1.36523660020212 1.36523001343537 1.3652300134141 s3=newton('exp(x)-tan(x)','exp(x)-(sec(x))^2',-4,1e-14) s3 = -4 -3.49353222947037 -3.13351894660936 -3.09641459850047 -3.09641230491352 -3.09641230491364 s4=newton2('exp(x)-tan(x)','exp(x)-(sec(x))^2',-4,1e-14) s4 = -4 -3.49353222947037 -3.13351894660936 -3.09641459850047 -3.09641230491352 -3.09641230491364 -3.09641230491364 fzero(@(x) x^2-10, 3) ans = 3.16227766016838 fzero(@(x) x^3 + 4*x^2 -10, 1) ans = 1.3652300134141 H=[.919 .175 0; .06 .825 .049; .021 0 .951]; x=[100000; 70000; 30000] x = 100000 70000 30000 H*x ans = 104150 65220 30630 H*ans ans = 107127.35 61556.37 31316.28 H*ans ans = 109222.3994 58746.14397 32031.45663 H*ans ans = 110655.96024335 56588.45411412 32755.58564253 transmat1 x = 109656.068283714 49317.3395585203 41026.5921577655 transmat1 ans = Column 1 111381.428643817 49833.8949287744 38784.6764274081 Column 2 111080.464536204 49696.298179811 39223.2372839853 Column 3 110779.799090238 49586.2124974315 39633.9884123303 Column 4 110484.222550979 49497.4786879995 40018.2987610211 Column 5 110197.05929475 49425.3699099484 40377.5707953016 Column 6 109920.537226116 49366.2547023622 40713.2080715216 transmat1 ans = Column 1 105761.057260542 48944.0760071557 45294.8667323022 Column 2 105759.62492369 48943.9746114188 45296.4004648908 Column 3 105758.29086187 48943.8801726215 45297.8289655087 Column 4 105757.048332267 48943.7922134349 45299.159454298 Column 5 105755.891054705 48943.7102892804 45300.398656015 Column 6 105754.813179898 48943.6339860833 45301.5528340191