blob: 40e49d44bcff180929024209dd4e764542d86264 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
clc;
clear;
m=input("Input value of m:")
h=input("Input value of h:")
function y=f(x)
y=x^m
endfunction
x1=1;
x2=x1+h;
fx1=f(x1);
fx2=fx1+m*(fx1^(m-1))*h;
if m==1 then
R=0;
else if m==2 then
R=2*(h^2)/factorial(2);
end
if m==3 then
R=(6*(x1)*(h^2)/factorial(2))+(6*(h^3)/factorial(3));
end
if m==4 then
R=(12*(x1^2)*(h^2)/factorial(2))+(24*(x1)*(h^3)/factorial(3))+(24*(h^4)/factorial(4));
end
end
disp(R,"Remainder:",fx2,"The value by first order approximation:")
disp(f(x2),"True Value at x2:")
|
