function r=myfun1
f=@(x,y)(exp(2.*x.^2+3.*y.^2));
r=integral2(f,0.5,1,0,1);
end
function r=myfun2
syms x;
f=@(x)(1./(sqrt(x)+1));
r=quad(f,1,4);
end
function C=myfun3
A=[1.15 1.12 1.12 1.06 1.16];
C=[1.15 1.12 1.12 1.06 1.16];
for i=2:20
B=power(A,i);
C=[C;B];
end
C
end
function [x0,y0]=myfu