close all
clear all
clc
x=-2:0.1:2;w=1
A=1;
B=1;
t0=-2.*pi;dt=0.01;tk=2.*pi;
w=1;
t=t0:dt:tk;


[W,CZAS]=meshgrid(w,t)
 y= A.*((-5).*x.^2+x.^4+4).*B.*cos(w.*CZAS);
figure(1)
 subplot(3,1,1)
 [X,CZAS]=meshgrid(x,t);
 mesh(CZAS,X,y)
function[CZAS,X,y]=Krzysztof_Kubacki(el1,el2,el3,el4,el5,el6,el7)
A=el1;
B=el2;
w=el3;
x=el4;
t0=el5;
dt=el6;
tk=el7;
t=el5:dt:el7
 end