r语言计算微积分例题及代码.docVIP

计算向量、矩阵差分 x=1:5 x diff(x) [1] 1 1 1 1 y=1:12 z=matrix(y,3,4,T) z diff(z) [,1] [,2] [,3] [,4] [1,] 4 4 4 4 [2,] 4 4 4 4 diff(x,diff=2) [1] 0 0 0 #求一元函数导数 f=expression((5*x)/(1+x^2)) D(f,x) f=expression(2*x^2) D(f,x) #求二元函数偏导数及梯度 f=deriv(z~x^2+y^2+x*y,c(x,y),func=T, hessian=F) f f(1,2) 例1梯形积分 trape=function(fun,a,b,tol=1e-6){ N=1;h=b-a T=h/2*(fun(a)+fun(b)) repeat{ h=h/2;x=a+(2*1:N-1)*h I=T/2+h*sum(fun(x)) if(abs(I-T)tol) break N=2*N;T=I} I} source(trape.R) f=function(x) exp(-x^2) trape(f,-1,1) 高精度积分 f=function(x) exp(-x^2) integrate(f,-1,1) 例2简单二重积分 library(