See_Star的博客

1 题目(1) 编制用Crank-Nicolson格式求抛物方程数值解的通用程序。{∂u∂t−a∂2u∂x2=f(x,t)(0<x<1,0≤t≤ T)u(x,0)=φ(x)(0≤x≤1)u(0,t)=α(t),u(1,t)=β(t)(0<t≤T)\left\{\begin{array}{ll} \frac{\partial u}{\partial t}-a \frac{\partial^{2} u}{\partial x^{2}}=f(x, t) & (0<

1 题目用Romberg求积法计算积分∫−1111+100x2dx\int_{-1}^{1} \frac{1}{1+100 x^{2}} d x∫−11​1+100x21​dx的近似值，要求误差不超过0.5×10−70.5\times10^{-7}0.5×10−7。2 程序程序1：Romberg通用程序function[result] = romberg(a,b,F,epsilon) %% 初始赋值 e = inf; xk = [a,b];% T = []; %

一、复化梯形公式求积分题目：利用复化梯形公式计算I(f)=∫15sinxxdxI(f) = \int_1^5{\frac{sinx}{x}dx}I(f)=∫15​xsinx​dx的近似值，精确至7位有效数字。Matlab程序：clc; clear allsyms x %% 已知条件Fx(x) = sin(x)/x;a = 1;b = 5;epsilon = 5e-8;e = 1;T = []; % 存放Tf的值h = b-a;Xk1 = [a, b]; % 存放已选的 x

敲了一下午终于把代码敲完了，留念。文章目录题目：三次样条插值程序1：Gauss列主元消去法程序2：三次样条插值程序3：主函数及输入数据题目：三次样条插值程序1：Gauss列主元消去法function [X] = Gauss_elimination(A, Y) %% 参数初始化 matrix = [A, Y]; [n,m] = size(matrix); % 矩阵大小 l = zeros(n,n); % 比例系数矩阵 X =

文章目录题目：牛顿迭代法程序1：牛顿迭代法通用函数程序2：求最大Delta程序3：观察结果题目：牛顿迭代法程序1：牛顿迭代法通用函数function [x] = newton(x0,epsilon,f,print_flag)digits(10) % 控制牛顿迭代法的运算精度，精度太大迭代速度过慢syms xf(x) = f;df(x) = diff(f);count = 0;e = 1; while abs(e) > epsilon x1 = vpa(x0