zhangxiaolu2015-CSDN博客

原创四阶龙格-库塔法求解常微分方程的初值问题-matlab通用程序

参考教材《数值分析》李乃成.梅立泉 clear clc format long m=input('请输入常微分方程的阶数m='); a=input('请输入x下限a='); b=input('请输入x上限b='); h=input('请输入步长h='); ym=input('令y(1,1)=y,y（2,1）=y’,y(3,1)=y’’...请输入ym=','s'); %输入

2015-01-15 18:46:36 17747

原创龙贝格积分-matlab通用程序

数值分析课程的大作业，教材《数值分析》李乃成.梅立泉 clear clc format long f=input('请输入原函数f=','s'); a=input('积分下限a='); b=input('积分上限b='); eps1=input('精度eps1='); T(1)=double((b-a)/2*(limit(sym(f),findsym(sym(f)),a)+limi

2015-01-15 18:42:23 17825 1

原创三次样条差值-matlab通用程序

数值分析课程的大作业，教材《数值分析》李乃成.梅立泉 clear x=input('请按照格式[x1,x2,x3...]格式输入y=f(x)函数已知点的横坐标xi='); %三次样条差值函数 y=input('请按照格式[y1,y2,y3...]格式输入y=f

2015-01-15 18:32:30 21387

原创共轭梯度法求解线性方程组-matlab通用程序

数值分析课程的大作业，教材《数值分析》李乃成.梅立泉 function x=Gongetidu2(A,b,x0,epsa) n=size(A,1); x=x0; r=b-A*x; d=r; for k=0:(n-1) alpha=(r'*r)/(d'*A*d); x=x+alpha*d; r2=b-A*x; if ((norm(r2)

2015-01-15 18:31:40 11034