sszhuanyi.m

来自「1．编写欧拉前差、后差、梯形公式。 2．编写二阶、三阶龙格库塔法通用程序。 」· M 代码 · 共 5 行

function [ft,fmt]=sszhuanyi(A,B,T)
syms s
ft=expm(A*T);                           %solve eAT

fmt=int(expm(A*(T-s))*B,s,0,T);         %solve int(0~T)eA(T-to)Bdto