(X, Y) Meshgrid (1:0.5:10,1:20) Z Sin (X) + Cos (Y) Surf (X, Y, Z)
(X, Y) Meshgrid (1:0.5:10,1:20) Z Sin (X) + Cos (Y) Surf (X, Y, Z)
(X, Y) Meshgrid (1:0.5:10,1:20) Z Sin (X) + Cos (Y) Surf (X, Y, Z)
5;
[X,Y]=meshgrid(x,y);
2
F=X.^2+Y.^2;
surf(X,Y,F);
xlabel('x');
ylabel('y');
[X,Y] = meshgrid(1:0.5:10,1:20);
Z = sin(X) + cos(Y);
surf(X,Y,Z)
t = 0:pi/10:2*pi;
figure
[X,Y,Z] = cylinder(sin(t)+cos(t));
surf(X,Y,Z)
axis square
x=-1:.1:1;y=0:.1:1.5;
[X,Y]=meshgrid(x,y);
2
F=(cos(X)+sin(Y));
surf(X,Y,F);
xlabel('x');
ylabel('y');
x=-1:.1:1;y=0:.1:1.5;
[X,Y]=meshgrid(x,y);
2
F=Y.^2+5*cos(pi*X);
surf(X,Y,F);
xlabel('x');
ylabel('y');
x=4:4.2:5.2;
y=log(x);
trapz(x,y)
clear all
clc
syms x y
f=[x.^2-y.^2 2*x*y];
disp('Along the curve y=x.^2')
a=subs(f,y,x.^2);
b=diff(x.^2,x);
c=b*a(2);
d=int(a(1),x,0,1);
e=int(c,x,0,1);
u=d+e
disp('Along the curve y=sqrt(x)')
p=subs(f,y,sqrt(x));
q=diff(sqrt(x),x);
r=q*p(2);
s=int(p(1),x,1,0);
t=int(r,x,1,0);
v=s+t
I=u+v
x=-2:0.5:2;
y2=sqrt(x);
y1=x.^2;
plot(x,y1,'r', x,y2,'g');
grid on
f=@(x,y) (3*x^2+1);
[x,y]=ode23(f,[1:0.5:2],2)
[x,y]=ode45(f,[1:0.5:2],2)
f=@(x,y) (x+y+x*y);
[x,y]=ode23(f,[0:0.025:1],1)
[x,y]=ode45(f,[0:0.025:1],1)
f=@(x,y) (x.^2+y.^2);
[x,y]=ode23(f,[0:0.1:0],1)
[x,y]=ode45(f,[0:0.1:0],1)
f=@(x,y) (2*exp(x)-3*y);
[x,y]=ode23(f,[0:0.25:0],0)
[x,y]=ode45(f,[0:0.25:0],0)
9–a
%Bisection Method%
f=@(x) x^3-x-1;
display('Equation is x^3-x-1= 0')
i=1;
while(i)
xl=input('Enter lower value:');
xu=input('Enter upper value: ');
e=input('Enter accuracy: ');
if f(xl)*f(xu)<0
i=0;
else
warning('Enter proper range');
end
end
if f(xl)<0
xn=xl;
xp=xu;
else
xn=xu;
xp=xl;
end
while (abs(xn-xp)>e)
xm=(xn+xp)/2;
if f(xm)<0
xn=xm;
else
xp=xm;
end
end
Root=xm
%Regula Falsi%
f=@(x) x^3-x-1;
display('Equation is x^3-x-1=0')
i=1;
while(i)
xl=input('Enter lower value:');
xu=input('Enter upper value: ');
e=input('Enter accuracy: ');
if f(xl)*f(xu)<0
i=0;
else
warning('Enter proper range');
end
end
if f(xl)<0
xn=xl;
xp=xu;
else
xn=xu;
xp=xl;
end
xm=xl;
while (abs(f(xm))>e)
xm=(xn*f(xp)-xp*f(xn))/(f(xp)-f(xn));
if f(xm)<0
xn=xm;
else
xp=xm;
end
end
Root=xm
syms x
f= x^3-x-1;
fdash=diff(f);
disp('The equation is: '),disp(f);
disp('The derivative of equation is: '),disp(fdash);
y=inline(f);
dy=inline(fdash);
x0=input('Enter approximate root: ');
e=input('Enter the accuracy: ');
while abs(feval(y,x0))>e
h=-feval(y,x0)/feval(dy,x0);
x0=x0+h;
end
root=x0
9–b
%Bisection Method%
f=@(x) 4*sin(x)-exp(x);
display(‘Equation is 4*sin(x)-exp(x)=0’)
i=1;
while(i)
xl=input('Enter lower value:');
xu=input('Enter upper value: ');
e=input('Enter accuracy: ');
if f(xl)*f(xu)<0
i=0;
else
warning('Enter proper range');
end
end
if f(xl)<0
xn=xl;
xp=xu;
else
xn=xu;
xp=xl;
end
while (abs(xn-xp)>e)
xm=(xn+xp)/2;
if f(xm)<0
xn=xm;
else
xp=xm;
end
end
Root=xm
%Regula Falsi%
f=@(x) 4*sin(x)-exp(x);
display('Equation is 4*sin(x)-exp(x)=0’)
i=1;
while(i)
xl=input('Enter lower value:');
xu=input('Enter upper value: ');
e=input('Enter accuracy: ');
if f(xl)*f(xu)<0
i=0;
else
warning('Enter proper range');
end
end
if f(xl)<0
xn=xl;
xp=xu;
else
xn=xu;
xp=xl;
end
xm=xl;
while (abs(f(xm))>e)
xm=(xn*f(xp)-xp*f(xn))/(f(xp)-f(xn));
if f(xm)<0
xn=xm;
else
xp=xm;
end
end
Root=xm
syms x
f=4*sin(x)-exp(x);
fdash=diff(f);
disp('The equation is: '),disp(f);
disp('The derivative of equation is: '),disp(fdash);
y=inline(f);
dy=inline(fdash);
x0=input('Enter approximate root: ');
e=input('Enter the accuracy: ');
while abs(feval(y,x0))>e
h=-feval(y,x0)/feval(dy,x0);
x0=x0+h;
end
root=x0
9–c
%Bisection Method%
f=@(x) 2*x-cos(x)-3;
display('Equation is 2*x-cos(x)-3=0’)
i=1;
while(i)
xl=input('Enter lower value:');
xu=input('Enter upper value: ');
e=input('Enter accuracy: ');
if f(xl)*f(xu)<0
i=0;
else
warning('Enter proper range');
end
end
if f(xl)<0
xn=xl;
xp=xu;
else
xn=xu;
xp=xl;
end
while (abs(xn-xp)>e)
xm=(xn+xp)/2;
if f(xm)<0
xn=xm;
else
xp=xm;
end
end
Root=xm
%Regula Falsi%
f=@(x)2*x-cos(x)-3;
display('Equation is 2*x-cos(x)-3=0')
i=1;
while(i)
xl=input('Enter lower value:');
xu=input('Enter upper value: ');
e=input('Enter accuracy: ');
if f(xl)*f(xu)<0
i=0;
else
warning('Enter proper range');
end
end
if f(xl)<0
xn=xl;
xp=xu;
else
xn=xu;
xp=xl;
end
xm=xl;
while (abs(f(xm))>e)
xm=(xn*f(xp)-xp*f(xn))/(f(xp)-f(xn));
if f(xm)<0
xn=xm;
else
xp=xm;
end
end
Root=xm
syms x
f=2*x-cos(x)-3;
fdash=diff(f);
disp('The equation is: '),disp(f);
disp('The derivative of equation is: '),disp(fdash);
y=inline(f);
dy=inline(fdash);
x0=input('Enter approximate root: ');
e=input('Enter the accuracy: ');
while abs(feval(y,x0))>e
h=-feval(y,x0)/feval(dy,x0);
x0=x0+h;
end
root=x0