微分方程的求解程序问题

猴王

这是一个微分方程的求解程序，其方程是这样的mv''+cv'+kv=f,我想求出这个方程的v''，请指教是不是求解出了什么问题：
function xdot=danziyoudukuangjia(t,x)
F=20*sin(25*t);
  M=31.2;
  C=0.01;
  K=67.51;
xdot=[x(2);(1/M)*(F-C*x(2)-K*x(1))];
执行文件：
clc;clear all;
close all;
t0=0;
   tf=60;
   x0=[0,0]';

   [t,x]=ode23('danziyoudukuangjia',[0:0.005:60],x0);
   d=x(:,1);
   v=x(:,2);
   
t=0:0.005:60;

F=20*sin(25*t);                                
M=31.2;
C=0.01;
K=67.51;
s1=(1/M)*(F-C*v.'-K*d.');

[ 本帖最后由 eight 于 2007-6-25 10:18 编辑 ]

xjzuo

程序本身没有什么问题.

猴王

这个结果与我师兄的结果之间还有差距，我想知道如果我通过对v(:,1)进行求导，这样通过对速度曲线求导从而得到加速度曲线的方式是否可行，应该怎样通过程序实现，请指教，请帮我写出程序的框架好吗，还有就是程序我知道没有问题，求解过程您看有没有问题

xjzuo

这个自己看看公式不就明白了?