matlab 微分方程求Dp/Dt=a+b*p最佳ab值（p,t已知多个值）

用dsolve('Dp=a+b*p','t')可以解得：p=-(a - c*exp(b*t))/b
然后求在最小二乘原则下的最佳的a、b值。
clear;
clc;
p=[1 4 5 6 10 11 ];
t=[86 86 49 130 114 103];
f=@(var) sum((p+(var(1) - var(3)*