Prova 3 Diego 2015/2 Nota 10

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

P3Frente.jpg
P3Verso.jpg
q1.sce
function y=f(x)
 y=x^3*log(x)
endfunction
dcentral=DC(f,12,0.001)
disp(dcentral)
d=derivative(f,12)
disp(d)
q10.sce
function F=f(y,t)
 F(1)=y(2)
 F(2)=-y(2)+cos(t)
endfunction
y3=[4;0]
[sol,t]=euler(f,3,y3,0.01,300)
plot(t,sol)
xgrid
q2.sce
function y=f(x)
 y=100*(sin(2*x^2)+cos(5*x))
endfunction
dprog=DP(f,2,0.0001)
disp(dprog)
d=derivative(f,2)
disp(d)
q3.sce
function y=f(x)
 y=93*sqrt(tan(x))
endfunction
ddupla=D2C(f,10,0.001)
disp(ddupla)
q4.sce
function y=f(x)
 y=exp(x-4)/sqrt(3*x^5)
endfunction
simps=simpson(f,3,7,1)
disp(simps)
q5.sce
function y=f(x)
 y=sin(cos(x^2))
endfunction
trap=trapezio(f,3,3.1,20)
disp(trap)
exata5=intg(3,3.1,f)
disp(exata5)
q6.sce
function y=f(x)
 y=sin(100*x)
endfunction
mediao=pmedio(f,-4,-3,100)
disp(mediao)
exata6=intg(-4,-3,f)
disp(exata6)
q7.sce
function y=f(x)
 y=log(1+sin(x))
endfunction
x=raizes_legendre(3)
gauss7=gauss(f,x,0,1)
disp(gauss7)
exata7=intg(0,1,f)
disp(exata7)
q9.sce
function F=f(y,t)
 F=-cos(y)
endfunction
y3=4
[y,t]=runge4(f,3,y3,0.1,10)
plot(t,y)
xgrid