请求思路 和相关解法，高手赐教

zk1985

Fgp1=0.172.*(-335725854766920783./137438953472000.*cos(100.*pi.*t)-1197484068866141075195003436757375./2050028705056589735232949714944./(1-sin(11./46.*sin(100.*pi.*t)).^2).^(3./2).*sin(11./46.*sin(100.*pi.*t)).^2.*cos(11./46.*sin(100.*pi.*t)).^2.*cos(100.*pi.*t).^2-1197484068866141075195003436757375./2050028705056589735232949714944./(1-sin(11./46.*sin(100.*pi.*t)).^2).^(1./2).*cos(11./46.*sin(100.*pi.*t)).^2.*cos(100.*pi.*t).^2+1197484068866141075195003436757375./2050028705056589735232949714944./(1-sin(11./46.*sin(100.*pi.*t)).^2).^(1./2).*sin(11./46.*sin(100.*pi.*t)).^2.*cos(100.*pi.*t).^2+335725854766920783./137438953472000./(1-sin(11./46.*sin(100.*pi.*t)).^2).^(1./2).*sin(11./46.*sin(100.*pi.*t)).*cos(11./46.*sin(100.*pi.*t)).*sin(100.*pi.*t))

Fgr1=0.186*Fgp1

Fp1=a1*exp(-((t-b1)/c1).^2)+a2*exp(-((t-b2)/c2).^2)+a3*exp(-((t-b3)/c3).^2)+a4*exp(-((t-b4)/c4).^2)+a5*exp(-((t-b5)/c5).^2);（其中a1，a2等都为常数）

B（角度）=asin(sin(0.2391.*(sin(100.*pi.*t))))
上述参数都是关于t时间的函数。

现在知道 ax..+by..=tgB(Fp1+Fgp1+Fgr1+1.9992)
             cx..+dy..=-0.5(Fp1+Fgp1+1.6856)
(a  b  c  d都是可以确定的常数。而x..表示x的二阶倒数。y..相同）
现在要绘出y，x与时间t的函数曲线，直接dsolve命令求的的式子太常，无法得到x，y的函数，不知道用什么方法可以得到，或者其他的数值分析方法，请各位路过的指点。

xjzuo

建议将你的问题完整帖出来.
如果是微分方程(目前没看出来),可以用ode45等求解．